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Abstract 



In this thesis, we present theoretical and statistical techniques broadly related 
to systems of dynamically-interacting particles. We apply these techniques to 
observations of dense star clusters in order to study gravitational interactions 
between stars. These include both long- and short-range interactions, as well as 
encounters leading to direct collisions and mergers. The latter have long been 
suspected to be an important formation channel for several curious types of 
stars whose origins are unknown. The former drive the structural evolution of 
star clusters and, by leading to their eventual dissolution and the subsequent 
dispersal of their stars throughout the Milky Way Galaxy, have played an 
important role in shaping its history. Within the last few decades, theoretical 
work has painted a comprehensive picture for the evolution of star clusters. 
And yet, we are still lacking direct observational confirmation that many of 
the processes thought to be driving this evolution are actually occuring. The 
results presented in this thesis have connected several of these processes to real 
observations of star clusters, in many cases for the first time. This has allowed 
us to directly link the observed properties of several stellar populations to the 
physical processes responsible for their origins. 

We present a new method of quantifying the frequency of encounters 
involving single, binary and triple stars using an adaptation of the classical 
mean free path approximation. With this technique, we have shown that dy- 
namical encounters involving triple stars occur commonly in star clusters, and 
that they are likely to be an important dynamical channel for stellar mergers 
to occur. This is a new result that has important implications for the origins 
of several peculiar types of stars (and binary stars), in particular blue strag- 
glers. We further present several new statistical techniques that are broadly 
applicable to systems of dynamically-interacting particles composed of sev- 
eral different types of populations. These are applied to observations of star 
clusters in order to obtain quantitative constraints for the degree to which dy- 
namical interactions affect the relative sizes and spatial distributions of their 
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different stellar populations. To this end, we perform an extensive analysis 
of a large sample of colour- magnitude diagrams taken from the ACS Survey 
for Globular Clusters. The results of this analysis can be summarized as fol- 
lows: (1) We have compiled a homogeneous catalogue of stellar populations, 
including main-sequence, main-sequence turn-off, red giant branch, horizontal 
branch and blue straggler stars. (2) With this catalogue, we have quantified 
the effects of the cluster dynamics in determining the relative sizes and spatial 
distributions of these stellar populations. (3) These results are particularly in- 
teresting for blue stragglers since they provide compelling evidence that they 
are descended from binary stars. (4) Our analysis of the main-sequence pop- 
ulations is consistent with a remarkably universal initial stellar mass function 
in old massive star clusters in the Milky Way. This is a new result with impor- 
tant implications for our understanding of star formation in the early Universe 
and, more generally, the history of our Galaxy. Finally, we describe how the 
techniques presented in this thesis are ideally suited for application to a num- 
ber of other outstanding puzzles of modern astrophysics, including chemical 
reactions in the interstellar medium and mergers between galaxies in galaxy 
clusters and groups. 
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Chapter 



1 



Introduction 



The vast majority of the stars in our Universe are found in galaxies, each of 
which is typically populated by billions of members. The hierarchy of stellar 
communities within galaxies can be understood if an analogy is made between 
stars and people. In this context, parallels can be drawn between galaxies 
and countries, star clusters and cities, as well as multiple star systems and 



families. It is now known that at 
in star clusters (e.g. 



Kroupa , 



2002 



east half of the stars in galaxies were born 



Kroupa fc Boily 



20021 ). spherical conglom- 



erates composed of anywhere from a few hundred to a few million members 
all orbiting their common center of mass (see Figure II. ip . Within these stellar 
nurseries, around 50% of all stars are thought to spend the majority of their 



lives in close pro ximity to one or more o 



star systems ( e 



Kroupa 



2001 



g. 



Durisen fc Sterzik 



Sollima et al. 



1994 



their siblings by 



Bate fc Bonnell 



form i ng multiple 



1997 



Bate 



1997 



20101 ). Most of these are binary star systems, 
which consist of two stars orbiting their common centre of mass under the in- 
fluence of their mutual gravitation attraction. Although less common, stable 



configurations of triple, quadruple and quintuple star systems also exist, along 
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with even highe r order multiples (e.g 



Latham 



2005 



2008 



Tokovininl . 



Eggleton fc Tokovinin 



20081 ). Despite the fact that stars typically spend at most 
a small fraction of their considerably long lives in even relative isolation, our 
understanding of the mechanisms via which the presence of their siblings can 



affect their evolution is far from complete. 



Figure 1.1 An image taken with the Hubble Space Telescope of NGC 7078 
(M15), a globular cluster in the Milky Way galaxy. Image credit: NASA, 
ESA, and The Hubble Heritage Team. 



Gravity mediates a wide variety of interactions between stars. Within 
the dense cores of star clusters, for instance, strong gravitational encounters 



often occur in which the distance of c 
radii of the stars (e.g . 



1975 



Hut fc Bahcall 



Spitzer fc Thuan 



osest approach is comparab 



1972 



Spitzer &: Shull 



1975 



e to the 



Heggie 



19831 ) . These encounters often lead to co mplex dances 



that can result in stars being ejected from their host clusters (e.g. 



Henon 



1969 
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Hut & Diorgovski. 


1992 


Kroupa & Boilvl. 


2002 


De Marchi. Paresce & Portegies Zwart 


2010 


), or even stellar collisions (e.g. 


Leonard 


■ - 


L989; 


Sills et al.. 


2001 


). Binary 



stars are also at the mercy of gravity and, if their orbital separation becomes 
sufficiently small, this can le ad to the trans fer of mass from the surface of 



one compo nent to the other (IMcCreal. 



cence (e.g. 



19641) . or even their complete coales- 



Andronov. Pinsonneault fc Terndrup 



20061 ) . Not only are gravi- 



tational interactions important for stellar mergers, they also drive the evo- 
lution of star clusters. Gravity slowly acts to dissolve clusters, continuously 
ejecting: stars so that they m ay join the rest of the Galactic population (e.g. 



Portegies Zwart et al. 



20011 ) . It follows that a complete understanding of the 
dominant physical processes operating in clusters, including how their stars 
are formed and the dynamical interactions that cause their stars to escape, is 
a key ingredient in piecing together the history of our Galaxy. 

Mergers are a particularly interesting and important outcome of stellar 
interactions. Observations have revealed numerous examples of curious stars 
and mysterio us astrophysical processes that are thought to be related to stellar 



merg e rs (e.g 



2010 



Sandagei. 



1953 



Miroshnichenko et al 



Webbink 



2007 



1984 



Paczvnski 



Farrell et al 



1986 



Troja. Rosswog fc Gehrelsl . 



20091 ). Many of these ob- 



jects are commonly found in dense star clusters, hinting at the complexity of 
their dynamical evolution and the potential importance of dynamics for stellar 
mergers. Although their origins remain unknown, some of these objects have 
served as important tools for furthering our understanding of the Universe. 
For instance, Type la supernovae have provide d the most robust const raints 
to date for the rate of Universal expansion (e.g. iPerlmutter et all Il999l ). De- 
spite their considerable astrophysical significance, we are far from completely 
understanding the physical mechanisms that cause mergers to occur. 
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Stellar mergers are but one example of this ubiquitous physical process, 
which occurs between objects on many scales. From colliding galaxies to fusing 
atoms, mergers occur on spatial and temporal scales ranging by many orders of 
magnitude. Familiar principles such as conservation of energy and momentum 
often apply, however, regardless of scale. This is certainly the case for the 
various forms of dynamical interactions that lead to mergers. Examples include 
chemical reactions in the interstellar medium, strong gravitational interactions 
involving massive black holes and even encounters between galaxies in galaxy 
clusters. It follows that the development of tools designed to further our 
understanding of stellar dynamics and mergers also have applications for a 
number of other physical sub-disciplines. 

In this thesis, we present statistical and theoretical techniques related 
to the physics of mergers. These are applied to observations of dense star 
clusters in order to study the mergers of stars. In this chapter, we will outline 
our motivation for conducting this research. In Section 11.11 we briefly review 
our current understanding of stellar evolution theory, and discuss the observed 
properties of star clusters and their stellar populations, including the various 
types of stellar exotica. In Section 11.21 we provide a brief description of the 
dominant physical processes that drive the dynamical evolution of star clusters. 
All of these issues are connected in Section 11.31 where we describe how they 
have motivated the development of the techniques that will be presented in 
this thesis, along with their application to observations of star clusters. 
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1.1 Stellar Populations in Star Clusters 
1.1.1 Single Star Evolution 

By serving as sites for star formation, star clusters have played a crucial role 
in shaping the present-day features of our Galaxy. And yet, this is but one 
example of their astrophysical significance. Star clusters are also ideal labo- 
ratories for learning about stellar evolution. Historically, their importance in 
this regard has stemmed from the often-adopted assumption that all of the 
stars in a given cluster were born from the same gas cloud at more or less 
the same time. If true, star clusters offer a large sample of stars with a wide 
range of masses, but a very narrow range in both age and chemical composi- 
tion. Therefore, observations of star clusters offer robust constraints for stellar 
evolution theories by providing direct tests for their predictions for a diverse 
spectrum of stellar masses of known age and composition. 

Colour-magnitude diagrams (CMDs) are one of the m ost important 



tools available to astronomers for studying stellar evolution (IHertzsprung 



19091 ). In these diagrams, the observed surface temperature is plotted against 
luminosity (or brightness) for every star in the cluster. The latter quantity is 
plotted backwards, and provides a proxy for colour - hot stars are blue and 
cool stars are red. This creates a characteristic appearance that is observed in 
all cluster CMDs. An example is shown in Figure 11.21 for the GC NGC 6205. 
This characteristic shape can be understood by considering two principles of 
stellar evolution. The first is called the Vogt-Russell Theorem and states that 
the mass, composition and age of a star uniquely determine its radi us, luminos 



ity and intern al structure, in addition to its subsequent evolution (IVogt 



1925 



Russell 



19251 ). The second is simply that the rate at which a star evolves, and 
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in so doing changes its brigh t ness a nd colour, is inversely proportional to some 



power of its mass (e.g. 



Ib en 



19911 ). It follows that, for a cluster composed of 
a large number of stars with a distribution of masses but similar compositions 
(at birth) and ages, the most massive stars are also the most evolved. There- 
fore, every evolutionary phase in the life of a star is typically represented in at 
least moderately old clusters, and this causes the distribution of stars in the 
CMD to adhere to the characteristic shape shown in Figure 11.21 Our under- 
standing of single star evolution is sufficiently complete that it is now known 
how this shape chang function of composition and age. In general, the 

majority of the features characteristic of CMDs are well understood, although 
exceptions do exist. Before getting to these in Section H. 1.21 we will review the 
life cycle of a typical star in our Galaxy and connect each evolutionary phase 
to its corresponding location in the CMD. 

Most of the star clusters in our Galaxy were born from gas clouds 
composed primarily of hydrogen and, t o a lesser extent, he l ium ( as well as 



Lada 



1985 



Pringle 



19891 ). Gravity 



heavier elements in trace amounts) (e.g. 
contracted the gas into clumps, which grew increasingly hot and dense as they 
accreted material from the surrounding medium. When the central tempera- 
ture of a clump becomes sufficiently high (on the order of 10 K), hydrogen is 
ignited at its centre. A star is born. This marks the beginning of the main- 
sequence (MS) phase, during which time hydrogen is fused into helium in the 
stellar core. 

Let us consider the life of a typical low-mass star in the Milky Way, 
from beginning to end. Most of the stars in our G alaxy have low masses 



roughly spanning the range 0.08 < m < 2.0 M (e.g. 



Kroupa 



20021). There 



are two reasons for this. First, the least massive stars are also the longest 



6 



Ph.D. Thesis — N. Leigh 



McMaster - Physics & Astronomy 




0.5 1 
Colour 



Figure 1.2 Colour- magnitude diagram for the Milky Way glob ular cluster NGC 
6205. The data used to create this plot were taken from ISarajedini et al. 
(120071 ). As described in the text, the x-axis can be read as temperature (or, 
equivalently, colour), which increases from right to left (red to blue). On the 
y-axis, brightness increases from bottom to top. Labels for blue straggler 
(BS), red giant branch (RGB), horizontal branch (HB), extended horizontal 
branch (EHB), asymptotic giant branch (AGB), Hertzsprung gap (HG), main- 
sequence turn-off (MSTO) and main-sequence (MS) stars are shown. Stars 
with large photometric errors have been omitted from this plot. 



lived. Second, ISalpeterl (jl955l ) showed that the initial distribution of stellar 
masses for stars in the solar-neighbourhood with masses in the range 0.4 — 10 
M Q can be described by a power-law form with an index a ~ 2.35. Today, 
we know that the distribution flattens below ~ 0.5 M Q , however the primary 
conclusion is the same: most of the stars in our Galaxy have masses lower 



than that of our Sun. What's more, most of these are MS stars since, for stars 
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of very low-mass, the time-scale for this evolutionary phase exceeds the age of 
the Universe. 

Nuclear reactions in the stellar core release heat, and are the pri- 



mary sc^rrceo 



[e.g. 



Clayton 



ener gy generation within stars, driving their high luminosities 



19681 ). The transport of energy outward within stars in turn 
causes their compositional and structural profiles to evolve. This occurs on 
a typically slow time-scale on the order of millions of years or longer (e.g. 



Kippenhahn fc Weigert 



1990 



Maeder 



20091 1 . Both the rate and outcome of 



the complex interplay of physical processes occurring within stars to drive 
their evolution depend somewhat sensitively on their mass. As a result, stellar 
evolution varies considerably between low- and high-mass stars. Below, we 
focus our description of this evolution to stars with masses < 2.0 M Q . 

The MS crosses the CMD from bottom right to top left, as shown in 
Figure [L~2l During the MS phase, the radius (and hence surface temperature) 



lb en. 


1991; 


, *- — 
Tout et al.. 


1996) 



and initial composition also contribute (e.g. 
a result, the lowest mass stars in a cluster occupy the bottom right end of the 
MS, with stellar mass increasing upward and blue-ward. Most stars spend the 
majority of their lives on the MS slowly converting hydrogen into helium in 
their cores. After slowly increasing their radius and luminosity by a factor o 



abou t 2.5-3, stars reach the end of the MS phase of their evolution ( lEggleton 



20061 ). This can roughly be defined as the point at which the hydrogen fuel 
within a central core containing about 10% of the star's mass has run out. 
Found at the left-most tip of the MS in the CMD, called the main-sequence 
turn-off (MSTO), the most massive MS stars are in the process of depleting 
their hydrogen fuel in and near their centres. At this point, loosely referred 
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to as the terminal-age main-sequence (TAMS), stars leave the MS and start 
moving to the right across the CMD, becoming increasingly red. 

The stellar evolution time-scales shorten considerably at this point, and 
the helium core begins to contract while the envelope expands. This causes 
a drop in the surface temperature of the star, causing it to move horizontally 
acro ss the CMD, crossing what is called the Hertzsprung gap (Figure 11.21) 



(e.g. 



Popper 



19801 ). The reason for this rapid reddening can be understood as 
follows. Once the surface temperature drops well below abo ut 10 4 K, the pr i- 



2006|). 



mary mechanism of energy transport in the envelope changes (lEggleton 
This is because the drop in temperature allows free electrons to recombine with 
free protons to form hydrogen atoms, which can then act to absorb out-going 
radiation. This presents a challenge for the star since energy still needs to 
be transported outward from its nuclear-burning centre, yet radiation can no 
longer leak out freely. To com pensate, convec t ion ta kes over as the dominant 



form of energy transport (e.g. 



Bohm-Vitense 



19581 ). As the convective base 



of the envelope deepens, stars of a given luminosity and mass converge to an 
approximately unique radius. The change in stellar radius that occurs during 
this short-lived phase of evolution can range from a factor of less than two fo r 



1991)). 



very low-mass stars to a factor of > 100 for very massive stars (jlben 
This marks the beginning of the red giant branch (RGB) phase of evolution, 
and the base of its corresponding sequence in the CMD (Figure fL2i) . 

Due to the unique envelope structure created by the deepening of the 
convective envelope, stars of a given m ass but different luminosities must li e 



on a roughly vertical locus in the CMD (IHayashi. Nishida fc Sugimoto 



19621 ). 



This structure also causes the stellar lumi nosity to in crease steeply as a func 



tion of the mass of the helium core (e.g. 



Ib en 



19681 ). As hydrogen burning 
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progresses in a shell immediately outside the core, the helium that is produced 
rains down onto it, slowly increasing its mass. This in turn causes the star 
to brighten by up to several orders of magnitude and, in so doing, ascend the 
RGB in the CMD. 

Eventually, the mass and temperature of the helium core reach ~ 0.47 
M and ~ 10 8 K, respectively, although the prec ise values depen d on both 



20061 ). It is at 



the total stellar mass and chemical composition ( lEggleton . 
this point that helium ignites, producing mainly carbon at first, although 
later the production of oxygen takes over. The horizontal branch (HB) in 
CMDs (Figure 11.2ft corresponds to the core-helium burning phase. It is not 
clear why HBs show such a large spread in their colours since, in principle, 
both the luminosity and surface temperature should be roughly constant for 
core helium-burning stars. This mysterious and ubiquitous feature of CMDs 
marks a significant gap in our understanding of stellar evolution theory. We 
will return to the curious HB morphologies observed in the CMDs of old star 
clusters in Section 11.1.2.11 

Within about 10 8 years of the onset of core helium-burning, low-mass 
stars will typic ally possess a core composed of carbon and oxygen with a mass 



of ~ 0.55 M (llben 



19741 ). The core is surrounded by a much less massive 
shell composed primarily of helium and next to no hydrogen. This, in turn, 
is surrounded by a n envelope of unprocessed material extending out to 30-50 



R Q from the core (IMaeder 



20091 ) . Two burning shells now power the star as 
it continues its rise in luminosity above the RGB. The sequence in the CMD 
corresponding to this phase of evolution is called the asymptotic giant branch 
(AGB). As this phase progresses, the radius and luminosity of the star continue 
to grow and severe mass-loss typically occurs due to powerful stellar winds. 
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The life of a typical low-mass star ends when the last of its envelope is 
burnt, leaving its luminosity to plummet by several orders of magnitude and 
its surface temperature to increase dramatically. The product is an incredibly 
dense remnant co mposed of carbon and oxygen called a C/O white dwarf 



(WD) dlben 



19741 ). This the case fo r stars wi t h ma sses up to ~ 2 M Q when 



mass-loss on the AGB is factored in (jEggleton 



20061 ). For stars more massive 



than this, stellar death can be considerably more dramatic, ending in the 
ignition of carbon and a subsequent thermonuclear explosion. As previously 
explained, these stars are short-lived and do not occur commonly in old star 



clusters, which wi 
reader to 



De the focus of this thesis. Therefore, we refer the interested 



Clavtonl (119681 ) and 



Maederl (120091 ) for descriptions of the evolution 



of massive stars. 



1.1.2 Additional Features of Colour-Magnitude Diagrams 

There remain several features characteristic of CMDs that cannot be explained 
by standard single star evolution. In order to account for their existence, it is 
necessary to invoke the aid of other physical processes known to be operating in 
star clusters, such as binary star evolution and stellar dynamics. As described 
in the subsequent sections, these processes are thought to be responsible for 
producing various types of exotic stellar populations and multiple star systems. 
These curious objects typically appear as outliers in CMDs, and do not fall 
on any single star evolution tracks. We will discuss the various mysterious 
features of CMDs and types of stellar exotica known to populate star clusters, 
and review the currently favoured hypotheses for their origins. 
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1.1.2.1 Horizontal Branch Morphology 



The horizontal branches of star clusters in the Milky Way have been observed 
to display a range of morphologies. That is, most CMDs have a horizontal 
sequence that extends blue- ward from the top of the red giant branch, and this 
is referred to as the horizontal branch. However, the length and appearance of 
this sequence varies considerably from cluster-to-cluster and, in some cases, it 
makes a sharp vertical transition to dimmer luminosities at its blue-most tip. 
CMDs are shown in Figure 11.31 for four different clusters, each of which has 
a distinct HB morphology. The source of this variance is unknown, however 
our curren t bes t gues s is mass-loss due t o stellar winds, most likely on the 



RGB (e.g. 



Iben 



1974 : 



Potter et al. 



20101 ). The more mass- loss that occurs, 



the smaller the mass of the remaining envelope. The surface temperature 
increases with decreasing envelope mass since this leaves more of the hot core 
exposed. Consequently, the amount of mass-loss is thought to determine the 
colours of HB stars, with stars that experience the greatest loss in mass ending 
up the bluest. In order to reproduce the observed surface te mperatures of th e 



hottest HB stars, the envelope must lose a mass of 0.2-0.3 M flEggleton 



2006|). 



This seems reasonable and agrees rather well with empirical estimates (e.g. 



Judge fc Stencel 



19911 ) . It is not yet clear, however, why stars of comparable 
mass in the same cluster would experience such different degrees of mass- 
loss. It is also unknown how or why this mass-loss distribution should change 
from cluster-to- cluster, a l thoug h the rate of mass-loss is thought to depend on 



metallicity (e.g. 



Reimers 



19751 ) 



Previous studies have confirmed that the observed differences in the 



HBs of Milky Way clusters are related to metallicity ( jSandage fc Wallerstein 



12 



Ph.D. Thesis — N. Leigh 



McMaster - Physics & Astronomy 




Figure 1.3 Colour-magnitude diagram sections for the Milky Way globular 
clusters NGC 104 (top left inset), NGC 6205 (top right inset), NGC 1261 
(lower left inset) and NGC 6934 (lower right inset). Note the distinct HB 
morphology characteristic of each cluster. Blue stragglers are also clearly visi- 
ble just brighter and bluer than the MS TO. The data used to create this plot 
were taken from lSarajedini et al.l (120071 ). 



19601 ). This does not tell the whole story, however, since at least one additional 
parameter is required to explain the spread in their colours (see Figure [L~3]) . 
Many cluster properties have been suggested as possible Second and Third 
Parameters. These include age, the cluster luminosity and central density. 
Unfor t unately, no definitive candid ates have as of yet been identified (e.g. 



Rood 



1973 



Fusi Pecci et al. 



19931 ). Notwithstanding, observations have re- 



vealed several peculiar trends for thi s curio us stellar population in individual 



clusters. For instance, 



Saviane et al. 



( 119981 ) presented evidence that blue HB 
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stars could be more ce ntrally concentrated than red HB stars in the cluster 



NGC 1851. Conversely, 



Cohen et al. 



( 119971 ) showed that blue HB stars could 



be centrally depleted relative to other stellar types in the cluster NGC 6205. 
To date, no clear evidence has been found linking the spatial distributions of 
HB stars to any global cluster properties. 

The bluest HB stars are often called extreme horizontal branch (EHB) 
stars. It has been suggested that RGB stars can shed their entire envelopes and 



still ignite helium as EHB stars (jd'Cruz et al 



19961 ) . Moreover, it is thought 



that EHB star s could skip the AGB phase and evolve into WDs directly (e.g. 



Maeder 



20091 ). EHB stars also go by the name of sub-dwarf B (sdB) stars. 



Both are core-helium burning stars with very small outer envel opes. Interest 



ingly , EHB stars in the field tend to have binary companions (IMaxted et al. 



200ll ). This could be interpreted as evidence that the presence of a binary 



companion is responsible for causing the dramatic loss of envelope mass. On 



the other hand, most EHB stars in c 



nary companion (e.g. 



Moni Bidin et al. 



usters 


rave 


)een i 


• 


2006 


• 


2009. 


2011) 



}een shown to lack a bi- 



20H|). Consequently, the 



mechanism responsible for (at least) their extreme mass-loss is still less clear. 
1.1.2.2 Blue Stragglers 



First discovered by ISandagd (119531 ). blue straggler stars occupy the region of 
the CMD that is just brighter and bluer than the MSTO (see Figure [L3]) . That 
is, they appear as an extension of the main-sequence. And yet, if all of the 
stars in a cluster were born at more or less the same time, then standard single 
star evolution predicts that this region of the CMD should be bare. Normal 
MS stars that are of sufficiently high mass to occupy this region should have 
long ago evolved onto the RGB and beyond. The presence of BSs in star 
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clusters therefore defies the predictions of single star evolution theory, and 
signifies another one of its mysteries. 

This puzzle can perhaps be solved if blue stragglers are produced via 
the addition of fr esh hydrogen to the cores of normal main-sequence stars (e.g. 



Sills et al. 



200ll ). This is the currently favoured origin for BSs. It can oc- 
cur via multiple channels, most of which involve the mergers of low-mass MS 
stars. Stars in binaries can be driven to merge if enough orbital angular mo- 
mentum i s lost. This can be mediated by dynami cal interactions with other 



star s (e.g. 



Leonard!. 



1989 



Ivanova fc Taam 



Chen fc Hani. 



2007 



Leonard &: Linnell 



20031 ). tidal dissipation (e.g. 



1992), magnetized stellar winds 



Cleary fc Monaghanl. 



2008 b]) or even an outer triple companion (e.g. 



19901: 



Perets Sz Fabrycky 



Fabrycky fc Tremaine 



20091 ). Alternatively, stars can collide directly, al- 



though this is a lso thought to usually be mediated by multiple star systems 



( Leonard 



19891 ). Stars in close binaries can transfer mass if their orbital sep- 
arations become sufficiently small for one of the components to over-fill its 
Roche lobe. This can also deliver fresh hydrogen to normal MS stars and in 
so doing produce BSs. 

Whatever the dominant BS formation mechanism(s) operating in dense 
star clusters, dynamical interactions are sufficiently common that they should 
play at least some role. For example, even if blue stragglers are formed as 
a result of binary evolution processes such as mass-transfer, the progenitor 
binaries themselves are expected to have expe rienced at least one dynamical 



inter a ction over the course of their l ifetim e (e.g. 



1989 



Davies. Piotto fc De Angeli . 



Hut fc Bahcalll . 



1983 



Leonard 



20041 ) . It follows that the study of BSs in 



star clusters offers an indirect means of probing the interplay between stellar 
evolution and stellar dynamics. 
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In order to explain recent observations of BS populations in star clus- 
ters, it has been suggested that several BS formation mechanisms operate si- 



tribu t ion for blue s t raggle rs (e.g. 



Ferraro et al.. 


1997 


1999. 


2004 


Lanzoni et al. 



2007 



Geller et al. 



20081 ). That is, the number of BSs is the highest in the 
central cluster regions, then falls off as the distance from the cluster centre 
increases but eventually rises again in the cluster outskirts. This has primar- 
ily been observed in globular clusters (GCs), which are particularly massive, 
dense and old star clusters usually found in the outer reaches of our Galaxy. 
One theory proposed to explain this result is that the blue stragglers found 

rereas those in the cluster 



2004T ). This 



in the cluster core were formed from collisions, w; 
outskirts were formed from binary mass-transfer (IFerraro et al. . 
hypothesis is motivated by the fact that the time-scale for collisions to occur 
between stars is very short in the core, but drops off considerably in the cluster 



19891 ). Alterna- 



outskirts where the stellar densities are much lower ( iLeonard 
tively, it has been suggested that the blue stragglers in the cluster outskirts 
could have been formed in the core but were lat er kicked out as a result of dy- 
namical encounters involving binary stars (e.g. 



Mapelli et al. 



20061). 



Sigurdsson fc Phinne 



1993 



A spectroscopic survey of the cluster NGC 188 performed by 



Mathieu &: Geller 



(200 9J) revealed t h at mo st, if not all, of its BSs have binary companions. 



Mathieu fc Geller! (J2009J) proposed that their results are consistent with the 



general picture that both mass-transfer and mergers are simult aneously pro- 



ducing BSs in NGC 188, which is also consistent with the results of 



Chen fc Han 



(l2008al ). The latter authors fit theoretical stellar evolution tracks to the ob- 



served colour-magnitude diagram for NGC 188. Based on their results, they 
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argue that most BSs in this cluster are too massive to have formed from mass- 
transfer alone. If true, this suggests that many BSs in NGC 188 must be 
the products of stellar me rgers. In support of this. 



argued that the results of 



Perets fc Fabrvckvl (l2009h 



Mathieu &: Gellerl (120091 ) can be explained if triple 



stars play an important ro le in BS formation by acting as catalysts for mergers 



(IPerets fc Fabrycky 



20091 ). Triple stars can evolve internally by transferring 



angular momentum between their inner and outer orbits. This can ca use the 



eccentricity of the inner binary to increase dramatically (IKozai . 



19621 ). Tidal 



friction can then act to reduce the orbital separation by removing orbital en 



ergy from the inner binary at each periastron passage ( iFabrycky fc Tremaine 



20071 ) (the term periastron refers to the point of closest approach for an ec- 
centric binary orbit). 

Blue stragglers have also been found in the field of our Galaxy. Obser- 
vations of these BSs sugges t that both mass - transf er and binary coalescence 



are occurring. For example. 



Di Stefano et al. 



(120101 ) recently reported the dis- 



covery of two BSs with white dwarf companions in the field of our Galaxy. The 
authors suggest that these BSs were formed from mass-transfer when the white 
dwarf progenitors evolve d to ascend the RG B and in so doing over-filled their 



Roche lobes. Moreover, iBrown et al.l (120 lOh discovered a ~ 9.1 M Q hyperve 



locity MS star in the field. The authors argue that it must be a blue straggler 
since its flight time from the Milky Way (MW) exceeds its MS lifetime. They 
conclude that this BS must be the product of a close binary that coalesced 
sometime after being ejected from the Galactic centre. This could explain 
its hypervelocity provided its ejection was caused by a dynamical interaction 
involving the central massive black hole. 

Several spectroscopic studies have also been performed to study indi- 
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vidual BSs in clusters. For example, at least two BSs in the globular clus- 
ter 47 Tuc have been foun d with masses exceeding tw i ce that of the main - 



sequence turn-off (MSTO) (IShara. Saffer fc Livio 



1997 



Knigge et al. 



20081 ). 



The large masses of these BSs suggest that they are the products of the merg- 
ers of two or more low-mass MS star s. Another similarly massive BS was 



reported by 



van den Berg et al. 



(I200ll ) in the old open cluster M67. This BS 



is though t to be a membe r of a triple star system that contains not one, but 



two BSs (ISandquist et al. 



20031 ). BSs have been identified in both binaries 
and t riples in several other cl usters as well, including the open cluster NGC 

201oL 



6819 (ITalamantes et al. 



1.1.2.3 Binary Stars 

Unresolved binary stars pollute the CMDs of star clusters, which become pep- 
pered with peculiar outliers that do not fall on any single star evolution tracks. 
This is because clusters are located at sufficiently great distances from us that 
the components of a given binary pair are indistinguishable and they appear 
as a single object. The combined light of the binary components blends to- 
gether, producing a combination of colour and brightness that ordinary single 
stars are never expected to display over the course of their evolution. This 
typic ally creates a second ary sequence in cluster CMDs located above the MS 



[e.g. 



Milone et al. 



20081 1 . The reason for this is that m ost cluster binaries 



are composed of normal MS stars (e.g. 



Geller et al. 



20081 ) . and the combined 



light of their components makes these binaries appear slightly brighter than 
ordinary single MS stars. 

There are many examples of objects that are almost certainly bonafide 
cluster members but appear in curious locations in the CMD. Examples in- 
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elude red stragglers (e.g. iKaluznyl . 120031 ). yellow str agglers fe.g.lLat 



laml 120051) 



2003|). Red 



and the sub-subgiant branch stars found in M67 (IMathieu et al. 
stragglers appear to the right (red-ward) of the RGB, yellow stragglers appear 
to the left (blue-ward) of the RGB and above the Hertzsprung gap, and sub- 
subgiant branch stars appear below the Hertzsprung gap. Most of these are 
thought to be unresolved binaries or triples and, in some cases, their curious 
CMD locations are speculated to be due to a recent episode of mass-transfer. 

1.2 Stellar Dynamics in Star Clusters 

Star clusters are in many ways analogous to stars themselves. They are both, 
in effect, self-gravitating spheres of interacting particles. In both cases, gravity 
plays a key role in driving the structural and compositional changes that char- 
acterize their evolution. Virial equilibrium is a central physical principle that 
deepens this analogy. This condition must be satisfied in order for any self- 
gravitating system to achieve a state of dynamic equilibrium, at which point 
the inward pull of gravity is balanced by the outward push of the pressure 
endowed to the system by the relative motions of its particles. The condition 
for virial equilibrium in a star cluster can be expressed in the following form: 



2T + W = 0, 



(i.i; 



where T is the total kinetic energy of the system and W is its total gravitational 
energy. Although additional physics must be factored in, such as the energy 
generated fro m nuclear-burning, a s imilar criterion for virial equilibrium exists 



for stars (e.g. 



Chandrasekhar 



19391 1. It follows that star clusters are virialized 
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systems composed of objects which are themselves virialized. 

Both short- and long-range gravitational interactions between stars act 
to evolve the structures and compositions of star clusters on time-scales ranging 
from millions to billions of years. In the following sections, we describe the 
primary dynamical processes responsible for this evolution. 



1.2.1 Two-Body Relaxation 



Although many of the details are still unclear, star clusters are thought to be 
born from massive gas clouds that fragmen t on several length scales t o form 



a large number of stars (e.g. 



Lada fc Lada 



1995 



McKee fc Ostriker 



2003). 



It ends up that the presently observed structures of star clusters are largely 
independent of the details of this collapse. This is because gravity will quickly 
act to mix the stars into a gra vitationally-bound sp herical distribution for al- 



most any initial configuration ( IHeggie fc Hut 



20031). Once settled, most of the 



stars in a cluster will t r avers e stable orbits, often with rosette-shapes, through- 



out it (IHeggie fc Hut 



20031 ) . Their paths will typically remain more or less 
undisturbed during a single orbit. Over longer time-scales, however, the mo- 
tions of the stars are affected by the cumulative effects of distant gravitational 
encounters as well as the odd close encounter. 

Like stars, star clusters evolve via the slow diffusion of heat from their 
centres to their outer edges. Gravitational encounters between pairs of stars 
act as the mechanism for heat transport in much the same way collisions be- 
tween pairs of molecules govern the flow of heat throughout an ideal gas. In 
the case of star clusters, however, it is the cumulative effect of many weak, and 



therefore distant, encounters that dominates, as opposed to the odd strong or 
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close encounter (IHeggie Sz Hutl . 120031 ). This process is called two-body relax- 
ation. The time-scale for it to occur, called the relaxation time, is considerably 
longer than the time required for a given star to cross the length of the cluster, 
called the crossing time. A useful measure of the relaxation time for the entire 
cluster comes from the half-mass relaxation time t r h, which is calculated from 
average quantities inside the half-mass radius r^ (defined as the distance from 
the cluster centre containing half its total mass). This is given by: 



trh 



o.issiv 1 / 2 ^ 72 

(Gm) l / 2 lnk ' 



1.2) 



where m and N are the average stellar mass and the total number of stars 
inside r^, respectively. The Coulomb logarithm InA is the factor by which 
small-angle encounters (i.e. encounters for which the angle between the initial 
and final velocity vectors of the deflected star is small) are more effective 
than large-angle encounters in a star cluster with given density and velocity 
distributions. 

For a typical Milky Way GC, t r h is on the order of a billion years. 
GCs are among the ol dest objects in the Univ erse, with ages usually ranging 



from ~ 9-12 Gyrs (e.g. 



De Angeli et al. 



20051 ). Consequently, we expect that 



most will have had sufficient time for two-body relaxation to have played a 
significant role in shaping their present-day features. Indeed, two-body relax- 
ation has been shown to dominate cluster evolution for a s ignific ant fraction 



of the lives of old MW GCs (e.g. 



Gieles. Heggie fc Zhao 



201lh . Other ef- 



fects also play their part. For example, mass- loss due to stellar evolution has 
also been shown to affect the dynamical evolution of star clusters, although 
its primary role is played during their early evolutionary phases when mas- 
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sive stars with power f ul ste l 



Chernoff Sz Weinberg 



1990 



ar winds are still present (e.g. lApplegatd . Il986 



Fukushige fc Heggie 



19951 ). The gas expelled by 



these winds can be significant. It escapes from the cluster, which expands in 
response to the loss in mass. This serves to delay the evolutionary progression 
induced by two-body relaxation. 

Equation 11.21 provides only a rough guide since the cumulative effects 
of two-body encounters depend on the stellar mass. This brings us to an- 



other centra 
(e.g. 



Henon 



concept of stellar dynamics : equipartition of kinetic energies 



1969 



Giersz fc Heggie 



19961 ). This is the tendency for all stars 



in a cluster to end up with comparable kinetic energies. This results from the 
fact that, during an individual gravitational encounter, the more massive star 
will impart a net positive acceleration to the less massive star, increasing its 
speed. This comes at the expense of the kinetic energy of the more massive 
star, which receives a corresponding net deceleration and slows down. In sys- 
tems where stars of widely differing mass occur, which is the case for real star 
clusters, the cumulative effects of these gravitational encounters cause stars 
with masses greater than the average stellar mass to lose kinetic energy. As 
a result, massive stars tend to slow down and sink to lower orbits within the 
cluster potential. At the same time, stars with masses less than the average 
stellar mass tend to gain kinetic energy and speed up, causing them to rise 
within the potential well of the cluster. In conjunction with close encounters 
occurring primarily within the dense cluster core, this process contributes to a 
steady stream of stars escaping from the system. The important point is that 
the stars that comprise clusters are all born with comparable velocities, but 
over time they tend toward a state in which they all have comparable kinetic 
energies. 
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The tendency for the most massive stars in a cluster to accumulate in 
the central regions and low-mass stars to be dispersed to wider orbits is called 
mass segregation. The time-scale for it to occur is typically short compared 
to the half-mass relaxation time. This can be understood as follows. The 
time-scale on which two-body encounters operate on an avera ge star o 
m has been shown to be well-approximated by Equation 11.21 (ISpitzer . 
In general, however, the ti me-scale on which this process o perates on a star of 



mass 



1983). 



mass M » m is given by ( iVishniac 



1978 



Spitzer 



1983) 



m 
M 



(1.3) 



This implies that the most massive stars in a cluster will segregate into the core 
on a time-scale much shorter than t^. It follows that most GCs will become 
fully mass segregated very early on in their 



2008 



McMillan. Vesperini fc Portegies Zwart 



etim es (e.g. 



Gaburov fc Gieles 



20071 ). This means that the spa- 



tial distributions of their stars become stratified according to their mass. In 
other words, the probability of finding a star in the cluster outskirts is inversely 
proportional to its mass. It is important to note, however, that we still see a 
range of masses at all cluster radii. 

Two-body relaxation is one of the primary mechanisms responsible for 
the transfer of heat outward within a cluster. This causes the inner cluster 
regions to contract, and the outer regions to expand. To see why this is the 
case, consider what occurs if energy is transferred from the inner part of an 
isolated system of gravitationally- interacting particles to the outer part. The 
inner part should cool. However, gravity then causes these stars to drop to 
lower orbits within the cluster potential and, in so doing, speed up. The net 
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effect of this is that the average velocity of stars in the inner part increases, 
causing its temperature to increase. The energy transferred to the outer part 
causes it to heat up as well, but it expands in response to the addition of heat. 
In general, the increase in the temperature of the outer part is larger when 
it is smaller in size. This creates a temperature gradient between the inner 
and outer parts. As a result, more heat flows outward and the temperature 
gradient is enhanced. This leads to a runaway effect that has been dubbed the 
gravothermal catastrophe. For typical MW GCs, a phenomenon known a s core 



collapse occurs within about 10-20 half-mass r elaxation times (IHenon 



1969 



Spitzer fc Shulll . 



1975 



Hut fc Djorgovski 



19921 ) . This marks an enhancement 



in the central density by several orders of magnitude. 

Theoretical models suggest that the time-scale for core collapse to occur 
is often longer than the age of the Universe. Therefore, most MW GCs should 
currently be in a phase of core contraction. This evolutionary phase will only 
come to an end once the central density become s sufficient l y high for encoun 



ters involving binaries to halt the pro cess (e.g. 



Hut 



1983 



Goodman fc Hut 



1993 



Fregeau. Ivanova fc Rasiol . 



20091 ). The net effect of these interactions is 



for single stars to steal energy from the orbits of the binaries. This imparts 
additional kinetic energy to the single stars at the expense of reducing the 
orbital separations of the binaries. In turn, this provides a heat source for the 
cluster, and ultimately halts the collapse of the core. 



1.2.2 Small-N Dynamics 

Our discussion of cluster evolution has naturally brought us to the issue of 
small-N gravitational dynamics. These are short-range interactions involving 
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only a few stars for which the distance of closest approach can be compara- 
ble to the stellar radii. During single-binary and binary-binary encounters, 
resonant interactions often occur in which the stars remain bound for many 
crossing times. Two objects (where an object can refer to a single, binary 
or even triple star) approach one another and a series of close, or strong, 
gravitational interactions ensue. A number of outcomes are possible. These 
include exchanges between the components of binaries or even their complete 
dissociation, the formation of triples, as well as stellar collisions and mergers. 

Numerous scattering experiments have been performed to explore the 
outcomes of binary-binary and, in par t icular , single-binary encou nters (e.g. 



McMillan 



1986 



Sigurdsson fc Phinney 



1993 



Fregeau et al. 



2004J ). Most of 



the earliest of these studies were performed in the point-particle limit. Conse- 
quently, they ignored the the stars' finite sizes, and so neglected the importance 
of taking into account t he dissipative effects of t idal interaction s and direct con- 



tact between stars (e.g. 



Hut fe Verbunt 



remed ied by, for instance, 



1983 



McMillan et al. 



Mikkola 



( 19871 ) and 



\ 983). This was later 



Cleary fc Monaghan 



(Il990l ). Encounters involving four or more stars require longer integration 



times to run the simulations to completion, and involve a large number of free 
parameters. As a result, few studies have been conducted to explore the out- 
comes of binary-binary encounters or interactions involving triple systems. To 
date, none of these have considered the finite sizes of the stars in a completely 
realistic way. 

It has been known for some time that encounters betwee n stars, and 



even direct collisions, can occur frequently i n dense stellar systems (jHills fc Day 



1976 



Hut fc Verbunt 



1983 



Leonard 



19891 . e.g.). In the cores of globular clus- 



ters (GCs), the time between collisions involving two single stars can be much 
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shorter than the cluster lifetime ( [Leonard 



19891 ). The time between encoun- 



ters involving binary stars can be considerably shorter still given their much 
larger cross-sections for collision. In globular and, especially, open clusters 
with high binary fractions, mergers are tho ught to occur frequently d uring res- 
onant interactions involving binaries (e.g. iLeonard fc Linnelll . Il992l ). What's 
more, collision products have a significant probability of undergoing a sec- 
ond or even third collision collision during a given single-binary or binary- 
binary interaction. This is because the impact causes the collision product to 
expand, increasing t he cross-section for a subsequent collision to occur (e.g. 



Fregeau et al. 



20041). 



Small-N dynamical interactions play a number of important roles in star 
cluster evolution. By acting as an important heat source for clusters, binaries 
become modified by dynamical interactions. Although the details are not yet 
clear, this can change the distribution of orbital parameters of binary pop ula- 



tions in clusters (e.g. 



Hut fc Bahcall 



1983 



Sigurdsson &: Phinney 



1993|). In 



turn, this could have important implications for binary star evolution by, for 
example, stimulating mass-transfer. Observational evidence has been found in 
support of this. Previous studies have reported evidence that the sizes of some 
binary populations thought to be undergoing m a ss-tra nsfer are correlated with 



the rate of stellar collisions (e.g. 



Poolev &: Hut 



2006|). 



1.3 Where To Go From Here 

We have a good working knowledge of both stellar evolution and stellar dy- 
namics in star clusters. But our understanding is far from complete. This 
brings us to the question: Where do we go from here? The interaction be- 
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tween stellar evolution and stellar dynamics in star clusters remains a largely 
unexplored frontier. What observational effects might we expect to find? How 
will these affect star cluster evolution, and what implications could they have 
for the history of our Galaxy? The search for answers to these questions will 
be the primary focus of this thesis. 

The techniques presented in this thesis are motivated by our current 
understanding of both stellar evolution and stellar dynamics, in addition to 
evidence that suggests that their interaction could account for the origins of 
several mysterious stellar populations. Gravitational dynamics is the common 
theme unifying all of these methods. They will be applied to observations 
of dense star clusters in order to study, among other things, stellar mergers. 
An emphasis will be placed on the development of statistical and theoretical 
tools that can be applied to a number of other astrophysical sub-disciplines 
with analogous dynamical processes and population statistics. These tools 
will be used in this thesis to further our understanding of the various channels 
via which stellar evolution and stellar dynamics interact in star clusters, as 
well as to quantify the implications of these effects for observations of stellar 
populations. 

In Chapter [2} we introduce a new adaptation of the classical mean free 
path approximation. With it, we compare the rates of gravitational encoun- 
ters occurring between single, binary and triple stars using observations of 
real star clusters. This has allowed us to outline a systematic methodology 
that can be used to constrain the dynamical origins of observed multiple star 
systems containing merger products. In Chapter [3J we introduce a statistical 
technique that can be used to compare the relative sizes of the different stellar 
populations in a large sample of star clusters spanning a diverse range of prop- 



27 



Ph.D. Thesis — N. Leigh 



McMaster - Physics & Astronomy 



erties. In order to study the effects had by the cluster dynamics on each of 
the different stellar populations, we apply our method to a large sample of 56 
Milky Way GCs. In Chapter HI we refine this technique and apply it to a new 
sample of 35 GCs compiled using much more sophisticated observations taken 
from the ACS Survey for Globular Clusters. Using this new data, we present a 
homogeneous catalogue for the different stellar populations, along with a sim- 
ple prescription to select stars out of the CMD belonging to each population. 
In Chapter [5l we present an analytic model for blue straggler formation in 
globular clusters. We compare its predictions to observed BS numbers taken 
from our stellar population catalogue using a new statistical technique. With 
this method, we constrain the dominant blue straggler formation mechanism 
operating in GCs. In Chapter [61 we present a new method of quantifying 
cluster-to-cluster differences in the stellar mass functions of a large sample of 
clusters. We apply our technique to the ACS data, and use it to constrain 
the degree of universality of the initial mass function for our sample. Given 
the very old ages of Milky Way globular clusters, this has important implica- 
tions for our understanding of star formation in the early Universe. Finally, in 
Chapter El we summarize our results and discuss directions for future research. 
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Chapter 



2 



An Analytic Technique for 
Constraining the Dynamical 
Origins of Multiple Star Systems 
Containing Merger Products 



2.1 Introduction 



It has been known for some time that encounters, and even direct col 



can occur frequently between stars in d ense stellar systems ( iHills fc Day 



lsions. 



1976 



Hut fc Verbunt 



1983 



Leonard 



1989 



e.g.). In the cores of globular clusters 



(GCs), the time between collisio ns involving tw o single stars can be much 



19891 ). The time between encoun- 



shorter than the cluster lifetime ( [Leonard! 
ters involving binary stars can be considerably shorter still given their much 
larger cross sections for collision. In globular and, especially, open clusters 
(OCs) with high binary fractions, mergers are thought to occur frequently 
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during resonant interactions involving binaries (e.g. iLeonard fc Linnelll . Il992l ). 
What's more, collision products have a significant probability of undergoing 
more than one collision during a given single-binary or binary-binary interac- 
tion since the initial impact is expected to result in shock heating followed by 
adia batic expansion, incre asing the cross section for a second collision to occur 



(e.g. 



Fregeau et al. 



2004|). 



Several types of stars whose origins remain a mystery are speculated 
to be the products of stellar mergers. Blue stragglers (BSs) in particular 
are thought to be produced via the addition of fresh hydrogen to the cores 
of low-mass main-sequence (MS) stars. Recent evidence has shown that, 
whatever the dominant BS formation mechanism(s) operating in both glob- 
ular and open cluster s , it is likely to in some way depend on binary stars 



(IKnigge. Leigh fc Sills 



2009 



Mathieu fc Oellerl . 



20091 ). The currently favored 



mechanis ms include colli sions during single-binary and binary-binary encoun- 



ters (e.g. 



Leonard, 



system (e.g. 



1989)), mass t ransfer between the components of a binary 



Chen &: Han 



2008a 



bj) and the coalescence of two stars in a close 



binary due to perturbations fro m an orbiting triple companion (e.g. 



2006 



Perets fc Fabrycky 



20091 ) 



Eggleton 



A handful of spectroscopic studies have revealed that in some GCs 
there exist BSs with masses exceeding twice that of the MS turn-off (e.g. 



Shara. Saffer fc Livio 



1997 



Knigge et al. 



20081 ). Such massive BSs must have 



been formed from the mergers of two or more low-mass MS stars since they 
are too massive to have been formed from mass transfer. In a few cases, this 



can also be argued 



or entire BS populations using photometry. For instance, 



Chen fc Hani (j2008bl ) performed detailed binary evolution calculations to study 



dynamical stability during mass transfer from an evolving giant star onto a MS 
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companion. Based on their results, it can arguably be inferred that most BSs 
in NGC 188 are sufficiently bright that they probably could not have formed 
from mass transfer alone. If true, this suggests that most of these BSs must 
be the products of stellar mergers. Regardless of the dominant BS formation 
mechanism(s) operating in dense star clusters, dynamical interactions should 
play at least some role. For example, even if blue stragglers are formed as a 
result of binary evolution processes such as mass transfer, the progenitor bina- 
ries themselves should have been affected by at least one dynamical interaction 
over the course of their lifetime. 

Numerous scattering experiments have been performed to explore the 
outcomes of binary-binary and, in par t icular , single-binary encou nters (e.g. 



McMillan 



1986 



Sigurdsson fc Phinney 



1993 



Fregeau et al. 



20041 ) . Most of 



the earliest of these studies were performed in the point-particle limit, ignoring 
altogether the often non-negligible i mplications of the stars' finite sizes (e.g. 



Hut fc Bahcall 



1983 



Mikkola 



19831 ) . Later, more realistic simulations clearly 



demonstrated the importance of taking into account the d issipative effects 



of tidal interactions and direct contact between stars (e.g. 



McMillan et al 



1987 



Cleary fc Monaghan 



19901 ). As a result of the increased number of free 
parameters for the encounters and the longer integration times required to run 
the simulations to completion, few studies have been conducted to explore the 
outcomes of binary-binary encounters or interactions involving triple systems. 

In this chapter, we introduce an analytic technique to constrain the 
most probable dynamical origin of an observed binary or triple system con- 
taining one or more merger products. Provided the observed system is found 
within a moderately dense cluster environment with binary and/or triple frac- 
tions of at least a few percent, the probability is often high that it formed from 
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a merger during an encounter involving one or more binary or triple stars. In 
Section 12. 2[ we present an equation for energy conservation during individual 
stellar encounters and outline the process for applying our technique. Specif- 
ically, we present a step-by-step methodology to evaluate whether or not an 
assumed dynamical history could have realistically produced an observed sys- 
tem and describe how to determine the most probable dynamical formation 
scenario. In Section |273| we apply our technique to a few observed binary and 
triple systems thought to contain merger products, in particular a triple system 
that is thought to contain two BSs and the peculiar period-eccentricity distri- 
bution of the BS binary population in NGC 188. We discuss the implications 
of our results in Section 12.41 



2.2 Method 



In this section, we present a general prescription for conservation of energy 
during stellar encounters. We will limit the discussion to typical interactions 
thought to occur in globular and old open clusters, although our technique 
can be generalized to any choice of parameter space. The types of encoun- 
ters of interest in this chapter will predominantly involve low-mass MS stars 
with relative velocities at infinity ranging fr om < 1 km/s to ~ 10 km/s (e.g. 



Leonard 



1989 



Sigurdsson fc Phinney 



19931 ) . Our technique describes how to 



isolate the most probable dynamical formation history for an observed binary 
or triple containing one or more merger products by providing an estimate 
for the time required for a given interaction to occur in a realistic cluster 
environment. 

We begin by assuming that an observed system was formed directly 
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from a dynamical interaction (or sequence of interactions). In this case, the 
observed parameters of the system provide the final distribution of energies for 
the system resulting from the interaction(s). After choosing an appropriate 
dynamical scenario (i.e. whether the objects involved in the interaction(s) are 
single, binary or triple stars), we can work backwards using energy conservation 
to constrain the initial energies going into the encounter. This provides an 
estimate for the initial orbital energies and therefore semi-major axes of any 
binaries or triples going into the interaction. This in turn gives the cross section 
for collision and hence the time required for the hypothesized interaction(s) to 
occur. 

Since the formation event must have happened in the last tbs years, 
where tbs is the lifetime of the merger product, a formation scenario is likely 
only if the derived encounter time-scale is shorter than the lifetime of the 
merger product (s). Conversely, if the derived encounter time-scale is longer 
than the lifetime of the merger product (s), then that dynamical formation 
scenario is unlikely to have occurred in the last tbs years. In general, the 
shorter the derived encounter time-scale, the more likely it is that one or 
more such encounters actually took place within the lifetime of the merger 
product (s). Finally, if the derived encounter time-scale is longer than tbs 
for every possible dynamical formation scenario, then a dynamical origin is 
altogether unlikely for an observed multiple star system containing one or 
more BSs. Either that, or the encounter time-scales must have been shorter 
in the recent past (or, equivalently, the central cluster density must have been 
higher) . 
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2.2.1 Conservation of Energy 

Consider an encounter in which at least one of the two bodies involved is a 
binary or triple star. Though a complex exchange of energies occurs, energy 
must ultimately be conserved in any dynamical interaction. The total en- 
ergy that goes into the encounter must therefore be equal to the total energy 
contained in the remaining configuration: 

Ni Ni 

y^fii(i,,a>i) + Wj(mi, X h Z h Tj) 

i i 
Ni S z 

Ui{nii, X u Z h Ti ) + T A M ch Mj)+ 

* 3 

Mi N f 

^2 e k(Vk, M k , a k ) = ^ (2.1) 

k ii 

Nf Nf 

^ ^ W^iiijTl'ii) Xa, Znj Tjj) 4~ ^ ^ UaiTflu, Xa, Za, Tjj) 

ii ii 

Sf M; 

+ Z~2 T iA M 3i) + Z~L e kk(^kk, M kk, a k k) + A, 

jj kk 

where Ni, Si and Mj are the total number of stars, objects (single, binary, 
triple or even quadruple stars) and orbits, respectively, that went into the 
encounter. Similarly, Nf, Sf and Mf are the total number of stars, objects 
and orbits remaining after the encounter. 

We let Qi = Qi(Ii,Ui) represent the bulk rotational kinetic energy in 
star i, which is a function of the star's moment of inertia Ii and angular rotation 
rate ojc 

Wi) = -hul (2.2) 
The moment of inertia is in turn a function of the density profile within a star, 
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which changes along with its internal struc ture and composition as the star 



evolves. The moment of inertia is given by (IClaret Sz Gimenez , 



1989|): 



8vr 



R 



3 

/3 2 mR 



o 

2_ D 2 



(2.3) 



where (3 is the radius of gyration. For example, a typical 1 M star in an 
old open cluster having an age of ~ 6 Gyrs has (3 = 0.241. For comparison, 



a 1 M Q st ar in a typical Mi 



(3 = 0.357 (jClaret fc Gimenez 



ky Wa y GC having an age of ~ 11 Gyrs has 



19891 ) . A large spread of rotation speeds have 



been observed for MS stars in o pen and globular c 



ranging from ~ 1 - 20 km s 1 (IMathieu fc Geller 



usters, with measured values 



20091 ). A 1 M star having 
a radius of 1 R Q and a rotation speed of 2 km s _1 with (3 = 0.241 has Q ~ 10 36 
Joules. 

We let Wi = Wi(rrii, Xi, Zi,Ti) represent the gravitational binding en- 
ergy of star i, where m; is the star's mass, Xi is its initial hydrogen mass 
fraction, is its initial metallicity and r, is its age. For a spherical mass with 
a density distribution p«(r), the gravitational binding energy is given by: 



Wi(mi,Xi,Zi,Ti) 



16tt 2 G 



Pi(r)r 4 dr 



(2.4) 



-Sim^X^Z^r^G-^, 



where the parameter 5 is chosen to reflect the structure of the star and is 
therefore a function of its mass, age and chemical composition. For instance, a 
typical 1 M star with (X, Z) = (0.70, 0.02) in an old open cluster with an age 
of ~ 6 Gyrs has 5 = 1.892. For comparison, an older but otherwise identical 
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star in a typical Milky W ay GC with an age of ~ 11 Gyrs has 5 = 6.337 



( jClaret Sz Gimenez 



1989|) ■ Roughly regardless of age, this gives \W\ ~ 10 42 
Joules for a 1 M Q star with a radius of 1 R Q . 

We let Ui = Ui(rrii,Xi, Zi,Ti) represent the total internal energy con- 
tained in star i (i.e. the star's thermal energy arising from the random motions 
of its particles). By solving the equations of stellar structure, the total inter- 
nal energy of a purely isolated single star is uniquely determined by its mass, 
initial composition and age. Stars are made up of a more or less virialized 
fluid so that, ignoring magnetic fields, the gravitational binding energy of 
a star in hydrostatic equilibrium is about twice its internal thermal energy 



(I Chandr asekhar 



19391 ). Using this version of the virial theorem, a 1 M star 
having a radius of 1 R has U ~ 5 x 10 41 Joules. 

The total translational kinetic energy of object j (single, binary, triple, 
etc. star) is represented by Ty. 



T j (M ch M j ) = -M j vl 



(2.5) 



where Mj is the total mass of the object, and Vj is its bulk translational speed. 
The translational velocities of stars in clusters for which energy equiparti- 
tion has been achieved as a result of two-body relaxation obey a Maxwell- 
Boltzmann distribution, with th e heav iest stars typically having the lowest 



velocities and vice versa (ISpitzer 



19871 ). According to the virial theorem, the 



root-mean-square velocity v rms of the distribution depends on the total mass 
of the cluster M c i. It follows that the velocities of stars in a fully relaxed 
cluster are approximately determined by their mass and the total mass of the 
cluster. Assuming that the typical velocity of a < m > star is roughly equal to 
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Vrms, where < m > is the average stellar mass in the cluster, energy equiparti- 
tion can be invoked in some clusters to approximate the translational kinetic 
energy, and hence velocity, of a star or binary of mass Mf 



We note that Equation 12.61 can only be applied in clusters for which the half- 
mass relaxation time is much shorter than the cluster lifetime. At the same 
time, the tidal truncation of the velocity distribution must not be significant. 
We will return to this in Section 12.2.21 

For a 1 M Q star with a translational speed of 1 km s _1 (typical of stars 
in old open clusters), we find T ~ 10 36 Joules. For comparison, the same star 
traveling at a speed of 10 km s -1 (typical of stars in globular clusters) has 
T ~ 10 38 Joules. We can put this into perspective by equating Equation 12.51 
with Equation 12 A\ which shows that a direct collision between two 1 M stars 
would require an impact velocity of ~ 1000 km s _1 in order to completely 
unbind the merger remnant. 

The total orbital energy of orbit k is denoted e k , and is given by: 



where fi k = m^/ (mi+m 2 ) is the reduced mass of the orbit, M k = m!+m 2 is 
the total mass and a k is the orbital semi-major axis. A binary composed of two 
1 M stars with a period of 1000 days has |e| = 10 38 Joules. For comparison, 
an otherwise identical binary with a period of 1 day has |e| = 10 40 Joules. 
Since most stable triples are observed to have outer periods of ~ 1000 days 




(2.6) 



e k (fi k ,M k ,a k ) 



Gfi k M k 
2a k 



(2.7) 
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with close inner binaries ( jTokovininl . 119971 ; iPerets fc Fabryckyl . l2009l ). it follows 
that the total orbital energies of stable triples will be dominated by the orbital 
energies of their inner binaries. 

After the encounter occurs, the energies that collectively define the 
state of the newly formed system include the rotational kinetic (fljj(/jj,Wij)), 
thermal (Uu(mu, X iU Zu, To}) and gravitational binding (Wu{m iU Xu, Z iU m)) 
energies of the stars that are left-over, as well as the total translational kinetic 
(Tjj(rrijj)) and orbital (ekki^kk, Mkk, a>kk)) energies of any left-over stars, bina- 
ries or triples. The internal and gravitational binding energies of the left-over 
stars once again depend on their mass, composition and age. The mass, com- 
position and evolutionary status of a star will decide how it responds to tidal 
interactions (and therefore how much tidal energy is deposited) since these 
are the principal factors that determine its internal s tructure, in pa r ticula r 



1996|). 



whether or not its envelope is radiative or convective (IPodsiadlowski 
Finally, we let A represent the energy lost from the system due to radiation 
and mass loss. We do not expect this term to be significant for the majority 
of the encounters we will consider since the rate of mass-loss from low-mass 
MS stars is small and the time-scales under consideration are relatively short 
compared to the lifetimes of the stars. Moreover, the velocity dispersions char- 
acteristic of the clusters we will consider are sufficiently small that we expect 
most collisions to have a relatively low impact velocity. Though significant 
mass loss can occur for high impact collisions as often occur in the Galactic 
centre, mass loss is only < 5% for the l ow impact velocit y collisions expected 



200ll ). Therefore, we will 



to occur in open and globular clusters ( pills et al 
henceforth assume A ~ 0. 

Depending on the initial parameters of the encounter, one or more terms 



47 



Ph.D. Thesis — N. Leigh 



McMaster - Physics & Astronomy 



in Equation 12 . 1 1 can often be neglected. For example, given that the rotational 
energies of typical MS stars in both open and globular clusters are several or- 



ders of magnitude smaller than the other terms in Equation 12. 11 we can neglect 
these terms for the types of encounters of interest. Even a 1 M star rotating 
at a rate of 100 km s _1 has a rotational energy of only ~ 10 39 Joules. This 
is considerably higher than ev en the highest rotation r ates observed for both 



BSs and normal MS stars (e.g. 



Mathieu fc Geller 



20091 ). From this, we expect 



typical rotational energies to be significantly smaller than the orbital energies 
of even moderately hard binaries. Furthermore, Equat ion 12 . 5 1 can be combined 
with Equation 12.61 and equated to Equation 12. 71 in order to define the hard-soft 
boundary for a given cluster. If a binary is hard, its initial orbital energy will 



outw eigh its initia 



ters ( IHeggie 



translational kinetic energy during typical stellar encoun- 



19751) . We do no t expect soft binaries to survive for very long 



in dynamically-active clusters ( IHeggie 



19751 ) (i.e. most binaries are hard), so 



that the translational kinetic energies of the stars, binaries and triples typi- 
cally found in OCs can usually be neglected when applying Equation 12. II This 
suggests Tj = 0. However, the final translational kinetic energies of the stars, 
binaries and triples left-over should not be neglected for encounter outcomes 
in which one or more stars are ejected with very high velocities. This can 
leave the remaining stars in a much more tightly bound configuration since 
stars ejected with high velocities carry off significant amounts of energy. Fi- 
nally, provided there are no very high impact collisions and tides dissipate a 
negligible amount of energy from the system, we also expect that: 



N r 



N, 



(2. 
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We can simplify our energy conservation prescription considerably for 
the majority of encounters occurring in old OCs. Neglecting the rotational 
kinetic energies of the stars, Equation 12.11 becomes: 



Mi M f 

£k(Hk, M k , a k ) = tkki^kk, M kk , a kk ) + 

k kk 
33 



(2.9) 



where we have assumed Tj = and A m is the amount of energy deposited in 
any merger products formed during the encounter. In other words, the term 
A m provides the required correction to Equation 12.81 resulting from internal 
energy being deposited in the merger product (s) as a result of collisions and 
tides. If no mergers occur or all of the orbital energy of the merging binary is 
imparted to the other interacting stars, A m = 0. 

Both the total linear and angular momenta must also be conserved 
during any dynamical interaction. This provides two additional constraints 
that must also be satisfied, however the total linear and angular momenta 
are vector quantities that depend on the angle of approach as well as the 
relative orientations of the objects. As a result, there are more free parameters 
to fit when trying to constrain the initial parameters of an encounter using 
conservation of linear and angular momentum. Therefore, it is considerably 
more difficult to extract information pertaining to the initial orbital parameters 
of an encounter using conservation of momentum than it is with conservation 
of energy. 

We have assumed that exchange interactions do not occur in applying 
our technique. Clearly, this assumption could be invalid for some systems. 
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If this is the case, then our method is also invalid. This suggests that our 
technique is ideally suited to clusters for which the encounter time-scales are 
comparable to the lifetime of a typical merger product. This maximizes the 
probability that BSs do not experience any subsequent dynamical interactions 
after they are formed. Of course, most clusters of interest are unlikely to satisfy 
this criterion. We can assess the probability that an exchange encounter has 
occurred for a particular system by calculating the different encounter time- 
scales and comparing to the lifetime of the merger product. If the time-scales 
are sufficiently short, then the possibility of an exchange interaction having 
occurred after the system's formation must be properly addressed. We will 
discuss this further in Section 12.41 

2.2.2 Generalized Approach 

In this section, we outline a step-by-step methodology to constrain realistic 
dynamical formation scenarios that could have resulted in the production of 
an observed stellar system containing one or more merger products. These 
steps are: 

1. We must first find qualitative constraints for the system's dynamical his- 
tory and, in so doing, converge on the most probable formation scenario. 
The choice of formation history should be guided by the observed proper- 
ties of the binary or triple system containing the merger product (s). The 
following guidelines can be applied to find the most probable scenario. 

First, the analytic rates for single-single (1+1), single-binary (1+2), 
single-triple (1+3), binary-binary (2+2), binary-triple (2+3) and triple- 
triple (3+3) encounters can be compared to obtain a rough guide as to 
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which of these encounter types will dominate in a given cluster. The 
total rate of encounters o f a given type in a cluster core is well approxi- 



mated by ((Leonard, 



19891 1 (see Appendix |A] for a more generalized form 



for this equation): 



(2.10) 



where Nq is the number of single, binary or triple stars in the core and 
no is the mean stellar, binary or triple number density in the core. The 
grav itational 



6 of 



Leonard 



y-focu sed cross section for collision a g f is given by Equation 
(119891 ) . Gravitationally-focused cross sections for the var- 
ious encounter types are provided in Appendix |A] along with the values 
assumed for their pericenters. 

In general, the number of single, binary and triple stars are given by, 



respectively, (1 - f b - f t )N c , f b N c an d f t N r , w 



total number of objects in the core ([Leonard 



r ere N r = 2/37rn r^ is the 



19891 ). fb is the fraction 



of objects that are binaries and ft is the fraction of objects that are 
triples. Assuming for simplicity that v re i is roughly equal for all types 
of encounters, the rates for two types of encounters can be compared to 
find the parameter space for which one type of encounter will dominate 
over another. These relations can be plotted in the fb — ft plane in 
order to partition the parameter space for which each of the various 
encounter types will occur with the greatest frequency, as illustrated in 
Figure 12.11 Given a cluster's binary and triple fractions, this provides a 
simple means of finding the type of encounter that will occur with the 
greatest frequency. Our results are in rough agreement with those of 
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Sigurdsson fc Phinneyl ( 119931 ) who found that single-binary interactions 
dominate over binary-binary interactions in clusters having core binary 
fractions /& < 0.1, and may dominate for up to 0.25-0.5 in some cases. 



E 0.6 



0.4 




0.-1 0.6 
Binary Fraction 



Figure 2.1 Plot showing the parameter space in the fb — ft plane for which 
each of the various encounter types dominate. Boundaries between regions are 
indicated by solid lines, each segment of which is obtained by equating two 
particular encounter rates using Equation [2.101 and the relevant cross sections 
derived using Equation 6 of lLeonardl (119891 ). We assume a t = 5ab and a>b = 90R 
in obtaining the relations between encounter rates. This is a reasonable choice 
for the ratio between the average binary and triple geometric cross sections 
given that the ratio between the outer and inner orbital semi- major axes of 
triples must be relative ly large (by a factor of > 10) in order for them to be 
stable dMardlind . [200lh . 



Second, the masses of the components of an observed binary or triple 
system provide a lower limit for the number of stars that could have 
gone into its formation. The minimum number of stars that must have 



52 



Ph.D. Thesis — N. Leigh 



McMaster - Physics & Astronomy 



merged to form a given collision product, labelled iV m j n) j, is equal to 
the integer nearest to and larger than the quantity m rem /mTo, where 
m rem is the mass of the merger remnant and vtlto corresponds to the 
mass of the main-sequence turn-off (MSTO). This assumes that nixo has 
not changed significantly since the dynamical formation of the system, 
which should be valid provided the merger products are significantly 
more massive than the turn-off so that their lifetimes are relatively short. 
The number of stars that went into an encounter JVj must therefore 
satisfy: 



where Mq is the number of normal stars (i.e. not formed from mergers) 
and M denotes the number of merger products. 

Third, an estimate for the average stellar mass, and therefore the masses 
of typical stars expected to undergo encounters, should be guided by 
a realistic stellar mass function for the host cluster. First of all, ob- 
servations have shown that a significant depletion of low-mass stars oc- 
curs in dynamically evolved clusters (i.e. those for which t rh <C t age , 



where t r h is the half-mass relaxation time and t age is the cluster age) 
since they are preferentially ejected from the cluster during close en- 



counters (e.g. 


von Hippel & Saraiedini 




1998 


Bonatto. Bica Sz Santos. 


2005; 


De Marchi. Paresce & Portegies Zwart 




2010 


). This suggests that 



OCs and the least massive Milky Way GCs should have stellar mass 
functions that, at least in their central cluster regions, appear eroded at 
the low- mass end. From this, we expect < m >< vtlto- Finally, most 
of the stars involved in dynamical encounters should have masses close 



M 




(2.11) 
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to or even slightly greater than the average stellar mass < m >. This is 
because gravitational focusing is strongest for the most massive objects, 
contributing to a shorter encounter time-scale. 

Finally, the most likely formation scenario will, strictly speaking, min- 
imize the total or cumulative encounter time-scale. However, the total 
time-scale required for a dynamical formation scenario that involves more 
than one encounter will typically be dominated by the second encounter. 
This is because, after an initial encounter has occurred between any two 
suitable objects to form a new stellar configuration, the time-scale for 
a second encounter to occur is given by the time required for the prod- 
uct of the initial encounter to experience a subsequent encounter. This 
increases the encounter time-scale by a factor Nq. In other words, Equa- 
tion [2TT0] provides the time required for any two of the specified objects 
to experience an encounter (two binaries, a binary and a triple, etc.). 
It can be multiplied by iVo to obtain an estimate for the time required 
for a specific object to experience an encounter. For example, to find 
the time required for a specific binary to experience an encounter with 
another single, binary or triple star, we must multiply by a factor fbN c . 
It follows that the time-scale for multiple encounters to occur is consid- 
erably longer than any of the single encounter time-scales. Therefore, 
unless the number of either binary or triple stars is very low, scenarios 
involving the fewest number of encounters are generally preferred since 
this tends to minimize the total time required for the encounter(s) to 
occur in a realistic cluster environment. 

The times between the various types of encounters can be derived using 
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Equation 12.101 and have been provided in Appendix [A] These time-scales 
can be multiplied by (1 — /&— ft)N c , fbN c or f t N c to find the time required 
for a particular single, binary or triple star, respectively, to encounter 
another object. T his yields time-scales t hat ar e in rough agreement with 



Equation 8-125 of 



Binney fc Tremaind (119871 ). As an example, we can 



multiply Equation IA.8I by /&iV c to find the time for a particular binary 
to experience an encounter with another binary. This gives: 



r 2 6 +2 = 2.7xl0 1 V f T 1 ( 



10 3 pc 
no 



.5Mr,\ ( 2 ' 12 ' ) 



/ Vrms \ / U-5M \ f *-AU \ 

\5km/sJ V< m >/ V a& / 



Provided the derived encounter time-scale is shorter than the lifetime 
of the merger product (s), this suggests that the encounter scenario in 
question could be realistic and is therefore a candidate formation his- 
tory. Conversely, if the derived encounter time-scale is longer than the 
lifetime of the merger product (s), then the dynamical formation history 
is unlikely to have actually occurred. In general, the shorter the derived 
encounter time-scale, the more likely it is that one or more such encoun- 
ters actually took place within the lifetime of the merger product (s). 

The preceding guidelines specify a narrow range of allowed formation 
scenarios. In particular, they constrain the number of stars involved in 
the encounter (s), the types of objects (i.e. single, binary or triple stars) 
involved, and the number of encounters that took place. In most cases, 
these guidelines will converge on a single qualitative formation history 
that is unique up to the possible initial distribution of energies that 
describe suitable interactions. 
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2. Next, we must assign an approximate value based on the observations to 
every parameter possible in Equation 12.11 Nearly all of the required in- 
formation pertaining to the final distribution of energies in Equation 12. II 
can be found from spectroscopy alone, though repeated measurements 
spread out over a sufficiently long time-line will typically be required to 
obtain orbital solutions and to detect outer triple companions whenever 
they are present. This gives the final orbital energies as well as the 
gravitational binding energies Wu of the stars according to Equation 12.71 



and Equation 12. 4[ respectively. Since merger prod ucts have 



Deen shown 



Sills et al. 



20011), the 



to typically be in hydrodynamic equilibrium (e.g. 
stars' internal energies Uu can then be approximated using the virial 
theorem. Alternatively, stellar models can be used together with pho- 
tometry. 

The measured broadening of spectral lines gives an estimate for the stars' 
rotation speeds (although there is a strong dependence on the angle of 
inclination of the stars' axis of rotation relative to the line of sight), 
which in turn provides their rotational kinetic energies flu according to 
Equation 12.21 In conjunction with proper motions, radial velocity mea- 
surements also provide an estimate for the systemic velocity of the final 
stellar configuration relative to the cluster mean, which in turn gives 
its translational kinetic energy Tjj according to Equation 12.51 Finally, 
the cluster velocity dispersion provides an estimate of the relative veloc- 
ity at infinity for a typical encounter, which in turn decides the initial 
translational kinetic energies Tj of stars or binaries involved in the en- 
counter. Under the assumption of energy equipartition for both single 
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and binary stars, the initial velocities of the impactors can be approxi- 
mated by Equation 12.61 The assumption of energy equipartition should 
be valid in clusters for which the half-mass relaxation time is consid- 
erably sh^rterjlmn_jlie_du^ case for most 



GCs ([Harris 



1996. 2010 update 



De Angeli et al. 



20051 ). On the other 



hand, this assumption is likely invalid for most open clusters. This is 
because the tidal truncation of the velocity distributions are significant 
in OCs since they are much less centrally concentrated. In this case, the 
initial velocities of the impactors can be approximated from the velocity 
dispersion, which is nearly independent of mass. 

3. We can now obtain quantitative constraints for the initial encounter(s). 
Once we have decided on a qualitative encounter scenario that could have 
produced the observed merger product (s), we can estimate the orbital 
energies of the initial binaries or triples going into the encounter using 
Equation 12.91 This gives us an equation that relates the initial orbital 
semi-major axes of all orbits going into the encounter (to each other). If 
only one orbit goes into the encounter, then we can solve for it explicitly. 
From this, we can constrain the initial collisional cross sections for real- 
istic encounters. If the derived cross section is significantly smaller than 
the average semi-major axes of all binaries and/or triples in the cluster, 
then we can infer that the time required for the encounter to occur is sig- 
nificantly longer than the corresponding time-scale given in Appendix [Al 
If we use our derived cross section found from Equation 12.91 instead of 
this average semi-major axis, this should give us an idea of the time 
required for that particular type of encounter to occur. Strictly speak- 
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ing, however, this is only a rough approximation since these are typical 
time-scales found using the average period (or, equivalently, semi-major 
axis and hence cross section). 

If we find that the derived encounter time-scale is longer than the life- 
time of a typical merger product, then the chosen formation history is 
unlikely to have actually occurred. In this case, it becomes necessary to 
re-evaluate the possible dynamical formation histories of the observed bi- 
nary or triple system, choosing the next most likely qualitative scenario 
for further quantitative analysis. These steps can be repeated until ei- 
ther a suitable formation history is found or the list of possibilities is 
exhausted so that the only remaining conclusion is that the observed 
binary or triple system is unlikely to have a dynamical origin. 



2.3 Results 



In this section, we apply our technique to two particular cases of observed 
binaries and triples containing merger products. The first is an observed 
triple system in the old open cluster M67 that is thought to contain two BSs 



( Ivan den Berg et al. 



2001 



Sandquist et al. 



20031 ) . The second is the period- 



eccentricity distribution of the BS binary pop ulation of the OC NGC 188 



20091 ). Af- 



which bears a remarkable resemblance to M67 flMathieu fc Gellerl . 
ter determining the most probable qualitative formation scenarios, we obtain 
quantitative constraints for suitable initial conditions that could have produced 
the observed orbital parameters. 
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2.3.1 The Case of S1082 



S108 2 is believed to be a trip le system in the old OC M67 (Ivan den Berg et al. 



2001 



Sandquist et al 



20031 ). The observations suggest that a distant triple 
companion orbits a close binary containing a BS and another peculiar star. 
The companion to the BS has a photometric appearance that puts it close 
to the MSTO in the CMD and yet, curiously, its derived mass is significantly 
greater than that of the turn-off. The outer companion is a BS in its own right, 
so that S1082 is thought to be composed of two BSs. Although both the inner 
and outer components of this suspected triple have systemic velocities that 
suggest they are both cluster members, it is important to not e that there is no 



direc t evidence proving a dynamical link between the two (ISandquist et al. 



20031 ). Assuming for the time being that a dynamical link does exist, we can 



apply the procedure outlined in Section [2.2.21 to the case of S1082: 

1. Before applying our technique, it is important to convince ourselves that 
a dynamical origin is possible for the observed system. This is certainly 
the case for S1082 since no known BS formation mechanism could have 
produced the observed stellar configuration without at least some help 
from dynamical interactions. The first step of our procedure is to find 
qualitative constraints and, in so doing, isolate the most probable en- 
counter scenario. 

First, we need to know the cluster binary and triple fractions in order 
to use Figure 12.11 to find the encounter type occurring with the greatest 
frequenc y. From th i s, we find that 2+2 encounters presently dominate 



in M67. 



Fan et al. 



(119961 ) showed that observations of M67 are consis- 



tent with a cluster multiple star fraction ~ 50%. More recent studies 
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report a 



Latham 



ower 



2005 



imit for fh + ft that is consist ent with their results (e.g. 



Davenport fc Sandquist 



2010l ). Radial velocity surveys 



and simulations of dynamical interactions suggest that old OCs like M67 
are likely to host a number of triples with outer periods < 4000 days (e.g. 



Latham 



2005 



the result of 



Ivanova et a 



FanetaL 



20081 ) . We assume ft/ fb ~ 0.1 and, using 



(11996|), this gives f b ~ 0.45 and f t ~ 0.05. Our as- 



sumed ratio ft/ fh is slightly lo wer than found for the field, or ft/ fb ~ 0.2 



( lEggleton fc Tokovinin 



2008|) 



Second, we need to constrain the number of stars that went into the 



encounter. The mass o 



(IMathieu &: Geller 



( iSandquist et al. 



the MSTO in M67 is estimated to be ~ 1.3 



2009|). The total mass of S1082 is ~ 5.8 M G 







20031 ) . From Equation 12.111 its formation must there- 



fore have involved at least 5 stars. 

At this point, we can conclude that a single 2+2 encounter could not 
have produced S1082 since this scenario involves only 4 stars and we 
know that at least 5 stars are needed. We must therefore consider ei- 
ther a single 2+3 or 3+3 encounter, or a multiple encounter scenario. 
In order to isolate the most probable of these possibilities, we must cal- 
culate and compare their encounter time-scales. To do this, we require 
estimates from the observations for a few additional cluster parameters . 



For M67, the core radius is r, 



Giersz. Heggie fc Hurley 



1.23 pc (IBonatto. Bica Sz Santosl . 



2005 



20081 ). From this and the central velocity dis- 



p ersion, we can calcula t e the central mass density using Equation 4-124b 



of 



Binney &: Tremaind ( 119871 ). which gives po ~ 10 19 M pc 3 . The cen- 
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tral stellar number density can then be approximated according to: 



Po M 
< m > L ' 



(2.13) 



where M/L is the cluster mas s-to-light ratio and s hould be around 1.5 



for an OC as old as M67 (e.g. 



De Grijs et al. 



2008|). We take f t ~ 0.05 



and assume most stable triples have outer periods of P ~ 1000 days so 
that a t ~ 3 AU (assuming all three stars have a mass of 1 M Q ). We 
also take fb ~ 0.45 and assume an average binary period of P ~ 100 
days so that ~ 0.6 AU (assuming both components have a mass of 
1 M Q ). We assume an average stellar mass of < m >^ 1.0 M^ , whic h 



is in reasonable agreement with the observations (IGirard et al. 



19891 ). 



Assuming that the average ma ss of merger remn ants is equal to 2 < m > 



and extrapolating the results of 



Sills et al. 



(120011 ) for solar metallicity and 



more massive parent stars, we will assume that the typical lifetime of a 
merger product is r BS ~ 1.5 Gyrs. 

A comparison of the relevant encounter time-scales suggests that the 
most probable dynamical formation scenario for S1082 is a single 2+3 
encounter, although a single 3+3 encounter is almost equally as probable. 
Given our assumptions, we find 8.9 x 10 8 years and 3.3 x 10 9 years for 
the times between 2+3 and 3+3 encounters, respectively, in the core of 
M67. From this, we expect approximately two and zero 2+3 and 3+3 
encounters, respectively, to have occurred within the last tbs years. 

From Equation 12.121 we find th.&t the time for a particular binary to 
encounter another binary is 7.3 x 10 10 years. Similarly, the time for a 
particular quadruple to encounter another binary is 8.9 x 10 9 years (using 
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Equation 12. 12l and assuming the quadruple has a mass of 4 < m > and its 
geometric cross section is twice as large as the average outer semi-major 
axis of triples). Since these time-scales are considerably longer than 
the cluster lifetime, this suggests that a formation scenario for S1082 
involving back-to-back 2+2 encounters is unlikely, even if the second 
encounter occurred sufficiently soon after the first that all four stars 
comprising the initial pair of interacting binaries are still gravitationally 
bound. If we replace one of the 3 binaries involved in this scenario with 
a triple system, the total encounter time remains longer than the cluster 
lifetime. 



2. Before more quantitative constraints can be found, we must refer to the 
literature in order to obtain estimates for every term in Equation 12.91 
The observations suggest that a binary system composed of a 2.52 M 
blue straggler (component Aa) an d a 1.58 MS star (comp onent Ab) 



has a period of P ~ 1.068 days ( 



van den Berg et al. 



200 lh . There is 



also evidence for a 1.7 M blue straggler companion (component B) 
that acts as a stable o uter triple with a period of P ~ 1188.5 days 



( iSandquist et al. 



20031 ). From this, we can calculate the final orbital 



energies (e^ and e#) of the triple: 



+ \e B \ ~ \e A \ ~ 10 41 J, 



(2.14) 



where we have used the fact that je^l ^> |e#| since |e#| ~ 10 39 J. It is 
important to note that these peculiar stars are often under-luminous, so 
the inferred mass of the ter tiary companion should be taken with caution 



( Ivan den Berg et al. 



200lh . 
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The s ystemic radial velocity of the system is 33.76±0.12 km/s ( jSandquist et al. . 



20031 ). Although this only provides us with an estimate for the systemic 
velocity of S1082 in one dimension, it is consistent with the cluster mean 
velocity of 33.5 km/s. Therefore, the available knowledge for the sys- 
temic velocity of S1082 is consistent with its final translational kinetic 
energy being negligible. From this, we take Ta,r ~ 0- The cent r al ve 



locity dispersion in M 67 is only 0.5 km s 1 (IMathieu fc Latham 



1986 



Mathieu et al. 



19901 ) . This suggests that the relative velocity at infinity, 
and therefore the typical impact velocities of collisions, should be small 
compared to the significant orbital energy of component A in S1082 (and 
any very hard binaries that went into the encounter). From this, we take 
A m ~0. 

3. We are now ready to obtain quantitative constraints for the formation 
of S1082. We can use Equation 12.91 since M67 is an old OC with a 
small central velocity dispersion so that the assumptions used to derive 
this equation are valid. At this point, we must consider the details of 
our chosen formation scenario more carefully in order to choose a set of 
initial conditions that will satisfy Equation 12.91 as well as reproduce the 
observed parameters of S1082. In doing so, we find that it is not possible 
to form S1082 with a single 2+3 encounter. This is because 5 stars with 
masses m < rnxo are insufficient to form both an inner binary with a 
total mass ~ 4.1 M and an outer tertiary companion with a mass 1.7 
M . That is, the total mass of the inner binary is larger than three times 
the mass of the MSTO so that its formation must have involved four or 
more stars. 
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A single 3+3 encounter is therefore the most probable formation scenario 
for S1082. Although there are a number of ways we can reproduce the 
observed parameters of S1082 with a 3+3 encounter, including the com- 
ponent masses, we must use Equation 12.91 to identify the most probable 
of these scenarios. To do this, we can solve for the initial orbital ener- 
gies of all binaries and/or triples going into the encounter, which will in 
turn constrain their initial semi-major axes and therefore cross sections 
for collision. In applying Equation 12.91 we are only concerned with the 
initial and final orbital energies of the inner binaries of the triples. This 
is because the orbital energies of any outer triple companions will be 
significantly outweighed by the orbital energies of their inner binaries. 
It follows that the contribution from the outer orbital energies of triples 
to the total energy of the encounter will be negligible. 

From Equation 12. 9[ we find that the formation of S1082 should have 
involved at least one hard binary with |ej| ~ 10 41 J whose components 
did not merge (with each other) during the encounter in order to account 
for the significant orbital energy of component A. This need not be the 
case, however, if one or more stars were ejected from the system with an 
escape velocity of > 100 km/s. Of course, this would require an increase 
in the total number of stars involved in the interaction and therefore 
a single encounter scenario involving 7 or more stars, and hence one or 
more quadruple systems. Although very few observational constraints for 
the fraction of quadruple systems in clusters exist, this seems unlikely. 

The orbital energies of the inner binaries of the two triples initially going 
into the encounter should both be on the order of 10 41 J. Although we 
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have found that only one hard binary with |ej| ~ 10 41 J is required 
triples are only dynamically stable if the ratio of their inne r to outer 
orbit al periods is large (roughly a factor of ten or more) (e.g 



Mardlingl . 



20011 ). We do not expect very wide binaries to survive for very long in 
dense cluster environments, which suggests that the inner binary of every 
triple is very hard. We expect from this that these inner binaries should, 
to within an order of magnitude, all have roughly the same orbital energy. 

The presence of outer triple companions is required in order for most very 
hard binaries to actually experience encounters within the lifetime of a 
typical merger product. Assuming masses of 1 M Q for both components, 
an orbital energy of 10 41 J corresponds to a period ~ 0.4 days, or a 
semi-major axis ~ 0.02 AU. Therefore, the cross section for collision for 
a 2+2 encounter in which both binaries are very hard is smaller than 
the average cross section for a 2+2 encounter (found from the observed 
binary period distribution) by a factor of ~ 30. This suggests that the 
time required for an encounter to occur between two very hard binaries 
is considerably longer than the cluster lifetime. This is not the case if 
the hard binaries have triple companions, however, since the outer orbit 
significantly increases the collisional cross section and hence decreases 
the encounter time-scale. 

Energy conservation also suggests that if S1082 did form from a single 
3+3 encounter, it is likely that the close inner binaries of the two triples 
collided directly. If two separate mergers then subsequently occurred 
during this interaction, this could have reproduced the observed orbital 
parameters of the close inner binary of S1082 (component A). The for- 
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mation of the outer triple companion is more difficult to explain via a 
single 3+3 encounter since it also appears to be the product of a stellar 
merger. Nonetheless, if the outer companions of both triples undergoing 
the encounter end up orbiting the interacting (or merged) pair of close 
inner binaries at comparable distances, it is possible that their orbits 
would overlap as a result of gravitational focussing. Although this seems 
unlikely, it could produce a blue straggler of the right mass and orbital 
period to account for component B. We will return to this possibility in 
Section E2Q 

Even though the analytic estimates presented here are only approxi- 
mate, they serve to highlight the low probability of a system such as S1082 
having formed dynamically in M67. Based on our results, the most likely for- 
mation scenario for S1082 is a single 3+3 encounter, although we expect very 
few, if any, to have occurred in the lifetime of a typical merger product. This 
need not be the case, however, provided M67 had a higher central density in 
the recent past since this would increase the total encounter frequency. This 
will be discussed further in Section 12.41 

2.3.2 The Period-Eccentricity Distribution of the Blue 
Straggler Binary Population in NGC 188 



Mathieu fc Gellerl (120091 ) found 21 blue stragglers in the old open cluster 
NGC 188. Of these, 16 are known to have a binary companion. Orbital 
solutions have been found for 15 of these known BS binaries. From this, 



Mathieu fc Gellerl (I2009h showed that the BS binary population in NGC 188 



has a curious period-eccentricity distribution, with all but 3 having periods of 
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~ 1000 days. Of these three, two have periods of < 10 days (binaries 5078 
and 7782). Interestingly, one of these short-period BS binaries has a non-zero 
eccentricity. The normal MS binary populatio n, on the other hand, sho ws no 



2009|). We 



sign of a period gap for 10 < P < 1000 days (IMathieu fc Gellerl . 
can apply the procedure outlined in Section 12.2.21 to better understand how 
we expect mergers formed during dynamical encounters to contribute to the 
BS binary population in NGC 188. Although the method described in Sec- 
tion 12.2.21 treats one system at a time, we will apply our technique to the BS 
binary population of NGC 188 as a whole. 

1. Before applying our technique, we must satisfy ourselves that a dynam- 
ical origin is possible for a large fraction of the observed BS population. 
Several examples of evidence in favour of a dynamical origin exist. For 
one thin g, most BSs in NGC 188 have been found to have binary com- 



panions (IMathieu fc Geller 



20091 ). This should not be the case if most 
BSs are the products of the coalescence of isolated binary systems. On 
the other hand, a binary companion should be exp ected if the BSs were 



Chen fc Hani (I2008bh 



formed from mass transfer. However, the results of 
could suggest that most of the BSs in NGC 188 were not formed from 
mass transfer alone, providing indirect evidence that many BSs were 
formed from mergers. A binary companion should also be expected if 
the coalesced binary had a tertiary companion to begin with. Although 
the Kozai mechanism can act to decrease the orbital separation of the 
inner binary of a triple, an additional pe rturbation is often requir ed in 



order to fully induce the binary to merge (IPerets Sz Fabrycky 



20091 ). All 



of this suggests that dynamical interactions likely played at least some 
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role in the formation of a significant fraction of the BS binary population 
in NGC 188. 

The first step of our techni que is to fin d the most commonly occur- 
ring type(s) of encounter (s). 



Geller et al 



( 120081 ) found a completeness- 



corrected multiple star fraction fb + ft ~ 0.27 out to a period of ~ 4000 
days. This represents a lower limit since it does not include binaries 
with P > 4000 days. Us ing the same lower limit for the ratio ft/ fb ~ 0.1 
found by 



Latham! (]2005f ). Figure [2TTI suggests that 1+2 encounters should 
currently dominate in NGC 188. 

Second, we must constrain the minimum number of stars required to form 
the observed systems. From this, we can show that a merger occurring 
during all but a 1+1 encounter could produce a BS in a binary since this 
is the only type of encounter that involves less than 3 stars. By dividing 
the total mass of an observed system by rriTo, we can find an estimate 
for the minimum number of stars required to have been involved in its 
formation. Since all but one BS binary contain only a single merger 
product with a mass > rriTOi realistic dynamical formation scenarios for 
these systems require 3 or more stars. Binary 7782, on the other hand, is 
thought to contain two BSs so that if its formation involved two separate 
mergers it must have required 4 or more stars. 

A more quantitative comparison of the different encounter types is re- 
quired since we do not yet know if enough encounters occurred in the 
last tbs years to account for all 16 observed BS binaries. To do this, 
we require estimates from the observations for a few additional clus- 
ter parameters in order to calculate and compare the relative encounter 
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time-scales. NGC 188 has been 
IP 2 2 pc~ 3 (ISandauist et al. 



Po 



pc (IBonatto. Bica fc Santos 



bund to have a central density of 



20031 ) and a core radius of r c ~ 1.3 



20051 ). We found an average stellar mass 



for the cluster of < m >~ 0.9 M . This was done by determining the 



members from the proper-motion study oi_ 
this, a theoretical isochrone taken from 



of 


Platais et al. 


(2003) 


Pols et al. 


( 


1998 


) was 



determine the cluster mass function and the average stellar mass was 
calculated. 

A sufficient number of suitable dynamical interactions should have oc- 
curred in the last tbs years for the formation of every BS binary in NGC 
188 to have been directly mediated by the cluster dynamics. Assuming 
once again that tbs ~ 1.5 Gyrs, the encounter time-scales derived in 
Appendix [Al suggest that a minimum of nine, eight, eight, three and one 
1+2, 2+2, 1+3, 2+3 and 3+3 encounters, respectively, occurred within 
the lifetime of a typical merger product. It follows that at least ~ 29 
dynamical encounters should have occurred in NGC 188 in the last tbs 
years. Of these 29 encounters, 12 should have involved tr iples. In de 



riving these time-scales, we have taken /& + f t ~ 0.27 from 



Geller at al. 



( 120091 ) and have adopted the same ratio ft/ fb ~ 0.1 as found for M67. 
Consequently, these represent lower limits for fb and ft, so that our de- 
rived encounter rates are also lower limits. We have also assumed an 
average outer semi-major axis for triples of at ~ 3 AU (corresponding to 
a period of ~ 1000 days for a binary composed of two 1 M Q stars). 



2. The next step is to apply our energy conservation prescription to the 
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observed BS population. This will allow us to constrain the orbital 
energies of typical binaries or triples expected to form BS binaries via 
dynamical interactions. First, we must refer to the literature in order 
to obtain estimates for every term in Equation 12.91 From the observed 
BS period-eccentricity distribution in NGC 188, we know that most BS 
binaries have periods of ~ 1000 days (we will call these long-period 
binaries), although there are a couple with very short periods of ~ 10 
days (which we will call short-period binaries). This provides an estimate 
for the final orbital energies e^k of BS binaries formed during dynamical 
interactions. Specifically, we find \ckk\ ~ 10 39 J and je^] ~ 10 40 J for the 
final absolute orbital energies of the long- and short-period BS binaries, 
respectively. 

Every BS in NGC 188 with a high cluster membership probability has 
both radial and proper motion velocities that, to within their respective 
uncertainties, are co nsistent with the observed centra l velocity dispersion 



of On = 0.41 km s 1 (IPlatais et al. 



2003 



Geller et al. 



2008 



Geller at al. 



20091 ). From this, we assume that the final translational kinetic energies 
of any binaries or triples formed from dynamical interactions will be neg- 
ligible, or Tjj ~ 0. As for M67, we assume A m ~ since we expect low 
impact velocities for collisions as a result of the low velocity dispersion 
in NGC 188. 



3. We are now equipped with estimates from the literature that will allow 
us to obtain quantitative constraints for specific encounter scenarios. 
In particular, we can use Equation 12.91 to constrain the initial orbital 
energies of all binaries and/or triples going into an encounter. We can 
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also constrain the specific details of interactions for which the encounter 
outcome reproduces the observed parameters of the BS binary. We will 
consider energy conservation separately for two different classes of BS 
binaries, namely short- and long-period. 

The short-period BS binaries have large (absolute) orbital energies. Equa- 
tion 12.91 suggests that this energy requires that at least one hard binary 
was involved in the encounter. Alternatively, an encounter which in- 
volved only softer binaries must have resulted in at least one star es- 
caping with a significant velocity (> 100 km/s). However, in order to 
produce a merger product in a short-period binary, encounters involving 
triples are the most favoured. Binary-binary (and especially 1+2) en- 
counters which involve at least one hard binary will have smaller cross 
sections for collision than those encounters involving wide binaries. But 
energy conservation suggests that 2+2 encounters involving wide bina- 
ries are unlikely to produce a merger product in a clos e binary. Stable 



20011 ) 



triples, on the other hand, contain a hard inner binary (IMardling 
which will naturally account for the large orbital energy of the resulting 
BS binary. Stable triples also have a large cross section for collision be- 
cause of the wide orbit of the tertiary companion. The times between 
1+2, 2+2 and 1+3 encounters are all comparable, suggesting that most 
encounters involving hard binaries are 1+3 encounters. 

The short-period BS binaries could have formed from a direct stellar 
collision that occurred within a dynamical encounter of a hard binary 
and another single or (hard) binary star (we call this Mechanism I). If at 
least one of the objects going into the encounter was a triple, then four or 
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more stars were involved in the interaction. Therefore, if binaries 5078 
and 7782 were formed from this mechanism, they could possess triple 
companions with sufficiently long periods that they would have thus far 
evaded detection. This is consistent with the requirements for both con- 
servation of energy and angular momentum. Interestingly, the presence 
of an outer triple companion could also contribute to hardening these 
BS binaries via Kozai cycles o perating in conjunction with tidal friction 



Fabrycky fc Tremaine 



20071 ). In this case, Equation 12.11 shows that 
tides can contribute to making a binary's orbital energy more negative 
by depositing internal energy into the component stars. 

Although 2+2 encounters should also contribute to the observed BS pop- 
ulation, most of these should occur between a short-period binary and a 
long-period binary. This is because most encounters resulting in mergers 
involve very hard binaries and binaries with long periods have large cross 
sections for collision and therefore short encounter times. Equation 12.91 
indicates that the minimum period of a BS binary formed during a 2+2 
encounter is usually determined by the orbital energy of the softest bi- 
nary going into the encounter. This is because the most likely merger 
scenario is one for which the hard binary is driven to coalesce by im- 
parting ene rgy and angular mom entum to other stars involved in the 



interaction ( iFregeau et al 



20041 ). Assuming most of this energy is im- 
parted to only one of the stars causing it to escape the system, we can 
take A m in Equation 12.91 to be very small and the orbital energy of the 
left-over BS binary will be comparable to the orbital energy of the initial 
wide binary going into the encounter. As more energy is imparted to the 
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star left bound to the merger product over the course of the interaction, 
its final orbital separation effectively increases. This in turn contributes 
to an increase in the final orbital energy of the left-over BS binary. Fi- 
nally, from Equation \A.8\ the period of the initial wide binary should 
be relatively long (roughly > 1000 days) since the contribution from the 
very hard binary to the total cross section for collision is negligible. Oth- 
erwise, the time required for such a 2+2 encounter to occur could exceed 

TBS- 

Now let us consider the long-period BS binaries in NGC 188. Based on 
Equation 12.91 we expect encounters involving 3 or more stars and only 
one very hard binary to typically produce long-period BS binaries if the 
hard binary is driven to merge during the encounter. This is because the 
hard binary merges so that its significant (negative) orbital energy can 
only be re-distributed to the other stars by giving them positive energy. 
Since the only other orbits involved in the interaction are wide, the left- 
over BS binary should also have a wide orbit. Alternatively, BSs formed 
during interactions involving more than one very hard binary should be 
left in a wide binary provided enough energy is extracted from the orbit 
of the binary that merges. In this case, a significant fraction of this 
energy must be imparted to the other stars in order to counter-balance 
the significant orbital energy of the other very hard binary. If not, the 
other short-period binaries are required to either merge or be ejected 
from the system. Otherwise Equation 12.91 suggests that the left-over BS 
binary should be very hard. However, it is important to recall that we 
are neglecting other non-dynamical mechanisms for energy extraction. 
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We will return to this last point in Section 12.41 

In order to help us obtain more quantitative estimates for the long-period 
BS binaries, consider two additional mechanisms for mergers during dy- 
namical interactions that involve both short- and long-period orbits. 
First, a merger can occur if a sufficient amount of orbital energy is ex- 
tracted from a hard binary by other interacting stars (called Mechanism 
Ha). Alternatively, a merger can occur if the encounter progresses in 
such a way that the eccentricity of a hard binary becomes sufficiently 
increased that the stellar radii overlap, causing the stars to collide and 
merge (called Mechanism lib). In the case of Mechanism lib, most of the 
orbital energy of the close binary must end up in the form of internal and 
gravitational binding energy in the merger remnant after the majority of 
its orbital angular momentum has been redistributed to the other stars 
(and tides have extracted orbital energy). In the case of Mechanism Ha, 
however, most of the orbital energy of the close binary must be imparted 
to one or more of the other interacting stars in the form of bulk kinetic 
motion. Consequently, one or more stars are likely to obtain a positive 
total energy and escape the system. This need not necessarily be the 
case for 2+3 and 3+3 encounters provided the second hardest binary 
orbit involved in the interaction has a sufficiently negative energy. 

Regardless of the type of encounter, Equation 12.91 shows that the extrac- 
tion of orbital energy from a hard binary in stimulating it to merge should 
increase the final orbital period of a BS binary. To illustrate this, we will 
consider 1+3 encounters since the predictions of energy conservation are 
nearly identical for the other encounter types of interest. Moreover, we 
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have shown that encounters with triples are the most likely to involve 
very hard binaries. Equation 12.91 can be re- written for 1+3 encounters: 

£12,4 — e i2,3 — fi2 x e i2 ; (2-15) 

where we have assumed that stars 1 and 2 comprise the initial hard inner 
binary of the triple, star 3 is the initial outer triple companion and star 4 
is the interloping single star. Stars 1 and 2 are assumed to merge during 
the encounter by exchanging energy and angular momentum with stars 
3 and 4. We further assume that enough energy is imparted to star 3 
that it escapes the system. We let / 12 represent the fraction of energy 
extracted from the orbital energy of the hard inner binary of the triple 
in the form of bulk kinetic motion by star 4. Since the remaining orbital 
energy of the close inner binary will end up in the form of internal and 
gravitational binding energy in the merger remnant, we have assumed 
(1 + fn) x ^12 ~ T 3 — A m in obtaining Equation 12.151 from Equation 12.91 

This formation mechanism could leave the BS as a single star if /12 
exceeds a few percent. As shown in Figure 12.21 the period of a BS 
binary formed during a 1+3 encounter via Mechanism Ha {P\2,±) is only 
slightly smaller than the period of the outer orbit of the triple initially 
going into the encounter (^12,3)- This assumes, however, that no energy 
is exchanged between the hard inner binary and star 4 (i.e. /12 = 0). The 
predictions from energy conservation for this case are therefore identical 
for Mechanism lib. In general, as the amount of energy extracted from 
the hard inner binary of the triple by star 4 increases, so too does the rate 
at which Pi2,4 increases with increasing Pi 2 ,3- If the amount of energy 
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extracted is > 5% of the orbital energy of the initial inner binary of the 
triple, P\2,a becomes a very steeply increasing function of Pi2,3- 




2000 4000 6000 8000 

Initial Outer Period (in days) 



Figure 2.2 Plot showing the typical periods of BS binaries formed from 1+3 
encounters in which the hard inner binary of the triple merges. As described in 
the text, the period of the BS binary formed during the interaction is denoted 
P\2A and corresponds to the y-axis, whereas the period of the outer orbit of the 
triple initially going into the encounter is denoted Pi.2,3 and corresponds to the 
x-axis. The lower straight line corresponds to the case where no energy was 
extracted from the inner binary of the triple by star 4 (i.e. fi2 = 0). As the 
amount of energy extracted increases, however, so too will the rate at which 
Pi2,4 increases with increasing Pvz,z- Cases where /12 = 0.005, /12 = 0.01, 
/12 = 0.02, f u = 0.03 and / i2 = 0.04 are shown as lines of increasing slope. 



Interestingly, the two general qualitative merger scenarios described above 
(Mechanisms I and II) naturally create a bi-modal period distribution 
similar to the period gap observed for the BS binaries if we assume that 
1+3 encounters produced these objects. To illustrate this, Figure 12.31 



76 



Ph.D. Thesis — N. Leigh 



McMaster - Physics & Astronomy 



shows a histogram of periods for 15 BS binaries formed during 1+3 en- 
counters via these two generic merger scenarios. In order to obtain this 
plot, we have used the observed period-eccentricity distribution f or the 



Geller at al. 



(120091 ) 



regular MS-MS binary population in NGC 188 from 
to obtain estimates for the orbital energies of any binaries and/or triples 
going into encounters. Specifically, in order to obtain periods for the 
outer orbits of triples undergoing encounters, we randomly sampled the 
regular period distribution, including only those binaries with periods 
satisfying 400 days < P < 4000 days. We have shown that the initial 
outer orbits of triples going into 1+3 encounters provide a rough lower 
limit for the periods of BS binaries formed via Mechanism II. Therefore, 
any BS binaries formed in this way could only have been identified as 
binaries by radial velocity surveys if the triple going into the encounter 
had a period < 4000 days (since this corresponds to the current cut-off 
for detection). All triples are taken to have a ratio of 30 between their 
inner and outer orbital semi-major axes. This ratio has been chosen to 
be arbitrarily large enough that the triples should be dynamically stable, 
however we will return to this assumption in Section 12.41 

We w ill adopt a ratio based on the observations of 



Mathieu fc Geller 



(120091 ) for the fraction of outcomes that result in each of these two pos- 
sible merger scenarios. In particular, if we assume that the three BS 
binaries with P < 150 days were formed via Mechanism I whereas the 
other 12 were formed via Mechanism II (either Ha or lib since energy 
conservation predicts similar periods for the left-over BS binaries), this 
would suggest that Mechanism II is ~ 4 times more likely to occur than 
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Mechanism I during any given 1+3 encounter. We will return to this 
assumption in Section I2~4l 



x 



y v v v 



X X X X X 



log P (in days) 

Figure 2.3 Histogram of the period distribution (in days) expected for BS 
binaries formed during 15 1+3 encounters. The parameter space assumed for 
the encounters is described in the text. The s tars show the observed B S binary 
period distribution in NGC 188 taken from iMathieu fc Gellerl ( 120091 ). where 
each star represents a single BS binary. 



As a result of the requirement for conservation of angular momentum, 
we would expect to see a large spread in the distribution of eccentricities 
for BSs in long-period binaries formed from encounters involving triples. 
This is because conservation of momentum requires that the total mo- 
mentum going into the encounter must be equal to the total momentum 
left over after the interaction. However, the total initial momentum 
depends not only on the initial orbital eccentricities of any binaries or 
triples going into the encounter, but also the relative orientations and 
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trajectories of the colliding objects. Since the relative orientations and 
trajectories are random, the final eccentricities of the BS binaries can 
take on a range of values. In other words, we cannot predict the final 
distribution of momenta for BS binaries formed from dynamical encoun- 
ters. However, the observed BS binaries in NGC 188 are observed to 
have a range of eccentricities and this is not inconsistent with a dynam- 
ical origin involving triples. 

BSs in short-period binaries formed via Mechanism I can also end up 
with just about any eccentricity immediately after the encounter for the 
same reasons outlined above. However, tidal effects become increasingly 
significant with decreasing orbital separation so that the hardest binaries 
should typically circularize the fastest. Althoug h theoretical es t imate s 



for the rate of tidal circularization are uncertain (IMeibom et al. 



2005|) 



the c ircularization cut-off period is estimated to be ~ 15 days in NGC 



188 flMathieu et al 



20041 ). From this, it is entirely possible that re- 
cently formed BS binaries with P ~ 10 days have not yet become fully 
circularized. 



2.4 Summary & Discussion 

In this chapter, we have presented a generalized analytic prescription for energy 
conservation during stellar encounters. Our method can be used to identify 
the most probable dynamical formation scenario for an observed binary or 
triple system containing one or more merger products. We have shown that, 
using the observed orbital parameters of the system, the allowed initial orbital 
semi-major axes of any binary or triple systems involved in its formation can 
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be constrained. The initial semi-major axes of the orbits in turn provide 
an estimate for the collisional cross section and therefore the time-scale for 
the encounter to occur in its host cluster. In order to apply our technique, 
repeated spectroscopic measurements of the binary or triple system containing 
the merger product (s) are needed in order to obtain its orbital solution and 
systemic velocity. However, the time-scales provided in Appendix |A] can still 
be applied if only the fraction of binary and triple stars are known, which 



can 



3e determined either spect roscopica 



20051 ) or photometrically (e.g. 



FanetaL 



l v (e.g . 



Mathieu et al. 



1990 



Latham 



19961 ). 



As we have shown, consideration of the requirement for energy conser- 
vation is ideal for identifying trends during stellar encounters, whereas numer- 
ical scattering experiments can require hundreds or even thousands of simula- 
tion runs before patterns will emerge. Some of these trends include: 



At least one short-period binary is usually required in a dynamical inter- 
action to produce another binary having a similarly short-period (pro- 
vided no stars are ejected with escape velocities > 100 km/s). This is 
because the orbital energy of a short-period binary is sufficiently neg- 
ative that it tends to considerably outweigh the other energy terms in 
Equation 12.11 for most of the encounters that typically occur in the cores 
of globular and esp ecially open clusters. This has been confirmed by 



Hurley et al. 



(|2005|). 



Previous studies have found that in order for triples to remain stable for 
many dynamical times, the ratio of their inner to outer orbit al periods 



must be relatively large (roughly a factor of ten or more) (e.g. 



Mardling 



200ll ). Based on our results, this has two important corollaries for stellar 
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mergers in dense cluster environments hosting a significant population 
of triples: 

1. The longest- lived triples will contain very hard inner binaries with 
a large |e|. This is important since stellar radii are in general more 
likely to overlap and hence mergers to occu r during resonant in 



teractions involvin g very hard binaries (e.g. 



Hurley et al. 



Fregeau et al. 



2004 



20051 ). We therefore expect stellar mergers to be com- 
mon during encounters involving stable triple systems. 

2. The longest-lived triples will contain wide outer orbits, creating a 
large cross section for collision. This suggests that the time-scale 
required for a stable triple system to encounter another object is 
typically short relative to the cluster age in dense environments. A 
significant fraction of encounters involving very hard binaries, and 
hence resulting in stellar mergers, will therefore involve triples in 
old open clusters such as M67 and NGC 188. 



van den Berg et al 



( 1200 lh and 



Sandquist et al. 



(120031 ) suggest that back- 



to-back binary-binary encounters, or even a single 3+3 encounter, could have 
formed S1082. We have improved upon these previous studies by estimating 
time-scales required for possible dynamical formation scenarios to occur. Since 
we have argued in Section l2~3l that the formation of S1082 must have involved 
at least 6 stars, it follows that only a 3+3 interaction could have reproduced 
the observed configuration via a single encounter. However, the derived 3+3 
encounter time-scale is sufficiently long that we expect very few, if any, 3+3 
encounters occurred within the lifetime of a typical merger product. Moreover, 
we have argued that the times for multiple encounters to occur are longer than 
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the cluster age. Although we cannot rule out a dynamical origin for S1082, our 
results suggest that it is unlikely (provided the derived encounter time-scales 
were not significantly higher in the recent past, which we will return to below). 
From this, it follows that a dynamical link between the close binary and third 
star is unlikely to exist. 

On the other hand, we have so far ignored the cluster evolution, and 
assumed that the currently observed cluster parameters have not changed in 
the last few Gyrs. N-body models suggest that the central density in M67 
co uld have been signi ficantly higher in the recent past. Specifically, Figure 5 



of 



Hurley et al. 



( 120051 ) indicates that the presently observed central density in 
M67 could have been higher within the past Gyr by a factor of > 2. If this 
was indeed the case, our previous estimates for each of the different encounter 
frequencies should increase by a factor of ~ 4, so that a significant number 
of dynamical encounters involving single, binary and triple stars should have 
occurred in M67 within the last tbs years. It follows that a dynamical origin for 
S1082 is not unlikely if the central density in M67 was recently larger than its 
presently observed value by a factor of > 2. This also increases the probability 
that a scenario involving multiple encounters created S1082, although we have 
shown that such a scenario is still likely to have involved one or more triples. 

Based on the preceding arguments, SI 082 offers an excellent example 
of how observations of individual multiple star systems containing BSs can 
be used to directly constrain the dynamical history of their host cluster. If 
a definitive dynamical link between components A and B is established, this 
would suggest that the central density in M67 was higher in the last 1-2 Gyrs. 
This is also required in order for the cluster dynamics to have played a role in 
the formation of a significant fraction of the observed BS population in M67. 
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Based on the current density, the encounter time-scales are sufficiently long 
that too few encounters should have occurred in the last tbs years for mergers 
during dynamical interactions to be a significant contributor to BS formation. 

We have obtained quantitative constraints for two generic channels for 
mergers during encounters involving triples - one in which a direct stellar 
collision occurs within a dynamical interaction of the hard inner binary of 
a triple and another single or (hard) binary star (Mechanism I) and one in 
which the hard inner binary of a triple is driven to coalesce by imparting 
energy and/or angular momentum to other stars involved in the interaction 
(Mechanism II). Our results suggest that these two general merger mechanisms 
could contribute to a bi-modal period distribution for BS binaries similar to 
that observed in NGC 188. These dual mechanisms predict a gap in period, 
with those BS binaries formed via Mechanism I having periods of a few to 
~ 100 days and those formed via Mechanism II having periods closer to ~ 1000 
days. Some 2+2, 2+3 and even 1+3 encounters could involve orbits with 
periods in this range, and Equation 12.91 confirms that the final period of a BS 
binary formed via Mechanism II will typically be determined by that of the 
second hardest binary orbit. Therefore, we might still expect some BS binaries 
to have periods that fall in the gap (100 days < P < 1000 days). Our results 
do indeed predict one such BS binary, as shown in Figure 12.31 

A number of assumptions went into obtaining Figure I2~3l many of which 
were chosen specifically to reproduce the observed BS binary period distribu- 
tion. Regardless, our assumptions were chosen to reflect encounter scenarios 
that are the most likely to result in mergers. These should involve triples with 
very hard inner binaries since these are the most likely to merge during en- 
counters. The triples should also have outer companions on very wide orbits 
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since these have the largest cross sections for collision. From this, we have as- 
sumed a ratio of 30 between the inner and outer semi-major axes of all triples. 
This ensures that all triples are dynamically stable. This also leads us to as- 
sume a minimum period of 400 days for the outer orbits of triples so that the 
corresponding minimum period of their inner orbits is not too small. These as- 
sumptions serve to show that encounters involving triples could produce both 
long-period and short-period BS binaries as well as a period gap. 

Our results predict that the short-period peak in Figure 12.31 is at a 
slightly longer period than the observations suggest. If we decrease the as- 
sumed ratio between the inner and outer orbital separations of triples, the 
short-period peak will move to even longer periods. However, if one or more 
stars were ejected with a high escape velocity or the dissipative effects of tides 
are considered (which are expected to be the most significant for encounters 
involving hard binaries) this would move the short-period peak to even shorter 
periods. In order to obtain the desired agreement with the observations at the 
short-period end of the BS period distribution, energy must have somehow 
been dissipated or removed from the hard inner binaries of the triples during 
(or even after) the encounter, or the inner binaries must have initially been 
even harder than we have assumed in obtaining Figure 12.31 

Our results could suggest that the hard BS binaries in NGC 188 (bina- 
ries 5078 and 7782) may have outer triple companions, perhaps with sufficiently 
long periods that they would have thus far evaded detection. This is consis- 
tent with the requirement for energy conservation since the orbital energy of 
the outer orbit of the triple is negligible compared to that of its inner binary. 
If binaries 5078 and 7782 do have outer triple companions, it is also possible 
that Kozai oscillations combined with tidal friction contributed to decreasing 
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their orbital periods ( lEggletonl . 120061 ) . Finally, the BS binaries in NGC 188 
are observed to have a wide range of eccentricities, which we have argued is 
not inconsistent with a dynamical origin involving triples. 

We have assumed that Mechanism II is more likely to occur than Mech- 
anism I during any given 1+3 encounter. This is a reasonable assumption since 
numerical scattering experiments of 1+2 and 2+2 encounters have shown that 
the coalescence of a hard binary is much more likely to occur during an en- 
counter than a dire ct collision between a hard binary and an interacting single 



or binary star (e.g. 



Fregeau et al. 



20041 ) . The ratio we have assumed for the 



frequencies with which Mechanisms I and II occur was specifically chosen in 
order to reproduce the observed numbers of short- and long-period BS bi- 
naries. The important point to take away is that the observed BS binary 
period-eccentricity distribution offers a potential constraint on the fraction of 
encounters that result in different merger scenarios. 

Based on our results, Mechanism II must occur ~ 4 times more often 
than Mechanism I in order to reproduce the observed BS period distribution 
from 1+3 encounters (or, equivalently, 2+3 encounters involving a very wide 
binary and 2+2 encounters between a short-period binary and a long-period 
binary). This can be tested by performing numerical scattering experiments 
of encounters involving triples. Therefore, our results highlight the need for 
simulations of 1+3, 2+3 and 3+3 encounters to be performed in order to 
better understand their expected contributions to BS populations in open and 
globular clusters. Once a preferred encounter scenario has been identified for 
an observed binary or triple containing one or more BSs, numerical scattering 
experiments can be used to further constrain the conditions under which that 
scenario will occur (or to show that it cannot occur). We have demonstrated 
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that a combination of observational and analytic constraints can be used to 
isolate the parameter space relevant to the dynamical formation of an observed 
multiple star system (or population of star systems) containing one or more 
merger products. This will drastically narrow the relevant parameter space 
for numerical scattering experiments. 



We have improve d upon the results of 



Perets fc Fabrvckv! (120091 ) and 



Mathieu fc Gellerl (]2009[ ) since we have shown that dynamical encounters in- 
volving triples could not only be contributing to the long-period BS binaries 
in NGC 188, but they could also be an important formation mechanism for 
short-period BS binaries and triples containing BSs. We have not ruled out 
mass transfer or Kozai-induced mergers in triples (primordial or otherwise) 



( Mathieu fc Geller 



2009 



Perets fc Fabrvckv . 



20091 ). or even various combina- 



tions of different mechanisms, as contributing formation channels to the BS 
binary population in NGC 188. For instance, a 1+3 exchange interaction could 
stimulate a merger indirectly if the resulting angle of inclination between the 
inner and outer orbits of the triple exceeds ~ 39°, ultimately allowing the 
triple to evolve via the Kozai mechanism so that the eccent ricity of the inn er 



20061). 



binary increases while its period remains roughly constant ( lEggletonl . 

There is evidence to suggest that mass transfer via Roche lobe over flow 
could play a role in the formation of at least some BSs. It is difficult to account 
for the near zero eccentricities of some of the long-period BS binaries without 
at least one episode of mass transfer having occurred. This is because none of 
the normal MS-MS binar ies with similar periods have such small eccentricities 



( IMathieu fc Geller 



20091 ) . On the other hand, it may not be unreasonable 
to expect that some collision products left in binaries undergo mass transfer 
since they are expected to expand adiabatically post-collision, and will sooner 
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or later evolve to ascend the giant branch. As a result of conservation of energy 
and angular momentum, the mass transfer process will usually ac t to increas e 



the orbital periods of these binaries provided it is conservative (liben . 



199ll ). 



Interestingly, t he cut-off period for Roche lobe overflow is ~ 1000 days for 



low-mass stars (jEggleton 



20061 ). which is in rough agreement with the long- 
period peak in the observed period-eccentricity distribution of the BS binary 
population in NGC 188. Therefore, mass transfer could also be contributing 
to the period gap observed for th e BS binaries 



According to the results of 



Geller at al. 



(|2009|), the number of giant-MS 



binaries with P < 1000 days is comparable to the number of BS binaries (A. 
Geller, private communication). It is unlikely that every giant-MS binary will 
form a BS from mass transfer, however, suggesting that at most a few of the 
long-period BS binaries in NGC 188 were formed via this mechanism. Finally, 
if the outer companion of a triple system evolves to over-fill its Roche lobe it 
could transfer mass to both of the components of the close inner binary. This 
mechanism could therefore also produce two BSs in a close binary, although 
it predicts the presence of an orbiting triple companion. For these reasons, a 
better understanding of triple evolution, as well as binary evolution in binaries 
containing merger products, is needed. 

The dissipational effects of tides tend to convert stars' bulk transla- 
tional kinetic energies into internal or thermal energy within the stars, leading 
to an increas e in the total gravitationa l bindi ng energy of the stellar config- 



uration (e.g. 



McMillan. Hut fc Makino 



19901 ). By increasing the terms U, 



in Equation 12.11 the initial orbital energies of any binaries going into an en- 
counter can increase accordingly in order to conserve energy. A higher orbital 
energy corresponds to a larger semi-major axis and hence cross section for 
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collision. This suggests that the derived encounter time-scales can be taken as 
upper limits in the limit that tidal dissipation is negligible. We expect tides to 
be particularly effective during encounters for which the total energy is very 
negative as a result of one or more very hard binaries being involved. 

We have argued in Section 12.2.21 that the average stellar mass is ex- 
pected to be comparable to (but slightly less than) the mass of the MSTO in old 
OCs and low-mass GCs. We have also argued that most encounters will involve 
stars having masses slightly larger than the average stellar mass. We might 
therefore expect to find that a high proportion of merger products have masses 
that exceed that of the MSTO in very dynamically-evolved clusters that have 
lost a large fraction of their low-mass stars. Consequently, a larger number of 
merger products could appear sufficiently bright to end up in the BS region of 
the cluster CMD in these clusters than in their le ss dynamically-evolved coun- 



terp arts. This is consistent wit 
and 



Knigge. Leigh & Sills! (120091 ) 



i the results of 
09) who found that 
cores of GCs scales sub-linearly with the core mass. In particular, since the 
cluster relaxation time increases with increasing cluster mass, it is the lowest 
mass GCs that should have lost the largest fraction of their low-mass stars. 
Therefore, if a larger fraction of merger products do indeed end up more mas- 
sive than the MSTO in these clusters, this could be a contributing factor to 
the observed sub-linear dependence on core mass. It is also interesting to note 
that, since BSs are among the most massive cluster members and many are 
thought to have a binary companion, BSs should be preferentially retained 
in clusters as they evolve dynamically compared to low-mass MS stars. This 
could also contribute to the observed sub-linear dependence on core mass. 

The exchange or conversion of energies that occurs during an encounter 
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takes place over a finite period of time, so it is important to specify whether 
or not the system has fully relaxed post-enc ounter when disc ussing the re- 
maining stellar configuration. For one thing, 



Sills et al. 



( 120011 ) showed that, 



although collision products may be in hydrodynamic equilibrium, they are not 
in thermal equilibrium upon formation and so contract on a thermal time-scale. 
Simulations also s ugges t that most merger products should be rapid rotators 



flSills et al. 



2002 



20051 ). However, at least in the case of blue stragglers, this 
is rarely supported by the observations. Some mechanism for angular momen- 
tum loss must therefore be operating either during or after the merger takes 
place in order to spin down the remnant. The time-scale considered must also 
be sufficiently short that subsequent dynamical interactions are unlikely to 
have occurred since these could affect the total energy and momentum of the 
system. 

With this last point in mind, N-body simulations considering BS forma- 
tion have shown that after they are forme d, BSs are often exchanged into other 



multiple star systems (IHurley et al. 



20051 ) . This suggests that for multiple star 



systems containing more than one BS, the BSs could have first been formed 
separately or in parallel, and then exchanged into their presently observed 
configuration. For the case of S1082, this would require at least 3 separate 
dynamical interactions. Given that the derived times between encounters are 
relatively long and the fact that the most likely formation scenario is usually 
that for which the number of encounters is minimized, the current state of 
M67 suggests that the probability of S1082 having formed from a scenario 
involving 3 encounters is low. Conversely, the derived encounter time-scales 
in NGC 188 are sufficiently short that many of the BS binaries could have 
experienced a subsequent dynamical interaction after their formation. BSs 
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tend to be more massive than normal MS stars, contributing to an increase 
in their gravitationally-focussed cross section for collision. This suggests that 
the encounter time-scale for multiple star systems containing BSs is slightly 
shorter than for otherwise identical systems composed only of normal MS stars. 
This contributes to a slight increase in the probability that a BS will expe- 
rience an exchange encounter after it is formed. Interestingly, it could also 
contribute to an increase in the probability that a close binary containing two 
BSs will form during an encoun ter between two di f ferent multiple star systems 



each containing their own BSs ( iMathieu fc Gellerl (120091 ) ; R. Mathieu, private 



communication). This is because it is the heaviest stars that will experience 
the strongest gravitational focussing and are therefore the most likely to ex- 
perience a close encounter, end up in a closely bound configuration, or even 
merge. 

Most exchange interactions will involve wide binaries for which the 
cross section for collision is large. Since wide binaries are typically relatively 
soft and the hardest binary involved in the interaction will usually determine 
the orbital energy of the left-over BS binary, most exchange interactions will 
leave the periods of BS binaries relatively unaffected. This need not be the 
case, of course, provided one or more stars are ejected from the system with 
a very high escape velocity. With these last points in mind, we have assumed 
throughout our analysis that all binaries and triples are dynamically hard. 
This is a reasonable assumption since the hard-soft boundary corresponds to a 
period of ~ 10 6 days in both M67 and NGC 188, for which the cross section for 
collision is sufficiently large that we do not expect such binaries to survive for 
very long. Nonetheless, considerations such as these must be properly taken 
into account when isolating a preferred formation scenario and predicting the 
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final distribution of energies. 

We have presented an analytic technique to constrain the dynamical 
origins of multiple star systems containing one or more BSs. Our results 
suggest that, in old open clusters, most dynamical interactions resulting in 
mergers involve triple stars. If most triples are formed dynamically, this could 
suggest that many stellar mergers are the culmination of a hierarchical build- 
up of dynamical interactions. Consequently, this mechanism for BS formation 
should be properly included in future N-body simulations of cluster evolution. 
A better understanding of the interplay between the cluster dynamics and the 
internal evolution of triple systems is needed in order to better understand the 
expected period distribution of BS binaries formed from triples. Simulations 
will therefore need to track both the formation and destruction of triples as 
well as their internal evolution via Kozai cycles, stellar and binary evolution, 
etc. On the observational front, our results highlight the need for a more 
detailed knowledge of binary and especially triple populations in clusters. 
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Chapter 



3 



Stellar Populations in Globular 
Cluster Cores: Evidence for a 
Peculiar Trend Among Red Giant 
Branch Stars 



3.1 Introduction 

Studying the radial distributions of the various stellar populations (red giant 
branch, horizontal branch, main-sequence, etc.) found in globular clusters 
(GCs) can provide useful hints regarding their dynamical histories. As clus- 
ters evolve, they are expected to undergo relatively rapid mass stratification 



as a consequence of two-body rela xation, with t 



ing to the central cluster regions (ISpitzer 



re heaviest stars quickly sink- 



19871 ). The shorter the relaxation 



time (typically evaluated at the half- mass radius), the quicker this process oc- 
curs. Clusters tend towards dissolution as two-body relaxation progresses and 
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they lose mass due to stellar evolution and the preferential escape of low-mass 
stars. External effects like tidal perturbations, encounters with giant molecular 
clouds, a nd passages through the Galactic disk serve only to speed up the pro- 



cess (e.g. 



Baumgardt &: Makino 



2003 



Kiipper et al. 



20081 ). Stellar evolution 



complicates this otherwise simple picture of GC evolution, however. Stars are 
expected to change in size and lose mass as they evolve, often dramatically, and 
this could significantly impact the outcomes of future dynamical interactions 
with other stars. For instance, a typical star in a GC is expected to expand 
by up to a few orders of magnitude as it ascends the red giant branch (RGB) 
and will shed up to a quarter of its mas s upon evolving from the tip of the 



RGB to the horizontal-branch (HB) (e.g. 



Caloi fc D'Antona 



2008 



Lee et al. 



1994h . 



Red giant branch stars have been report ed to b e defic ient in the cores 



Bailvnl (119941 ) found that the 



of some Milky Way (MW) GCs. For instance, 
morphology of the giant branch in the dense core of 47 Tuc differs markedly 
from that in the cluster outskirts. In particular, there appear to be fewer 
bright RGB stars in the core as well as an enhanced asymptotic giant branch 
(AGB) sequence. While a similar deficiency of bright giants has been observed 
in the cores of the massive GCs NGC 2808 and NGC 2419, better agreement 



between the observations and theoretical luminosi ty functions obtaine c 



the Victoria- Regina isoc h rones was found for M5 ([Sandquist fc Martel 



Sandquist fc Hess 



2008 



Sandquist fc Martell (120071 ) speculate that the giant 



with 



2007 



star observations in NGC 2808 could be linked to its unusually blue hori- 
zontal branch if a fraction of the stars near the tip of the RGB experience 
sufficiently enhanced mass loss that they leave the RGB early. Alternatively, 



Beer &; Daviesl (120041 ) suggest that RGB stars could be depleted in dense stel- 
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lar environments as a result of collisions between red giants and binaries. 

While stellar populations have been studied and compared on an indi- 
vidual cluster basis, a statistical analysis in which their core populations are 
compared over a large sample of GCs is ideal for isolating trends in their dif- 
ferences. Though a handful of studies of this nature have been performed (e.g. 



Piotto et al. 



20041 ). we present an alternative method by which quantitative 
constraints can be found for the relative sizes of different stellar populations. 
Specifically, a cluster-to-cluster variation in the central stellar mass function 
can be looked for by comparing the core masses to the sizes of their various 
stellar populations. Since stellar evolution is the principal factor affecting their 
relative numbers in the core, we expect the size of each stellar population to 
scale linearly with the core mass. If not, this could be evidence that other fac- 
tors, such as stellar dynamics, are playing an important role. In this chapter, 
we present a comparison of the core RGB, main-sequence (MS) and HB popu- 
lations of 56 GCs. In particular, we use star counts for each stellar population 
to show that RGB stars are either over-abundant in the least massive cores or 
under-abundant in the most massive cores, and that this effect is not seen for 
MS stars. We present the data in Section I3~2l and our methodology and results 
in Section 13.31 In Section 13.41 we discuss the implications of our results and 
explore various possibilities for the source of the observed discrepancy between 
RGB and MS stars. Concluding remarks are presented in Section 13.51 



3.2 The Data 



Colour-magnitude diagrams (CMDs) taken from 



Piotto et al 



(120021) are used 



to obtain star counts for the RGB, HB, MS and blue straggler (BS) populations 
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in the cores of 56 GCs. We apply the same selection criterion as outlined 
in Leigh, Sills & Knigge (2007) to derive our sample as well as to define 
the location of the main-sequence turn-off (MSTO) in the (F439W-F555W)- 
F555W plane. An example of this selection criterion, applied to the CMP of 



NGC 362, is shown in Figure 13.11 We in clude all stars in the 



Piotto et al. 



(120021 ) database. Since 



Piotto et al. 



( 120021 ) took their HST snapshots with the 



centre of the PC chip aligned with the cluster centre, a portion of the cluster 
core was not sampled for most GCs. We have therefore applied a geometric 
correction factor to the star counts in these clusters in order to obtain numbers 



2007 



2008|). 



that are representative of the entire core (ILeigh. Sills fc Knigge 
The total number of stars in the core is found by summing over all stars 
brighter than 1 mag below the MSTO and then multiplying by the appropriate 
geometric correction factor. 

Errors on the number of stars for each stellar population were calculated 
using Poisson statistics. Core radii, distance moduli, extinction corrections, 
central luminosity densities and central surface b rightne s ses were taken from 



the Harris Milky Way Globular Cluster catalogue ([Harris 



1996. 2010 updateh . 



Calibrated appa rent magnitudes in the F555W, F439W and Johnson V bands 



were taken from 



Piotto et al 



(120021 ) . 



3.3 Results 

This chapter focuses on the core RGB, MS and HB populations of 56 GCs, 
comparing their numbers to the core masses. Note that we are focusing on 
the total number of stars in the core as a proxy for the core mass instead of 
the total luminosity in the core in order to avoid concerns regarding cluster- 
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Figure 3.1 Colour-magnitude diagram for NGC 362 in the (F439W-F555W)- 
F555W plane. Boundaries enclosing the selected RGB, HB and MSTO popu- 
lations are shown. 



to-cluster variations in the central stellar mass function and selection effects. 
Given that a single bright HB star can be as luminous as 100 regular MS 
stars, a small surplus of bright stars could have a dramatic impact on the 
total luminosity. Therefore, the total number of stars in the core is a more 
direct and reliable estimate for the core mass than is the core luminosity. 

Upon plotting the logarithm of the number of core RGB stars versus the 
logarithm of the total number of stars in the core and performing a weighted 
least-squares fit, we find a relation of the form: 

hg(N RGB ) = (0.89 ± 0.03) log(iV core /10 3 ) + (2.04 ± 0.02) (3.1) 
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The sub-linear slope is either indicative of a surplus of RGB stars in the least 
massive cluster cores or a deficiency in the most massive cores. Errors for lines 
of best fit were found using a bootstrap methodology in which we generated 
1,000 fake data sets by randomly sampling (with replacement) RGB counts 
from the observations. We obtained lines of best fit for each fake data set, fit a 
Gaussian to the subsequent distribution and extracted its standard deviation. 
In order to avoid the additional uncertainty introduced into our RGB number 
counts from trying to distinguish AGB stars from RGB stars, as well as the 
difficulty in creating a selection criterion that is consistent from cluster-to- 
cluster when including the brightest portion of the RGB, stars that satisfy 
the RGB selection criterion shown in Figure 13.11 are referred to as RGB stars 
throughout this chapter. Note that it is the brightest portion of the RGB that 
should be the most affected by stellar evolution effects such as mass-loss. If we 
extend our selection criterion to include the entire RGB, however, our results 
remain unchanged. 

Interestingly, MS plus sub-giant branch stars (hereafter collectively re- 
ferred to as MSTO stars, the selection criterion for which is shown in Figure 
1) show a more linear relationship than do RGB stars and appear to dominate 
the central star counts. If we count only those stars having a F555W mag 
within half a magnitude above and below the turn-off, we obtain a relation of 
the form: 

log(AWo) = (1-02 ± 0.01) log(iV core /10 3 ) + (2.66 ± 0.01) (3.2) 

A nearly identical fit is found when counting only those stars having a F555W 
mag between the turn-off and one magnitude fainter than the turn-off. 
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We also tried plotting the logarithm of the number of core helium- 
burning stars (labeled HB in Figure 13. ip versus the logarithm of the number 
of stars in the core, yielding a relation of the form: 

hg(N HB ) = (0.91 ± 0.10) log(iV core /10 3 ) + (1.58 ± 0.05) (3.3) 



Note the large uncertainty associated with the fit, indicating that the slope is 
consistent with both those of the RGB and MSTO samples. We will discuss 
this stellar population further in Section 13.41 

The number of MSTO, RGB and HB stars are shown as a function of the 
total number of stars in the core in Figure I3~2l Interestingly, the blue stragglers 
in our sample also scale sub-linearly with core m ass, albeit more dramatically , 



2009|). 



obeying a relation of the form N BS ~ M"^.g (IKnigge. Leigh fc Sills! . 
Note that N core can be used interchangeably with M core . In this case, we obtain 
a fit of: 



log(N BS ) = (0.47 ± 0.06) log(iV core /10 3 ) + (1-22 ± 0.02) (3.4) 

In an effort to explore the influence of selection effects, we re-did our 
plots having removed from our sample clusters denser than log p > 10 5 L Q pc~ 3 
since we are the most likely to be under-counting stars in the most crowded 
cluster cores where blending of the stellar light is the most severe. This cut 
also removes from our sample the post-core collapse (PCC) clusters for which 
the core radii are poorly defined since King models are known to provide a 
poor fit to the observed surface brightness profiles in these clusters. Similarly, 
we applied a cut in the central surface brightness, removing from our sample 



103 



Ph.D. Thesis — N. Leigh 



McMaster - Physics & Astronomy 




o - 



-a -l o l 
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Figure 3.2 The number of RGB (red circles), MSTO (black squares) and HB 
(blue triangles) stars found within the cluster core plotted versus the total 
number of stars in the core brighter than 1 mag below the MSTO, along with 
the corresponding lines of best fit for the RGB and MSTO samples. 

clusters satisfying So < 15.1 V mag arcsec -2 . Finally, since clusters having 
both high surface brightnesses and small cores are the most likely to suffer 
from selection effects, we also tried adding to the aforementioned cut in S a 
cut in the angular core radius, removing clusters with r c < 0.05'. In all cases, 
the sub-linear power-law index reported for the RGB remains unchanged to 
within one standard deviation of our original result. Selection effects do not 
appear to be the source of the observed sub-linearity, though it is clear that 
its effects must properly be accounted for in future studies. 

In order to assess the effects of age-related cluster-to-cluster variations 
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in the stellar mass function, as well as to test our assumption that the number 
of stars in the core provides a reliable estimate for the core mass, we have 
obtained MSTO mass es for most of the GCs in our sample. We fit theoretical 



isochrones provided in 



Pols et al. 



( 119981 ) to the cluster CMDs, using the bluest 



point along the MS of a given isochrone as a pro xy for the MSTO m ass. 



Isochrones were calculated us ing the metallici t ies o 



cluster ages were taken from 



De Angeli et al 



Piotto et al. 



(12002! ) and 



(120051 ) using the Zinn & West 



(1984) metallicity scale. Core masses were estimated by multiplying the mass 
corresponding to the MSTO (to-msto) by the number of stars in the core 
brighter than 1 mag below the turn-off. This is a r easonable assumption g iven 



2005J ) and 



the very small dispersion in the ages of MW GCs (iDe Angeli et all 
the fact that we are only considering stars brighter than 1 mag below the TO. 
Consequently, the range of stellar masses upon which we are basing our number 
counts is very small. Our results remain entirely unchanged upon using M core 
~ N CO re m MSTO as a proxy for the core mass instead of pure number counts. 

In order to further check the sensitivity of our results to our estimate 
for the core masses, we re-did all plots shown in Figure 13.21 using various 
approximations for the total core luminosity instead of pure number counts. 
Core luminosities are calculated in the Johnson V band directly from the 
stellar fluxes which are summed over all stars in the core and then multi- 
plied by the appropriate geometric correction factor. We also adopted L core 



-irr^ Pn, where Pn is the central luminosity density in L Q pc 3 taken from 



Harris! ( 11996. 2010 update! ) . Additionally, since the number of core RGB stars 
is in reality a projected quantity, we tried plotting Nrgb versus L core = 7rr^S , 
where So is the central surface brightness in L Q pc~ 2 , so that we are con- 
sistently comparing two projected quantities. Finally, we can adopt slightly 
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more realistic estimates for the total core luminosity by integrating over King 
de nsity profiles. We fit si n gle-m ass King models calculated using the method 



of 



Sigurdsson fc Phinneyl (119951 ) to the surface brightness profiles of the ma- 



jority of the clusters in our 



s ample using the concentration parameters of 



:is 



McLaughlin &: van den Maxell (120051 ) and the central luminosity densities of 



Harris! (11996. 2010 update! ). We then integrated the derived luminosity den- 



sity profiles numerically in order to estimate the total stellar light co ntained 



Harris 



within the core. Aft er removing clusters labelled as post-core collapse in 
( 1996. 2010 update! for which King models are known to provide a poor fit, we 



once again compared the integrated core luminosities to the number of RGB 
stars in the core. For all four of these estimates for the total core luminosity, 
we find that our fundamental results remain unchanged, with the power-law 
index for RGB stars remaining sub-linear at the 3-cx level. Therefore, we con- 
clude that our result is robust to changes in choices of cluster and population 
parameters. 



3.4 Discussion 



We have shown that the number of RGB stars in globular cluster cores does 
not directly trace the total stellar population in those cores. In particular, the 
number of RGB (but not MSTO) stars scales sub-linearly with core mass at the 
3-cr level. Given that the MS lif etime is exp ected to be a factor of 10-100 longer 



than that of the RGB sample (jlben 



199l|) , the ratio Nmsto/^rgb indicates 



that the relative sizes of these stellar populations are in better agreement 
with the expectations of stellar evolution theory in the most massive cores. 
This suggests that our results are consistent with a surplus of RGB stars in 
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the least massive cores. We discuss below some of the key considerations in 
understanding the evolution of GC cores and the stars that populate them in 
an effort to explain our result. 



3.4.1 Stellar evolution 



Could this trend be a reflection of a stellar evolution process? The evolution 
and distribution of stellar populations can be thought of as the sum of many 
single stellar evolution tracks, which depend only on a star's mass and compo- 



(IHarris 



1996. 2010 UDdate 


) and 


small ( 


De Angeli et al.. 



20051 ). there is no reason to expect the RGB 
lifetime to depend on the cluster mass. On the other hand, recent studies 
suggest that the chemical self-enrichment of GCs during their early evolution- 
ary stages could help to explain many of t he population differences observed 



among them (e.g. 



Caloi fc D'Antona 



20071 ). In particular, many of the most 
ium and this i s exp ected to 



Romano et al. 



20071). While 



massive GCs are thought to be enriched in he. 
reduce the time scale for stellar evolution (e.g. 
this scenario predicts a deficiency of RGB stars in the most massive cores, it 
would contribute to depressing the slope of the RGB sample relative to that 
of the MSTO population in Figure 13.21 



3.4.2 Single star dynamics 

Two-body relaxation is the principal driving force behind the dynamical evo- 
lution of present-day GCs, slowly steering them towards a state of increased 
mass stratification as predominantly massive stars fall into the core and typi- 
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cally low-mass stars are ejected via dynamical encounters. The relaxation time 



19871 ) and the varia nce in the relative 



increases with the cluster mass (jSpitzer 
ages (and hence MSTO masses) of MW GCs is quite small (iDe Angeli et al. 



20051 ). Therefore, it is the least massive clusters that should show signs of 
being the most dynamically evolved. This assumes that cluster-to-cluster vari- 
ations in the initial mass function and the degree of initial mass segregation 
are small. In general, however, proportionately fewer massive stars should 
have had sufficient time to migrate into the most massive cores, while fewer 
low-mass stars should have been ejected out. While qualitatively correct, this 
effect should contribute little to the observed difference between the core RGB 
and MSTO populations sin ce RGB stars are only slightly more massive (e.g. 



De Marchi fc Pulone 



2003). 



Stars expand considerably as they ascend the RGB. Both the increase 
in collisional cross-section and the change in the average stellar density could 
have an important b earing o n the o utcomes of dynamical interactions involving 



RGB stars. Indeed, iBailynl ( 1l994l ) suggests that interactions between giants 
and other cluster members in the core could strip the outer envelope of the 
giant before it has a chance to fully ascend the RGB. Since our adopted RGB 
selection criterion does not include the brightest giants, we ar e only cons idering 



19911 ). This 



giants that are larger than MSTO stars by a factor of ~ 10 ( llben . 
small degree of expansion will have only a minor effect on the collision rate. 
Any scenario that relies on dynamical encounters to explain a depletion of 
RGB stars should be operating in very dense cores. Our results are consistent 
with some of the densest clusters in our sample having a surplus of giants, 
however. 
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3.4.3 Binary effects 



Stripping of the envelopes of large stars could also be mediat ed by a 



Dinary 



19941 ). 



companion as the expanding giant overfills its Roche lobe (IBailyn 
While this process should preferentially occur in the centres of clusters where 
binaries will congregate as a result of mass segregation, two-body relaxation 
progresses more slowly in the most massive clusters. Binaries should therefore 
sink to the cluster core more quickly in the least massive clusters, contributing 
to an increase in the core binary fraction at a rate that decreases with increas- 



ing cluster mass. Observationa . 



most notably by 



Sollima et al- 



ias been found in s upport of this 



Milone et al. 



(120081 ) who found an 



evid ence 
(120071 ) and 

anti-correlation between the cluster mass and the core binary fraction. Any 
mechanism for RGB depletion that relies on binary stars should therefore op- 
erate more efficiently in the least massive cores where the binary fraction is 
expected to be the highest. Our results are consistent with a surplus of giants 
in the least massive cores, however. This therefore argues against a binary 
mass-transfer origin for RGB depletion in massive GC cores. For similar rea- 
sons, it seems unlikely that our result can be explained by collisions between 
RGB stars and binaries. If, on the other hand, RGB stars are more com- 
monly found in binaries than are MS stars, perhaps as a result of their larger 
cross-sections for tidal capture, binary stars could still be contributing to the 
observed trend. Note that in the cluster outskirts where the velocity disper- 
sion has dropped considerably from its central value, individual encounters are 
more likely to result in tidal capture. Since mass segregation should deliver 
binaries to the core faster in the least massive clusters, a larger fraction of their 



RGB stars could have hitched a ride to the core as a binary companion. How- 
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ever, both the average half-mass relaxation time o f MW GCs and the RGB 



lifetim e tend to be on the order of a Gyr ( {Harris 



1996. 2010 update 



lb en 



199ll ). This does not leave much time for giants to be captured into binaries 



and subsequently fall into the core before evolving away from the RGB. 
3.4.4 Core helium-burning stars 

The fit for the HB sample is consistent with those of the RGB and MSTO 
samples at the 3-sigma level so that we are unable to draw any reliable conclu- 
sions for this stellar population. The high uncertainty stems from a number 
of outlying clusters. Selection effects and contamination from the Galaxy are 
likely to be playing a role in this, in addition to our formulaic selection crite- 
rion which may not be as suitable to the varying morphology of the HB as it 
is to other stellar populations. That is, the creation of a purely photometry- 
based cluster-independent selection criterion may not be possible for HB stars. 
Given that stellar evolution effects are expected to be the most dramatic at 
the end of the RGB lifetime, an interplay with the cluster dynamics could also 
be contributing. In particular, if the central relaxation time is shorter than 
the HB lifetime, significant numbers of HB stars could be ejected from the 
core via dynamical encounters as a result of having lost around a quarter of 
their mass upon evolving off the tip of the RGB. Moreover, since stars expand 
considerably as they ascend the RGB, many of the dynamical arguments pre- 
sented in Section 13.4.21 may more strongly affect the size of the HB sample if 
they are the direct evolutionary descendants of RGB stars. Since at most a 
handful of studies have b een performed compa ring the radial HB and RGB 



distributions in GCs (e.g. 



Iannicola et al. 
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any firm constraints can be placed on the source of the poor fit found for the 
HB sample. 

3.4.5 An evolutionary link with blue stragglers? 

The addition of a small number of extra RGB stars to every cluster is one 
way to account for the observed sub-linear dependence on core mass since 
the fractional increase in the size of the RGB population will be substantially 
larger in the least massive cores. In log-log space, the result is a reduction 
of the slope. Since blue straggle r s will evolve into RGB and eventually HB 



stars ( iSills. Karakas &: Lattanzio 



20091 ). evolved BSs could be the cause of 



a surplus of RGB (and possibly core helium-burning) stars in these clusters. 
This scenario also predicts that MSTO stars should scale slightly more linearly 
with core mass since there should be a smaller contribution from evolved BSs, 
as we have shown. Given the fits for the RGB and BS samples presented in 
Section [3^31 and their corresponding uncertainties, we find that the addition of 
evolved BSs to the RGB populations could inflate the slope enough that the 
dependence on core mass becomes linear. Upon subtracting the BS sample 
from the RGB sample, we find that the new fit is consistent with being linear: 

log(N RGB ) = (0.94 ± 0.04) log(iV core /10 3 ) + (1.97 ± 0.02) (3.5) 

The slope becomes larger if we have under-estimated the number of BSs, per- 
haps as a result of selection effects, our adopted selection criterion or a larger 
population size having existed in the past. 
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3.5 Summary 

In this chapter, we have performed a cluster-to-cluster comparison between 
the number of core RGB, MSTO & HB stars and the total core mass. We 
have introduced a technique for comparing stellar populations in clusters that 
is well suited to studies of both cluster and stellar evolution, in addition to 
the interplay thereof. Using a sample of 56 GCs taken from Piotto et al.'s 
2002 HST database, we find a sub-linear scaling for RGB stars at the 3-er 
level, whereas the relation is linear for MSTO stars. While the preferential 
self-enrichment of massive GCs, two-body relaxation, and evolved BSs could 
all be contributing to the observed sub-linear dependence, further studies with 
an emphasis on selection effects are needed in order to better constrain the 
source of this curious observational result. 
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Chapter 




Dissecting the Colour-Magnitude 
Diagram: A Homogeneous 
Catalogue of Stellar Populations in 
Globular Clusters 



4.1 Introduction 



Colour-magnitude diagrams (CMDs) are one of the most important tools avail- 
able to astronomers for studying stellar evolution, stellar populations and star 
clusters. And yet, there remain several features found in CMDs whose ori- 
gins are still a mystery. Examples include horizontal branch (HB) morphol- 



ogy, the presence 



glers (BSs) 



Potter et al. 



of exten d ed horizontal branch ( 



Sandage 



1953 



Zinn fc Barnes 



HB) stars, and blue strag- 



1996 



Peterson et al. 



2003 



20101 ) . Previous studies have shown that the observed differ- 



ences in the HBs of Milky Way globular clusters (GCs) are related to metal- 
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licity ( jSandage fc Wallersteinl . Il960l ). however at least one additional param- 
eter is required to explain the spread in their colours. Many cluster proper- 
ties have been suggested as possible Second and Third Parameters, including 



age, central density and c 
have been identified (e.g. 



uster l uminosity, although no defini tive candidates 



Rood . 



1973 



Fusi Pecci et al. 



19931 ) . An explana- 



tion to account for the existence of BSs has proved equally elusive. Many BS 
formation mechanisms have been p roposed, including stellar collision s (e.g . 



Leonard. 



1989 



Si 



Mathieu fc Geller 



Is fc Bailvn 



19991 ) and binary mass-transfer (IMcCrea 



1964 



20091 ). However, no clear evidence has yet emerged in favour 



of a dominant formation mechanism. 

In short, we still do not understand how many of the physical processes 
operating within star clusters should affect the appearance of CMDs (e.g. 



Fusi 



Pecci et al 



1992 



Buonanno et al. 



1997 



Ferraro et al. 



1999 



Beccari et al. 



20061 ). In general, the importance of these processes can be constrained by 
looking for correlations between particular features in CMDs and cluster prop- 
erties that serve as proxies for different effects. For example, the central density 
can be used as a rough proxy for the frequency with which close dynamical en- 
counters occur. Similarly, the cluster mass can be used as a proxy for the rate 
of two-body relaxation. Once the relevant effects are accounted for, CMDs 
can continue to provide an ideal tool to further our understanding of stellar 
evolution, stellar populations and star clusters. 

It is now clear that an important interplay occurs in clusters between 
stellar dynamics and stellar evolution. For example, dynamical models have 
shown that star clusters expand in response to mass-loss driven by stellar 
evolution, particularly du ring their early evolutionary phases when massive 



stars are still present (e.g. 



Chernoff &; Weinberg 



1990 



Portegies Zwart et al. 
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19981 ; iGieles et al.l . |2010| ). Mass- loss resulting from stellar evolution has also 
been proposed to cause horizontal branch stars to exhibit more extended radial 
distributions relative to red giant branch and main-sequence turn-off (MSTO) 



stars in globular clusters 
avera ge HB lifetime (e.g. 



laving short central re l axation times relative to the 



Sigurdsson fc Phinney 



1995 



Leigh. Sills &: Knigge 



20091 ). This can be understood as follows. Red giant branch (RGB) stars 
should be more mass segregated than other stellar populations since they are 
among the most massive stars in GCs. HB stars, on the other hand, are much 
less massive since RGB stars experience significant mass loss upon evolving into 
HB stars. Consequently, two-body relaxation should act to re-distribute HB 
stars to wider orbits within the cluster potential relative to RGB and MSTO 



stars (ISpitzer fc Shull 



19751 ). provided the average HB lifetime is shorter than 
the central relaxation time. Studies have shown that the radial distributions 
of the HB populations in some GCs cou l d diff er from those of other stellar 



populations. For instance, 



Saviane et al. 



( 119981 ) presented evidence that blue 



HB stars could be more centrally concentrated than red HB and s ub-giant 



Cohen et al 



fll997h showed 



branch stars in the GC NGC 1851. Conversely, 
that blue HB stars could be centrally depleted relative to other stellar types 
in the GC NGC 6205. To date, no clear evidence has been found linking the 
spatial distributions of HB stars to any global cluster properties. 

Peculiar trends have also been reported for the radial distributions of 
RGB stars. For example, a deficiency of bright r ed giants has been observed 



in the GC NGC 1851 (e.g. 



Iannicola et al 



20091 ) 



Sandquist fc Martell (120071 ) 



discussed the possibility that this deficiency could be the result of strong mass 
loss on the RGB. Alternatively, some authors have suggested that dynamical 
effects could deplete red giants. For instance, giants could experience collisions 
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more frequen tly than other stellar p opulations due to their larger cross- sect ions 



for collision f lBeers fc Davies . 



20041 ) 



One important example of the interplay that occurs in clusters be- 
tween stellar evolution and stellar dynamics can be found in the study of 
blue stragglers. Found commonly in both open and globular clusters, BSs 
are thought to be produced by the addition of hydrogen to the cores of low- 



mass main-sequence (MS^ 



MSTO in cluster CMDs flSandage . 



stars, and the refore appear as an extension of the 



19531 ) . This can occur via multiple chan- 



nels, most of which involve the mergers of low-mass MS stars since a significant 
amount of mass is typically required to reproduce the observed locations of 



BSs in CMDs (e.g. 



Sills fc Bailyn 



19991 ) . Stars in close binaries can merge if 



enough orbital angular momentum is lost, which can be mediated by dynam- 



or even an outer triple companion (e.g. 


Leonard & Linnell, 


1992: 


Li & Zhang. 


2006; 


Perets &; Fabrvckv. 


2009; 


Dervisoglu. Tout & Ibanoglu. 


2010 


). Alterna- 



tively, MS stars can collide directly, altho ugh this i s also 



mediated by 



Leigh fc Sills 



multip le star systems (e.g. 



Leonard . 



1989 



thought to usu a 



Fregeau et al- 



ly be 



2004 



20101 ) . Finally, BSs have also been hypothesized to form by 



mass-transfer from a n evolving pri mary onto a normal MS companion via 



Roche lobe overflow (McCrea 



196J). 



Whatever the dominant BS formation mechanism(s) operating in dense 
star clusters, it is now thought to somehow involve mult iple star systems. This 



was s hown to be the case in even the dense cores of GCs (IKnigge. Leigh Sz Sills 



200 9^ where co llisi ons between single stars are thought to occur frequently 



( Leonard 



19891). In 



Knigge. Leigh &; Sills! (120091 ). we showed that the numbers 



of BSs in the cores of a large sample of GCs correlate with the core masses. 
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We argued that our results ar e consistent with what i s expected if BSs are 
descended from binary stars. iMathieu Sz Gellerl (120091 ) also showed that at 



least 76% of the BSs in the old open cluster NGC 188 have binary companions. 
Although the nature of these companions remains unknown, it is clear that 
binaries played a role in the formation of these BSs. Dynamical interactions 
occur frequently enough in dense clusters that they are also expected to be at 
least partly responsible for the observed properties of BSs. It follows that the 
current properties of BS populations should reflect the dynamical histories of 
their host clusters. As a result, BSs could provide an indirect means of probing 
the physical processes that drive star cluster evolution. 

In this chapter, we present a homogeneous catalogue for red giant 
branch, main-sequence turn-off, horizontal branch and blue straggler stars in 
a sample of 35 Milky Way (MW) GCs taken f rom the Advanced Cam era for 



20070 . With 



Surveys (ACS) Survey for Globular Clusters (jSarajedini et al. . 
this catalogue, we investigate two important issues related to stellar popu- 
lations in G Cs. First, we test th e obse rvational correlation found for BSs 



presented in 



Knigge. Leigh fc Sillsl (]2009[ ) by re-doing the study with newer 



and more accurate photometry. The larger spatial coverage considered in our 
new sample offers an important additional constraint for the origin of this cor- 
relation. Second, we perform the same statistical comparison for RGB, HB 
and MSTO stars in order to study their radial distributi ons. This will allow us 



to test some of the results and hypotheses introduced in 



Leigh. Sills fc Knigge 



(120091 ). where we first presented this technique for studying stellar populations. 
In particular, we found evidence for a surplus of RGB stars in low-mass GC 
cores relative to MSTO stars. However, we concluded that the study needed 
to be re-done with better photometry. The ACS data are of sufficiently high 
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quality to address this issue. 

In Section 14.21 we present our selection criteria to determine the num- 
bers of BS, RGB, HB and MSTO stars located in the central cluster regions. 
The spatial coverage of the photometry extends out to several core radii from 
the cluster centre for most of the clusters in our sample. For these clusters, 
we have obtained number counts within several circles centred on the cluster 



centres provided in 



Goldsbury et al. 



(120 101 ) for various multiples of the core 



radius. This catalogue is presented in Section 14.31 In this section, we also 
present a comparison between the sizes of the different stellar populations and 
the total stellar masses contained within each circle and annulus. Finally, we 
discuss our results for both BSs and the other stellar populations in Section l4~4l 



4.2 Method 

In this section, we present our sample of CMDs and define our selection cri- 
teria for each of the different stellar populations. We also discuss the spatial 
coverage offered by the ACS sample, and describe how we obtain estimates for 
several different fractions of the total cluster mass from King models. 



4.2.1 The Data 



The data used in this study consist s of a sample of 35 MW GCs taken from the 

~0 



ACS Survey for Globular Clusters (jSarajedini et al. 



20071 ) J 1 ! The ACS Survey 



provides unprecedented deep photometry in the F606W (~ V) and F814W (~ 
I) filters that is nearly complete down to ~ 0.2 M . In other words, the CMDs 
extend reliably from the HB all the way down to about 7 magnitudes below 



lr The data can be found at http://www.astro.ufl.edu/ ata/publicjistgc/. 
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the MSTO. We have confirmed that the photometry is nearly complete above 
at least 0.5 magnitudes below the MSTO for every cluster in o ur sample. This 



was d one using the results of artificial star tests taken from 



Anderson et al. 



(120081 ). and confirms that the photometric quality of the stellar population 



catalogue presented in this chapter is very higho We have also considered 
foreground contamination by field stars, and it is negligible. 

Each cluster was centred in the ACS field, which extends out to sev- 
eral core radii from the cluster centre in most clusters and, in a few cases, 
beyond even 15 core ra dii. Coordinates for the cluster centres were taken from 



Goldsbury et al. 



(120101 ) . These authors found their centres by fitting a series 



of ellipses to the density distributions within the inner 2' of the clust er cen 



tre, and computing an average value. The core radii were taken from 



Harris 



(11996. 2010 updatel ). 



4.2.2 Stellar Population Selection Criteria 

In order to select the number of stars belonging to each stellar population, 
we define a series of lines in the (F606W-F814W)-F814W plane that act as 



boundaries enclosing each of the di fferent stellar 



fit theoretical isochrones taken from 



Potter et al. 



oopu 



at ions. To do this, we 



(120071 ) to the CMDs of every 



cluster in our sample. Each isochrone was generated using the metallicity and 
age of the cluster, and fit to its CMP using the c orresponding distance modulus 



and extinction provided in 



Potter et al. 



(120 10l ). The MSTO was then defined 



using our isochrone fits by selecting the bluest point along the MS. This acts 
as our primary point of reference for defining the boundaries in the CMP for 



2 Artificial star tests were obtained directly from Ata Sarajedini via private communica- 



tion. 
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the different stellar populations. Consequently, the selection criteria provided 
in this chapter are a s ignifi cant improvement upon the criteria presented in 



Leigh. Sills fc Kniggd (120071 1 . and our new catalogue for the different stellar 
populations is highly homogeneous. 

Two additional points of reference must also be defined in order for 
our selection criteria to be applied consistently from cluster-to-cluster. First, 
the selection criteria for the HB are determined by fitting a line through the 
approximate mid-point of the points that populate it in the CMD. This line 
is then used to define upper and lower boundaries for the HB. Theoretical 
isochrones become highly uncertain at the HB, so it is necessary to specify 
this additional criterion by eye. Second, the lower boundary of the RGB is 
defined for each cluster as the point along its isochrone corresponding to a 
helium core mass of 0.08 M Q . We do not include RGB stars brighter than the 
HB since the tilt of the upper RGB varies significantly from cluster-to-cluster, 
presenting a considerable challenge for the consistency of our selection criteria. 
Moreover, the distinction between RGB and asymptotic giant branch stars in 
the CMD is often ambiguous. 

Example selection criteria for each of the different stellar populations 
are shown in Figure 14.11 Formal definitions for the boundaries in the cluster 
CMD that define the BS, RGB, HB and MSTO populations are provided in 
Appendix [B] We note that the sizes of our selection boxes have been chosen 
to accommodate the photometric errors, which contribute to broadening the 
various evolutionary sequences in the CMD. With superior photometry, the 
sizes of our selection boxes could therefore be reduced. This would further 
decrease contamination from field stars in our samples. 
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1 

F606-F614 



Figure 4.1 Colour- magnitude diagram for the Milky Way globular cluster NGC 
6205. Boundaries enclosing the parameter space in the (F606W-F814W)- 
F814W plane that define each of the different stellar populations are indicated 
with solid lines, as described in the text. Absolute magnitudes are shown, con- 
verted from apparent magnitud es using the distance moduli and extinctions 



provided in iDotter et al.l ( 120101 ) . Labels for blue straggler, red giant branch, 
horizontal branch and main-sequence turn-off stars are indicated. Stars with 
large photometric errors have been omitted from this plot. 



4.2.3 Spatial Coverage 



The ACS field of view extends out to several core radii from the cluster centre 
for nearly every cluster in our sample. Consequently, we have obtained esti- 
mates for the number of stars contained within four different circles centred 
on the central cluster coordinates provided in 
these numbers only for clusters for which the indicated circle is completely 
sampled by the field of view. The radii of the circles were taken to be integer 



Goldsburv et al.l (J2010|). We list 
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multiples of the core radius, and we focus our attention on the inner four core 
radii since the field of view extends beyond this for only a handful of the clus- 
ters in our sample. An example of this is shown in Figure 14.21 The numbers 
we list are cumulative, so that entries for each circle include stars contained 
within all smaller circles. 




2S0.3B 250.4 250/42 250.44 250.4G 

RA 



Figure 4.2 RA and Dec coordinates for all stars in the GC NGC 6205. Circles 
corresponding to one, two, three and four core radii are shown. 



4.2.4 King Models 



In order to obtain accurate estimates for the total stellar mass contained 
within eac h circle, we generated s i ngle- m ass King models calculated using the 



method of 



Sigurdsson fc Phinneyl (119951 ) to obtain luminosity density profiles 
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for the majority of the clusters in o ur sample. The profiles were obtaine d, us- 



ing the concentration parameters of 



th e central luminosity densities o f 



m 



McLaughlin fc van den Marell (120051 ) and 



Harris! ( 11996. 2010 update! ) for each cluster 



McLaughlin fc van den Marell (120051 ) that overlaps with our sample. We 



then integrated the derived luminosity density profiles numerically in order 
to estimate the total stellar light contained withi n each circle. After remov - 



1996. 2010 update! ) 



ing clusters with high concentration parameters (lHarris 
for which King models are known to provide a poor fit, we multiplied the 
total stellar light by a mass-to-light ratio of 2 in order to obtain estimates 
for the total stellar mass contained within each circle. Calculating the total 
stellar mass contained within each circle from King models requires a number 
of assumptions that we will discuss fully in Section 14.41 



4.3 Results 

In this section, we present our catalogue along with the results of our compar- 
isons between the sizes of the different stellar populations and the total stellar 
mass contained within each circle and annulus. 



4.3.1 Catalogue 

The numbers of BS, RGB, HB and MSTO stars found within several different 
circles are shown for all clusters in Table 14.11 along with the total number 
of stars with magnitudes brighter than 0.5 mag below the MSTO. Number 
counts are only shown whenever the spatial coverage is complete within the 
indicated circle. 
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4.3.2 Population Statistics 



How can we use our catalogue to learn which, if any, cluster properties affect 
the appearance of CMDs? One way to accomplish this is by plotting the size 
of a given stellar population in a particular circle versus the total stellar mass 
contained within it. From this, lines of best fit can be found that provide 
equations relating the size of each stellar population to the total stellar mass 
contained within each circle. As described below, this is ideal for probing the 
effects of the cluster dynamics on the appearance of CMDs. 



The rate of two-body re 



the half-mass relaxation time ( jSpitzer 



axation 



or a c luster can be approximated using 



1983): 



t rh = 1.7 x ltf[r h {pc)f 2 N l l 2 [m/M & y l l 2 



years, 



(4-1) 



where is the half-mass radius, N is the total number of stars within 
and m is the average stellar mass. The half-mass radii of MW GCs are re- 
markably similar independent of mass, and simulations have shown that 



changes by a 



( Murray 



2009 



'actor o 



Henon 



at most a few over the course of a cluster's lifetime 



19731 ). The GCs that comprise our sample show a 



range of mass es spanning roughly 3 orders of magnitude, and have compara- 



bly old ages (IDe Angeli et al. 



20051 ). Therefore, Equation 14.11 suggests that 
the degree of dynamical evolution (due to two-body relaxation) experienced by 
a cluster is primarily determined by the total cluster mass for the GCs in our 
sample. In particular, more massive clusters are less dynamically evolved, and 
vice versa. Consequently, if the size of a given stellar population is affected by 
two-body relaxation, the effects should be the most pronounced in the least 
massive clusters in our sample. Additionally, the rate of (direct) stellar colli- 
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sions increases with increasing cluster mass (e.g. 



Davies. Piotto fc De Angeli 



20041 ). This suggests that, if a given stellar population is affected by collisions, 
the effects should be the most pronounced in the most massive clusters in our 
sample. Therefore, by comparing the size of each stellar population to the total 
stellar mass contained within a given circle, the effects of the cluster dynamics 
can be quantified. This technique also ensures a normalized and consistent 
comparison since it accounts for cluster-to-cluster differences in the fractional 
area sampled by the ACS field of view. That is, we are consistently comparing 
the same structural area for each cluster. The validity and implications of all 
of these assumptions will be discussed further in Section 14.41 

Plots showing the number of stars belonging to each stellar population 
as a function of the total stellar mass contained within each circle are shown 
in Figure 14.31 Uncertainties for the number of stars belonging to each stellar 
population were calculated using Poisson statistics. We also plot in Figure 14.41 
the number of stars belonging to each stellar population as a function of the 
total stellar mass contained in each annulus outside the core. That is, we 
considered the populations for each annulus individually, as opposed to con- 
sidering every star with a distance from the cluster centre smaller than the 
radius of the outer-most circles. Recall that we have neglected clusters for 
which our theoretical King models provide a poor description of the true den- 
sity distributions. This was the case for clusters in our sample having a high 
concen t ration parameter, most of which are labelled as post-core collapse in 



Harrisl ( 11996. 2010 updatel ). 
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Figure 4.3 The logarithm of the number of stars belonging to each stellar 
population is shown for each circle as a function of the logarithm of the total 
stellar mass. From left to right and top to bottom, each frame corresponds 
to number counts contained within a circle having a radius of r c , 2r c , 3r c 
and 4r c . Blue corresponds to blue stragglers, red to red giant branch stars, 
green to horizontal branch stars and black to main-sequence turn-off stars. 
Estimates for the total stellar mass contained within each circle were found 
using single-mass King models, as described in Section I4.2.4I 
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Figure 4.4 The logarithm of the number of stars belonging to each stellar 
population is shown for each annulus as a function of the logarithm of the 
total stellar mass. The annulus and colour corresponding to each inset and 
stellar population, respectively, are the same as in Figure 14.31 



131 



Table 4.2 Lines of Best Fit for log(M^ rde /10 3 ) Versus log(N pop ) 



Circle 


BS 


RGB 


HB 


MSTO 


< r c 


log(N BS ) 


(0.39 


± 


log(N RGB ) = 


(0.95 ± 


log(N HB ) 


(0.95 


± 


log(N J/STO ) = 


(0.90 


± 




0.05)log(M <I . c /10 3 ) 


+ (1.22 ± 


0.11)log(M <I> /10 3 ) 


+ (1.36 ± 


0.06)log(M <rc /10 3 ) 


+ (0.92 ± 


0.07)log(M <rc /10 3 ) 


+ (2.03 ± 




0.05) 






0.11) 




0.06) 






0.07) 






< 2r c 


l°g(N BS ) 


(0.36 


± 


log(N RGS ) = 


(0.87 ± 


l°g(N HB ) = 


(0.85 


± 


log(N J/STO ) = 


(0.82 


± 




0.05)log(M <2rc /10 3 ) 


+ 


(1.26 


0.08)log(M <2rc /10 3 ) 


+ (1.27 


0.06)log(M <2rc /10 3 ) 


+ 


(0.84 


0.06)log(M <2re /10 3 ) 


+ 


(1.98 




± 0.08) 






± 0.12) 




± 0.09) 






± 0.09) 






< 3r c 


log(N BS ) 


(0.47 


± 


log(N HGB ) = 


(0.80 ± 


log(N HB ) 


(0.79 


± 


log(N A j STO ) = 


(0.82 


± 




0.04)log(M <3 ,. c /10 3 ) 


+ 


(1.02 


0.06)log(M <3rc /10 3 ) 


+ (1.26 


0.08)log(M <3rc /10 3 ) 


+ 


(0.83 


0.11)log(M <3rc /10 3 ) 


+ 


(1.86 




± 0.08) 






± 0.12) 




± 0.14) 






± 0.20) 






< 4r c 


log(N BS ) 


(0.45 


± 


log(N HGS ) = 


(0.75 ± 


log(N Hi3 ) = 


(0.75 


± 


log(NA/STO) = 


(0.78 


± 




0.05)log(M <4l . c /10 3 ) 


+ 


(1.01 


0.07)log(M <4rc /10 3 ) 


+ (1.28 


0.09)log(M <4rc /10 3 ) 


+ 


(0.83 


0.12)log(M <4rc /10 3 ) 


+ 


(1.86 




± 0.12) 






± 0.15) 




± 0.19) 






± 0.25) 







Table 4.3 Lines of Best Fit for log(M ann?t ; MS /10 3 ) Versus log(N pop ) 



Annulus 


BS 


RGB 


HB 


MSTO 


r c < r < 2r c 
2r c < r < 3r c 
3r c < r < 4r c 


log(N BS ) (0.27 ± 
0.08)log(M,. c<r<2rc /10 3 ) + (1.04 
± 0.13) 

log(N BS ) = (0.39 ± 
0.09)log(M 2l . c<r<3l . c /10 3 ) + 
(0.52 ± 0.16) 

log(N BS ) = (0.15 ± 
0.21)log(M 3l . c<r< 4 I . c /10 3 ) + 
(0.77 ± 0.35) 


log(N R Gs) = (0.80 ± 
0.06)log(M rc<r<2rc /10 3 ) + (1.22 
± 0.08) 

log(N HGB ) = (0.79 ± 
0.11)log(M 2l . c<r<3l . e /10 3 ) + 
(0.97 ± 0.17) 

log(N RGB ) = (0.63 ± 
0.11)log(M 3l . c<r< 4 I . c /10 3 ) + 
(1.04 ± 0.18) 


log(N HB ) (0.77 ± 
0.06)log(M rc<r<2rc /10 3 ) + (0.79 
± 0.09) 

log(N BB ) = (0.68 ± 
0.10)log(M 2rc<r . <3rc /10 3 ) + 
(0.67 ± 0.16) 

l°g(Njf B ) = (0.69 ± 
0.18)log(M 3rc<r<4l . c /10 3 ) + 
(0.45 ± 0.31) 


log(N MSTO ) = (0.76 ± 
0.05)log(M rc<r . <2rc /10 3 ) + (1.91 
± 0.08) 

log(N A/STO ) = (0.78 ± 
0.12)log(M 2rc<r<3rc /10 3 ) + 
(1.64 ± 0.20) 

log(N A/STO ) = (0.69 ± 
0.16)log(M 3rc<r<4rc /10 3 ) + 
(1.60 ± 0.27) 
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We performed a weighted least-squares fit for every relation in Fig- 
ure 14.31 and Figure 14.41 Slopes and y-intercepts for these lines are shown in 
Table 14.21 and Table 14.31 respectively. Uncertainties for the slopes and y- 
intercepts were found using a bootstrap methodology in which we generated 
1,000 fake data sets by randomly sampling (with replacement) number counts 
from the observations. We obtained lines of best fit for each fake data set, fit a 
Gaussian to the subsequent distribution and extracted its standard deviation. 

As shown in Table I4T21 the power-law index is sub-linear for BSs within 
the core at much better than the 3— a confidence level, and it is consistent with 
the sl ope obtained in our earlier analysis presented in 



Knigge. Leigh &: Sills 



( 20091 ). The slopes for the BSs are also sub-linear at better than the 3 



a confidence level for all circles outside the core. This is also the case for 
all annuli outside the core, as shown in Table 14.31 Note, however, that the 
uncertainties for the BS slopes are very large for all annuli outside the core, 
whereas this is not always the case for corresponding circles outside the core. 
This is the result of the fact that the number of BSs drops off rapidly outside 
the core in several clusters so that the corresponding Poisson uncertainties, 
which are given by the square-root of the number of BSs, are significant. The 
rapid decline of BS numbers with increasing distance from the cluster centre 
in these clusters has also contributed to an increased degree of scatter in the 
relations for annuli outside the core relative to the corresponding relations for 
circles outside the core. 

The slopes are consistent with being linear for all other stellar pop- 
ulations in the core within their respective 3 — a confi dence intervals. This 



agrees with the results of our earlier analysis presented in 



Leigh. Sills fc Knigge 



( 120091 ) when we performed the comparison using the total core masses. The 
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slopes are also consistent with being linear for all circles outside the core for 
both HB and MSTO stars. The power-law indices are sub-linear at the 3 — o 
confidence level only for RGB stars, and this is only the case for circles outside 
the core. The power-law index is nearly unity for the core RGB population, 
yet the associated uncertainty is very large. Upon closer inspection, the dis- 
tribution of power-law indices obtained from our bootstrap analysis for RGB 
stars in the core is strongly bi-modal, with comparably-sized peaks centred at 
~ 0.82 and ~ 1.0. This bi-modality is most likely an artifact of our bootstrap 
analysis caused by a chance alignment of data points in the log M core -log Nrgb 
plane. Upon performing the comparison for only those stars found within par- 
ticular annuli, our results suggest that the slopes are consistent with being 
sub-linear for all stellar populations at the 3 — a confidence level in only the 
annulus immediately outside the core (i.e. r c < r < 2r c ). The slopes are 
consistent with being linear for RGB, HB and MSTO stars in all other annuli. 

We also tried performing the same comparisons using the total number 
of stars in each circle and annulus as a proxy for the total stellar mass. In this 
case, the slopes are sub-linear for BSs within all circles and annuli at the 3 — a 
confidence level. Once again, the uncertainties are very large for all annuli 
outside the core, whereas this is not the case for corresponding circles outside 
the core. The slopes are consistent with being linear within the 1— a confidence 
interval for all other stellar populations in all ci rcles and annuli. Our results 



Leigh. Sills fc Kniggel (120091 ) 



are therefore inconsistent with those presented in 
for the core RGB populations, in which we found that RGB numbers scale sub- 
linearly with the number of stars in the core at the 3 — a confidence level. We 



will discuss the implications of these new results in Section 14.4 
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4.3.3 Blue Stragglers and Single-Single Collisions 



Knigge. Leigh fc Sillsl (120091 ). we 



As a check of our previous results reported in 
also looked for a correlation between the observed number of BSs in the cluster 
core and the number predicted from single-single (1+1) collisions. The results 
of this comparison are shown in Figure H~5l We define the predicted number of 
BSs formed from 1+1 collisions as AT 1+1 = tbs /t~i+i, where tbs is the average 
BS lifetime and t±+i is the average time between 1+1 c ollisions in the cluster 



core. We adopt the same definition for t\ + \ as used in 



Knigge. Leigh fc Sills 



(120091 ). and assume tbs = 1-5 Gyrs as well as an average stellar mass and 
radius of 0.5 M Q and 0.5 R , respectively. We also adopt a constant mass- 
to-light ratio of M/L = 2 for all cl usters. Central luminosity den sities and 



veloci ty dispersions were taken from 



Harrisl (11996. 2010 update! ) and 



Webbink 



(119851 ). respectively. 

Upon performing a weighted line of best fit 



pie that overlaps with the catalogue of 



'or ev ery cluster in our sam- 



Webbinkl (119851 ). we find a power-law 



index of 0.15 ± 0.03 (the uncertainty was found using the bootstrap method- 
ology described in Section 14.3.21) . For the subset of dense clusters having a 
central luminosity density satisfying log po > 4, we find a power-law index of 
0.36 ± 0.14. As before, we find no significant correlation with collision rate, 
even for the subset of dense clusters. Although we do find a weak dependence 
of BS numbers on collision rate for the entire sample, this is not unexpected 
since the collision rate and the core mass are themselves correlated, and our 
results suggest that there exists a strong correlation between BS numbers and 
the core masses. 
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Figure 4.5 The logarithm of the number of BSs predicted to have formed in 
the core from single-single collisions N 1+1 versus the logarithm of the observed 
number of BSs in the core Ng^. Filled circles correspond to clusters having 
central luminosity densities satisfying log po < 4, whereas open triangles cor- 
respond to dense clusters for which log p > 4. The adopted definition for 
Ni + i has been provided in the text. 

4.4 Summary & Discussion 

In this chapter, we have presented a catalogue for BS, RGB, HB and MSTO 
stars in a sample of 35 GCs. Our catalogue provides number counts for each 
stellar population within several different circles centred on the cluster cen- 
tre. The radii of the circles were taken to be integer multiples of the core 
radius, and we have focussed on the inner four core radii since the field of 
view extends beyond this for only a handful of the clusters in our sample. 
Particular consideration was given to our selection criteria for the different 
stellar populations in order to ensure that they were applied consistently from 
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cluster-t o-cluster. In particular , we h ave improved upon our previous selection 



criteria ( jLeigh. Sills fc Knigge 



20071 ) by fitting theoretical isochrones to the 
cluster CMDs. This provides an unambiguous definition for the location of 
the MSTO, which acts as the primary point of reference for the application of 
our selection criteria. As a result, our new catalogue is highly homogeneous. 

We have used our catalogue to quantify the dependence of the size 
of each stellar population on the total stellar mass enclosed within the same 
radius. As described in Section 14.3.21 this provides a means of quantifying 
the effects, if any, had by the cluster dynamics in shaping the appearance of 
CMDs above the MSTO. Below, we summarize the implications of our results 
for each of the different stellar populations. 



4.4.1 Blue Stragglers 



We have confirmed our previous result that the numbers of BSs in the core s 



2009|). 



of GCs scale sub-linearly with the core masses (IKnigge. Leigh fc Sillsl . 
That is, we find proportionately larger BS populations in low-mass GCs. There 
exist several possibilities that could explain the origin of this sub-linear de- 
pendence. First, we previously suggested that this could be an artifact of an 
anti-correlation between the binary fraction and the cluster (or core) mass 



(IKnigge. Leigh &: Sills 



20091 ) . This assertion stems from the fact that, if BSs 



have a binary origin, we expect their numbers to scale with the core mass 
as Nbs ~ fbM core , where /& is the binary fraction in the core. As before, 
we find that N BS ~ M? o ^~ ' 5 . Our resul t could therefore be e xplain ed if 
Second, we suggested in 



, -(0.5-0 

Jb ~ Mcore 



Leigh. Sills fc Knigge! ((20091) that 



the fact that the least massive GCs in our sample should be more dynamically 
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evolved than their more massive counterparts could be contributing to the ob- 
served sub-linearity for BSs. In particular, the low-mass clusters in our sample 
should have experienced a significant depletion of their very low- mass stars as a 



result of stellar evaporation induced by two-body relaxation (e.g. 


Soitzer. 


1987; 


Heggie & Hut. 


2003; 


De Marchi. Paresce & Portegies Zwart. 


2010 


). In turn, 



this could contribute to a higher fraction of merger products having masses 
that exceed that of the MSTO in low-mass GCs. As a result, more merger 
products would appear brighter and bluer than the MSTO in these clusters' 
CMDs, leading to more merger products being identified as BSs. Finally, mass 
segregation could also be contributing to the observed sub-linear dependence 
for BSs. Again, this is the result of the fact that the rate of two-body relax- 
ation, and therefore dynamical friction, is in general the fastest in low-mass 
clusters. BSs are among the most massive stars in GCs, so they should rapidly 
migrate into the core via dynamical friction in clusters for which the half-mass 
relaxation time is shorter than the average BS lifetime. It follows that propor- 
tionately more BSs could have drifted into the core via dynamical friction in 
low-mass GCs. This would also contribute to the observed sub-linear depen- 
dence of BS numbers on the core masses. 

This last hypothesis can be tested by comparing our scaling relations 
for progressively larger circles outside the core. If mass segregation is indeed 
the cause of the observed sub-linear dependence of BS numbers on the core 
masses, then we might expect the power-law index to systematically increase 
as we consider progressively larger circles. That is, we could be including 
more BSs that have not yet migrated into the core via dynamical friction, 
particularly in the most massive clusters in our sample. However, our results 
suggest that the power-law index remains roughly constant for all circles. This 
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is the case for both the comparison to the total stellar masses as well as to the 
total number of stars contained in each circle. The fact that these relations 
are comparably sub-linear within all circles can be interpreted as evidence 
that mass segregation is not the dominant effect contributing to the observed 
sub-linear dependence of BS numbers on the total stellar masses (or the total 
number of stars). We note that for many of the clusters in our sample, the 
spatial coverage is comparable to or exceeds the half-mass radius. This is a 
sufficiently large fraction of the total cluster area for our comparison to be 
sensitive to the effects of mass segregation. 



On the other hand, several GCs are known to show evidence 
modal BS radial distribution (e.g. 



Mapelli et al. 



2006 



Lanzoni et al. 



or a bi- 



2003). 



That is, in these clusters the number of BSs is highest in the central cluster 
regions and decreases with increasing distance from the cluster centre until 
a second rise in BS numbers occurs in the cluster outskirts. This secondary 
outer peak has been shown to occur at a distance from the cluster centre that 
exceeds 20 core radii in several cases. Consequently, the spatial coverage pro- 
vided by the ACS data likely does not extend sufficiently far in most clusters 
to detect any bi-modality in the BS radial distribution. Nonetheless, if applied 
to a statistically-significant sample for which the spatial coverage is complete 
out to the tidal radius, the technique we have presented in this chapter could 
provide a powerful constraint for the origin of the bi-modal BS radial distri- 
bution observed in several MW GCs by addressing the role played by mass 
segregation. 

We also tried to correlate the number of BSs observed in the cluster core 
with t he number predicted from single-single collisions. As in[ 



Knigge. Leigh &: Sills 



( 20091 ). we find that BS numbers depend strongly on core mass, but not on 

140 



Ph.D. Thesis — N. Leigh 



McMaster - Physics & Astronomy 



collision rate. This also proved to be the case for the subset of dense clusters 
satisfying log po > 4. Our previous interpretation that our results provide 
strong evidence for a binary, as opposed to a collisional, origin for BSs in GCs 
therefore remains the same. 



4.4.2 Red Giant Branch Stars 



The technique used in this chapte r to compare the sizes of the different stellar 



populations was first presented in 



Leigh. Sills fc Knigge! fl2009f ). In that study 



we introduced our method and applied it to a sample of 56 GCs taken from 



Piotto et al. 



(120021 ) in order to study their RGB populations. Our results were 
consistent with a sub-linear dependence of RGB numbers on the core masses. 
In particular, we found evidence for a surplus of RGB stars relative to MSTO 
stars in the cores of low-mass GCs. We considered several possible causes for 
this result, but concluded that our analysis should ideally be repeated with 
superior photometry in order to properly assess effects such as completeness. 
Given the high-quality of the ACS data, we are now in a position to reassess 
our previous result for RGB stars. 

Upon applying our technique to the ACS sample, we find that the num- 
bers of RGB stars scale linearly with the core masses to within one standard 
deviation. This is also the case for our comparison to the total number of stars 
in the core. This suggests that we should reject our previous c onclusions for 



this stellar population reported in 



Leigh. Sills fc Kniggel (J2009|). Specifically 



if we take a strict 3 — a limit as our criterion for whether or not the slopes are 
sub-linear at a statistically significant level, then the core RGB slope reported 
in this chapter is consistent with being linear whereas this was not the case 
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for the core RGB slope reported in lLeigh. Sills fc Kniggd (120091 ) . However, if 
we take a more stringent criterion for statistical significance, then there is no 
inconsistency between our new and old results for RGB stars, and both slopes 
are consistent with being linear. 

We also tried com paring the RGB catalogue pr esented in this chapter 



with the one presented in 



Leigh. Sills fc Kniggel (12009[ ). This showed that the 



old RGB numbers are slightly deficient relative to the new numbers at the 
high-mass end. Although this difference is not sufficiently large to completely 
account for the difference in slopes (for the comparison with the core mass) 
found between our new and old RGB catalogues, it works in the right direc- 
tion and is likely a contributing factor. If our uncertainties are also factored 
in, then our new and old slopes agree to within one standard deviation (due 
mainly to the large uncertainty for the new slope). The source of the disagree- 
ment between our old and new RGB catalogues is unclear, and we cannot say 
whether or not incompleteness (in the old data set) is the culprit. The results 
of our artificial star tests have at least confirmed that incompleteness is not 
an issue for our new catalogue, howev er it could certainly have co ntributed 
to the lower RGB numbers reported in iLeigh. Sills fc Kniggel ( 120091 ). Indeed, 



the central cluster density tends to be higher in more massive clusters, which 
should negatively affect completeness. It is not clear, however, why this would 
have affected RGB stars more than MSTO stars in the old data set. Given all 
of these considerations, we feel that our new results show that this issue needs 
to be looked at in more detail before any firm conclusions can be drawn. 

The evidence in favour of RGB numbers scaling linearly with the core 
masses is interesting. For one thing, it suggests that two-body relaxation does 
not significantly affect RGB population size relative to other stellar popula- 
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tions of comparable mass in even the dense central regions of GCs. This is 
not surprising, since two-body relaxation is a long-range effect for which the 
stellar radius plays a negligible role. Second, it suggests that collisions do not 
significantly deplete RGB stars relative to other stellar populations despite 
their much larger radii. This is because the collision rate increases with in- 
creasing cluster mass, so we would expect RGB stars to appear preferentially 



depleted in mas s ive c 



Beers fc Davies 



2004 



usters if they are significantly affect ed by collisions (e.g. 



Davies. Piotto fc De Angeli . 



20041 ) . Third, it suggests 



that the sub-linear relation found for BSs does not contribute to a sub-linear 
relation for RGB stars despite the fact that B Ss should eventually evolve t o 



Leigh. Sills fc Knigge! (J2009|). 



occupy our RGB selection box, as discussed in 
This is likely the result of the relatively small sizes of the BS populations in 
our sample when compared to the numbers of RGB stars since the rate at 
which evolved BSs ascend the R GB is thought to be comparable t o the RGB 



lifetimes of regular MSTO stars (jSills. Karakas fc Lattanzio 



20091 ). Alterna- 



tively, this result could, at least in part, be explained if a smaller fraction of 
BSs end up sufficiently bright and blue to be identified as BSs in the CMDs 
of massive GCs. In other words, it could be that a larger fraction of BSs are 
hidden along the MS in massive clusters, as discussed in Section 14.4.11 In 
this case, the contributions to RGB populations from evolved BSs could be 
comparable in all clusters, in which linear relationship between RGB 

numbers and the core masses would be expected. Finally, evolved BSs would 
be expected to have a negligible impact on RGB population size if the average 
BS lifetime is considerably longer than the lifetimes of RGB stars. This effect 
is difficult to quantify , however, given that BS lifetimes are poorly constrain ed 



in the literature (e.g. 



Sandquist. Bolte fc Hernquist 



1997 



Sills et al. 



200 lh . 
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4.4.3 Horizontal Branch Stars 

Our results suggest that HB numbers scale linearly with the core masses. This 
can be interpreted as evidence that two-body relaxation does not significantly 
affect the radial distributions of HB stars in GCs relative to the other stellar 
populations above the MSTO. One reason to perhaps expect that two-body 
relaxation should affect the spatial distributions of HB stars stems from the 
fact that RGB stars are among the most massive stars in clusters, and they 
undergo significant mass loss upon evolving into HB stars. Consequently, 
the progenitors of HB stars should be heavily mass segregated. HB stars 
themselves, however, have relatively low-masses so that two-body relaxation 
and strong dynamical encounters should act to re-distribute them to wider 
orbits within the cluster potential. The HB lifetime is roughly constant at 10 8 



years (jlben 



199ll ) and it is comparable to or exce eds the core relaxation time s 



1996. 2010 updatel ). 



for most of the low-mass clusters in our sample (lHarris 
Therefore, we might expect the HB populations in these clusters to exhibit 
more extended radial profiles relative to more massive clusters. This would 
contribute to a sub-linear relationship between the numbers of HB stars and 
the core masses. Our uncertainties are sufficiently large that this possibility 
cannot be entirely ruled out, however our results are consistent with a linear 
relationship between HB numbers and the core masses. 

4.4.4 Additional Considerations 

Recent observations have reveal ed the presence of mu ltiple stellar populations 



in a number of MW GCs (e.g. 



Pancino et al. 



20031 ). The majority of these 



cases have been reported in very massive clusters. Moreover, their existence 
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is thought to be related to the chemical properties of GCs, in particular an 
observed anti-correlation between their sodium and oxygen abundances. In 
turn, these chemical signa tures have been argu ed to be linked to the cluster 



metallicity, mass and age ( iCarretta et al. 



20101). 



We identified clusters in our sample currently known to host multiple 
stellar populations, but none of these were clear outliers in our plots. Conse- 
quently, the effects had on our results by multiple stellar populations remains 
unclear. It is certainly possible that multiple stellar populations have con- 
tributed to the uncertainties for the weighted lines of best fit performed for 
the relations in Figure 14.31 and Figure 14.41 It is difficult to quantify the possi- 
ble severity of this effect, however, given the limited evidence linking multiple 
stellar populations to cluster properties. 

Although the uncertainties are sufficiently large that the slopes are con- 
sistent with being linear at the 3 — a confidence level for all stellar populations 
when performing the comparisons with the total stellar mass, the reported 
slopes are typically less than unity within the 1 — a, and often even the 2 — a, 
confidence interval. This does not appear to result from the fact that we have 
obtained our estimates for the total stellar masses by numerically integrat- 
ing 3-dimensional density distributions and are comparing to number counts, 
which are projected quantities. To address this, we also tried obtaining the 
total stellar masses by numerically integrating 2-dimensional surface bright- 
ness profiles so that we are consistently comparing only projected quantities. 
Despite this, our results remain unchanged and the new slopes agree with 
the old ones to within one standard deviation for all stellar populations. An- 
other possibility to account for this trend that is perhaps worth considering 
is a systematic dependence of the mass-to-light ratios of clusters on their to- 
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tal mass. There are two ways this could have affected our analysis. First 
stellar remnants have been shown to affect th e dynamical evolution of clus- 
ters, and therefore the sizes o f their cores (e.g 



Trenti. Vesperini &: Pasquato 



Lee. Fahlman fc Richer 



1991 



20101 ) . It follows that, if the number of stellar 
remnants depends on the cluster mass, then this could contribute to an addi- 
tional underlying dependence of the core radius on the cluster mass. This could 
perhaps arise as a result of the fact that the ratio of the rate of stellar evolu- 
tion to the rate of dynamical evolution is larger in more massive clusters, since 
in general the rate of two-body relaxation decreases with increasing cluster 
mass whereas the rate of stellar evolution is independent of the cluster mass. 
Coupled with their deeper gravitational potential wells, this could contribute 
to more massive clusters retaining more stellar remnants. Second, variations 
in the mass-to-light ratios of clusters can also occur as a result of changes in 



the a verage stellar mass (not including stellar remnants) (IKruijssen &: Mieske 



20091 ). That is, we can approximate the total stellar mass contained in the 



core as: 

M core rsj ^Ln—poT* ~ mN core , (4.2) 

where M/L is the mass-to-light ratio, po is the central luminosity density, m 
is the average stellar mass and N core is the total number of stars in the core. 
Based on our results, M core oc N^ e , where we have used the total number of 
stars in the core with magnitudes brighter than 0.5 mag below the MSTO as a 
proxy for N core . This could suggest that m oc iV~°g oc M"^ 1 . In other words, 
the average stellar mass in the core decreases weakly with increasing core 
mass. This could in part be due to the fact that more massive clusters should 
be less dynamically evolved than their less massive counterparts, and should 
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therefore be less depleted of their lo w-mass stars due to s t ellar evaporation 



induced by two-body relaxation (e.g. 


Ambartsumian . 


1938 


1958: 


HenoJ. 


1960; 


De Marchi. Paresce & Porteries Zwart. 



Spitzer fc Harm 



2010). Similarly, 

usters, 



mass segregation should also tend to operate more rapidly in low - mass c 
which acts to migrate preferentially massive stars into the core (ISpitzerl. 



Spitzer fc Hart 



1995 



1971 



Farouki &: Salpeter 



Meylan fc Heggie 



1982 



Shara et al. 



1995 



19691: 



King. Sosin fc Cool 



19971 ). Alternatively, differences in the stellar mass 
function in the core could result from variations in the degree of primordial 
mass segregation, or even variations in the initial stellar mass function. 

We have assumed throughout our analysis that the core mass is a suit- 
able proxy for the total cluster mass. We have checked that these two quanti- 
ties are indeed correlated, however this does not tell the whole story since we 
are also using the total cluster mass as a proxy for the degree of dynamical 
evolution. The central concentration parameter, defined as the logarithm of 
the ratio of the tidal to core radii, describes the degree to which a cluster is 
centrally concentrated. Previous studies have shown that there exists a weak 



corre 
(e.g. 



ation between the concentration parameter an d the total cluster mass 



Djorgovski fc Meylan . 



1994 



McLaughlin 



2000|) . In order to better use 



our technique to reliably probe the effects of the cluster dynamics on the sizes 
and radial distributions of the different stellar populations, the concentration 
parameter should ideally be accounted for when applying our normalization 
technique in future studies. It is not yet clear how the concentration parame- 
ter can be properly absorbed into the normalization, however its effect on our 
analysis should be small given the weak dependence on cluster mass. 

The assumption that the degree of dynamical evolution experienced 
by a given cluster depends only on its mass is also incorrect. Two-body re- 
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laxation has been shown to dominate cluster evolution for a signifi cant frac- 



tion of the lives of old MW GCs (e.g. 



Gieles. Heggie fc Zhao 



20111 ). however 



other effects can also play a significant role. For example, stellar evolution 
is known to affect the dynamical evolution of star clusters, altho ugh its pri 



mary 



1986 



role is played during their 



Chernoff &: Weinberg 



1990 



early evolutionary phases 



Fukushige &: Heggie 



e.g. 



Applegate 



19951). Tidal effects 



from the Galaxy have also been shown to play an important role in deciding 
the dynamical fates of clusters by i ncreas ing the rate of mass loss across the 



tidal boundary (e.g. 



Heggie fc Hut 



20031 1 . Consequently, clusters with small 



perigalacticon distances should appear more dynamically evolved than their 
total mass alone would suggest. This effect can be significant, and has likely 
contributed to increasing the uncertainties found for the comparisons to the 
total stellar mass. Therefore, tidal effects from the Galaxy should also ideally 
be absorbed into our normalization technique in future studies. This can be 
done by using the perigalacticon distances of clusters as a rough proxy for 
the degree to which t ides from the Galaxy should have affected their internal 



dynamical evolution (IGieles. Heggie fc Zhao 



201l|). 



Interestingly, tides could also help to explain why the uncertainties for 
the comparisons to the total number of stars in each circle are considerably 
smaller than for the comparisons to the total stellar mass. We have used the 
number of stars with magnitudes brighter than 0.5 mag below the MSTO as 
a proxy for the total number of stars. Consequently, we are comparing stars 
within a very narrow mass range, so that all populations of interest should have 
been comparably affected by two-body relaxation (except, perhaps, for HB 
stars) independent of tidal effects from the Galaxy. In other words, tides should 
affect all stars above the MSTO more or less equally, and this is consistent with 
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our results. It is also worth mentioning here that our King models consider 
only a single stellar mass. This assumption is not strictly true and could also 
be contributing to increasing the uncertainties found for the comparisons to 
the total stellar masses. 

An additional concern is that we do not know if the clusters in our 
sample are currently in a phase of core contraction or expansion. This has a 
direct bearing on the recent history of the stellar density in the core, and there- 
fore the degree to which stars in the core should have been affected by close 
dynamical interactions. These effects are independent of two-body relaxation 
and occur on a time-scale that is ty pically much shorter than the half-mass 



relaxation time (IHeggie fc Hut 



20031 ) . The effects could be significant in clus- 
ters that were recently in a phase of core-collapse but have since rebounded 
back out of this highly concentrated state. This could occur, for example, as a 
result of binary formation induced by 3-body inter actions com b ined with their 



subse quent hardening via additional encounters (IHut 



1983 



Heggie fc Hut 



20031 ). In general, binaries play an important role in the dynamical evolution 
of clusters, and could have affected our results in a number of ways. This is a 
difficult issue to address even qualitatively given how little is currently known 
about the binary populations in globular clusters. Theoretical models suggest, 
however, that the time-scale for core contraction is often longer than a Hub- 
ble time, and that this evolutionary phase will only come to an end once the 
central density becomes suf ficiently high for hardening enco unters involving 



binaries to halt the process ( iFregeau. Ivanova Sz Rasio 



20091). It follows that 



the cores of most MW GCs are expected to currently be in a phase of core 
contraction. This process is ultimately driven by two-body relaxation, so that 
our assumption that the total cluster mass provides a suitable proxy for the 
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degree of dynamical evolution is still valid. 

In summary, our results suggest that effects related to the cluster dy- 
namics do not significantly affect the relative sizes of the different stellar pop- 
ulations above the MSTO. This is the case for at least RGB, HB and MSTO 
stars. BSs, on the other hand, show evidence for a sub-linear dependence of 
population size on the total stellar mass contained within the same radius. 
Whether or not the cluster dynamics is responsible for this sub-linearity is 
still not clear. Notwithstanding, our results have provided evidence that mass 
segregation is not the dominant cause for this result, although it will be nec- 
essary to redo the comparison performed in this study with a larger spatial 
coverage in order to fully address this question. Further insight into the origin 
of the sub-linearity found for BSs will be provided by reliable binary fractions 
for the clusters in our sample, which are forthcoming (Sarajedini 2010, private 
communication) . 
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Chapter 



5 



An Analytic Model for Blue 
Straggler Formation in Globular 



Clusters 



5.1 Introduction 



Commonly found in both open and globular clusters (GCs), blue stragglers 
(BSs) appear as an extension of the main-sequence (MS) in cluster colour- 
magnitude diagrams (CMDs), occupying the regio n that is j ust b righter and 



bluer than the main-sequence turn-off (MSTO) fjSandage 



1953|). BSs are 



thou ght to be produc e d via the additi o n of hydrogen to low-mass MS stars 



(e.g. 



Sills et al 



2001 



Lombardi et al 



20021 ). This can occur via multiple 



channels, most of which involve the mergers of low-mass MS stars since a signif- 
icant amount of mas s is typically required to reproduce the observed locations 



of BSs in CMDs (e.g. 



Sills fc Bailyn 



19991 ). Stars in close binaries can merge if 



enough orbital angular momentum is lost, which can be mediated by dynam- 
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ical interactions with other stars, mag netized stellar winds, tidal dissipation 



or even an outer triple companion (e.g. 


Leonard & Linnell, 


1992 




Li & Zhane. 


2006; 


Perets & Fabrvckv. 


2009; 


Dervisoglu. Tout & Ibanoglu 


2010 


). Alterna- 



tively, MS stars can collide directly, alt hough this is a 



mediated by multip 



Fregeau et al. 



2001 



e star systems ( e.g. 



Leigh fc Sills 



Leonard 



1989 



so thought to usua 



Leonard &; Livio . 



201ll ). First proposed by 



McCreal ( 119641 ) 



ly be 



1995 



BSs have also been hypothesized to form by mass-transfer from an evolving 
primary onto a normal MS companion via Roche lobe overflow. 

Despite numerous formation mechanisms having been proposed, a satis- 
factory explanation to account for the presence of BSs in star clusters eludes us 
still. Whatever the dominant BS formation mechanism(s) operating in dense 
clusters, it is now thought to somehow involve multiple star systems. This was 



shown to be the case in even the dense cores of GCs (ILeigh. Sills fc Knigge 



2007 



2008 



Knigge. Leigh &: Sills 



are th ought to occur frequently ([Leonard 



200 9J) where co llision s between single stars 



1989|). In 



Knigge. Leigh &; Sills 



(120091 ). we showed that the numbers of BSs in the cores of a large sample of 
GCs correlate with the core masses. We argued that our results are consistent 
with what is expected if BSs are descended from binary stars since this would 
imply a dependence of the form Nbs ~ fbM core , where Nbs is the number of 
BSs in the core, fb is the binary fracti on in the core and M rnr p is the total 



Mathieu fc Gellerl (j2009[ ) also showed 



stellar mass contained within the core, 
that at least 76% of the BSs in the old open cluster NGC 188 have binary 
companions. Although the nature of these companions remains unknown, it 
is clear that binaries played a role in the formation of these BSs. 

Blue stragglers are typically concentrated in the dense cores of glob- 
ular clusters where the high stellar densities should result in a higher rate 
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of stellar encounters (e.g. iLeonardl . Il989f ). Whether or not this fact is di- 
rectly related to BS formation remains unclear, since mass s egregation also 



acts t o migrate BSs (or their p rogenitors) into the core (e.g. 



Saviane et al 



1998 



Guhathakurta et al. 



19981 ). Additionally, num erous BSs have been ob 



Andrievsky et al. 



20001) 



served in more sparsely populated open clusters (e.g. 
and the fields of GCs where collisions are much less likely to occur and mass- 
tran sfer within binary sys tems is thought to be a more likely formation scenario 



(e.g. 



Mapelli et al 



20041 ) . 



Several studies hav e provided evidence 



distribution in some GCs (IFerraro et al. 



; hat BSs show a bimodal sp atial 



1997 



1999 



Lanzoni et al. 



20071 ). In 



these clusters, the BS numbers are the highest in the central cluster regions 
and decrease with increasing distance from the cluster centre until a second 
rise occurs in the cluster outskirts. This drop in BS numbers at intermediate 
cluster radii is often referred to as the "zone of avoidance" . Some authors have 
suggested that it is the result of two separate formation mechanisms occurring 
in the inner and outer regions of the cluster, with mass-transfer in primordial 



binaries dominating i n the latter and stel 



(IFerraro et al 



2004 



Mapelli et al 



ar collisions dominating in the former 



20061 ). Conversely, mass segregation could 



also give rise to a "zone of avoidance" for BSs if the time-scale for dynamical 
friction exceeds the average BS lifetime in only the ou tskirts of GCs that 



exhibit this radial trend (e.g. 



Leigh. Sills fc Knigg 



201ll ). 



Dynamical interactions occur frequently enough in dense clusters that 



they are exp ected to 
of BSs (e.g. IStryker 



pe at l east partly responsib le for the observed properties 



1993 



Leigh fc Sills 



201ll ). It follows that the current 



properties of BS populations should reflect the dynamical histories of their 
host clusters. As a result, BSs could provide an indirect means of probing the 
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physical processes t hat drive star clust er evolution (e.g. 



Hurley et al. 



2005 



Leigh fc Sills 



201lh . 



Heggie fc Hut 



2003 



In this chapter, our goal is to constrain the dominant BS formation 
mechanism(s) operating in the dense cores of GCs by analyzing the principal 
processes thought to influence their production. To this end, we use an analytic 
treatment to obtain predictions for the number of BSs expected to be found 
within one core radius of the cluster centre at the current cluster age. Predicted 
numbers for the core are calculated for a range of free parameters, and then 
compared to the observed numbers in order to find the best-fitting model 
parameters. In this way, we are able to quantify the degree to which each of 
the considered formation mechanisms should contribute to the total predicted 
numbers in order to best reproduce the observations. 



In Section 15.21 we describe the BS catalogue used for comparison to 



our model predictions. In Section 15. 3[ we present our analytic model for BS 
formation as well as the statistical technique we have developed to compare its 
predictions to the observations. These predictions are then compared to the 
observations in Section 15.41 for a range of model parameters. In Section 15. 5[ 
we discuss the implications of our results for BS formation, as well as the 
role played by the cluster dynamics in shaping the current properties of BS 
populations. 



5.2 The Data 



Leigh. Sills fc Knigge! (120111 ). In 



The data used in this study was taken from 
that chapter, we presented a catalogue for blue straggler, red giant branch 
(RGB), horizontal branch (HB) and main-sequence turn-off stars obtained 
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from the colour-magnitude diagra ms of 35 Milky Way G Cs taken from the 



ACS Survey for Globular Clusters (jSarajedini et al. 



2003). The ACS Survey 



provides unprecedented deep photometry in the F606W (~ V) and F814W 
(~ I) filters that extends reliably from the HB all the way down to about 
7 magnitudes below the MSTO. The clusters in our sample span a range 



of tota 



( [Harris 



masses (by nearly 3 orders of magnitude) and central concentrations 



1996. 2010 update! ) . We have confirmed that the photometry is nearly 



complete in the BS region of the CMD for every cluster in o ur sample. This 



was d one using the results of artificial star tests taken from 



Anderson et al. 



( 12008h . 



Each cluster was centred in the ACS field, which extends out to sev- 
eral core radii from the cluster centr e in most of the clusters in our sample. 
Only the core populations provided in lLeigh. Sills fc Kniggd ( 120111 ) are used in 



this chapter. We have taken estimates for the co re radii and central luminos 



Harris ( 



1996. 2010 update!) . 



Webbinki (09851). Esti 



ity densities for the clusters in our sample from 
whereas central velocity dispersions were taken from 
mates for the total stellar mass containe d within the core were obtain ed from 
single- mass King models, as described in iLeigh. Sills fc Kniggd ( 1201 ll ). All of 



the clusters in our sample were chosen to be non-post-core collapse, and have 
surface brightness profiles that provide good fits to our King models. 



5.3 Method 

In this section, we present our model and outline our assumptions. We also 
present the statistical technique used to compare the observed number counts 
to our model predictions in order to identify the best-fitting model parameters. 
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5.3.1 Model 



Consider a GC core that is home to Nbs,o BSs at some time t = t . At a 
specified time in the future, the number of BSs in the core can be approximated 
by: 



N BS = N BS ,o + N coll + N bm + N in - N out - N e . 



(5-1) 



where N co ^ is the number of BSs formed from collisions during single-single 
(1+1), single-binary (1+2) and binary-binary (2+2) encounters, Nf,^ is the 
number formed from binary evolution (either partial mass-transfer between 
the binary components or their complete coalescence), N in is the number of 
BSs that migrate into the core due to dynamical friction, N out is the number 
that migrate out of the core via kicks experienced during dynamical encounters, 
and N e „ is the number of BSs that have evolved away from being brighter and 
bluer than the MSTO in the cluster CMD due to stellar evolution. 

We adopt an average stellar mass of m — O.65M and an average BS 
mass of rriBs = 2m = 1.3M Q . The mass of a BS can provide a rough guide 
to its lifetime, altho ugh a range of lifetimes are s t ill po ssible for any given 



mass. For instance, 



Sandquist. Bolte &: Hernquist 



(119971 ) showed that a 1.3 



M Q blue straggler will have a lifetime of around 0.78 Gyrs in unmixed mod- 
els, or 1.5 7 Gyrs in completely m i xed m odel s. Combined with the r esults of 



Sills et al 



fl200l[ ) 



Lombardi et al. 



((2OO2J) and 



Glebbeek fc Polsl (120081 ). we ex- 



pect a lifetime in the range 1-5 Gyrs for a 1.3 M BS. As a first approximation, 
we choose a like ly value of tbs = 1-5 Gyrs for the average BS lifetime (e.g. 



Sills et al. 



20011 ). The effects had on our results by changing our assumption 
for the average BS lifetime will be explored in Section 15.41 and discussed in 
Section 15.51 
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We consider only the last tbs years. This is because we are comparing 
our model predictions to current observations of BS populations, so that we 
are only concerned with those BSs formed within the last few Gyrs. Any BSs 
formed before this would have evolved away from being brighter and bluer 
than the MSTO by the current cluster age. Consequently, we set Nbs,o = N e „ 
in Equation 15.11 We further assume that all central cluster parameters have 
not changed in the last tbs years, including the central velocity dispersion, 
the central luminosity density, the core radius and the core binary fraction. It 
follows that the rate of BS formation is constant for the time-scale of interest. 
This time-scale is comparable to the half-mass relaxation time but much longer 
than the central relaxation time for the majority of the clusters in our sample 



(Harris 



1996. 2010 updatel ). This suggests that core parameters such as the 
central density and the core radius will typically change in a time tbs since 
the time-scale on wh ich these parameters vary is the central relaxation time 



(IHeggie fc Hut 



20031 ). Therefore, our assumption of constant rates and cluster 
parameters is not strictly correct, however it provides a suitable starting point 
for our model. We will discuss the implications of our assumption of time- 
independent cluster properties and rates in Section 15.51 

In the following sections, we discuss each of the remaining terms in 
Equation 15.11 



5.3.1.1 Stellar Collisions 

We can approximate the number of BSs formed in the last tbs years from 
collisions during dynamical encounters as: 

Ncou = fi+iN 1+1 + f 1+2 N 1+2 + / 2+2 iV 2+2 , (5.2) 
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where Ni+i, Ni + 2 and N2+2 are the number of single- single, single-binary and 
binary-binary encounters, respectively. The terms fi+i, fi+2 and f2+2 are the 
fraction of 1 + 1, 1+2 and 2+2 encounters, respectively, that will produce a BS 
in the last tbs years. We treat these three variables as free parameters since 
we do not know what fraction of collision products will produce BSs (i.e. stars 
with an appropriate combination of colour and brightness to end up in the BS 
region of the CMD), nor do we know what fraction of 1+2 and 2+2 encounters 
will result in a stellar collision. Numerical scattering experiments have been 



performed to study the outcomes of 1+2 and 2+2 encounters (e.g. 



Hut 



1983 



McMillan 



1986 



Fregeau et al. 



20041 ). however a large fraction of the relevant 
parameter space has yet to be explored. 

In terms of the core radius r c (in parsecs) , the central number density 
no (in pc -3 ), the root-mean-square velocity v m (in km s _1 ), the average stellar 
mass m (in M Q ) and the average stellar radius R (in R ), the mean ti me-scale 



between single-single collisions in the core of a GC is ( iLeonardl . 



1989|): 



r 1+1 = l.lxl0 10 (l-/ 6 ) 



\5km/s 



0.5M m \ /0.5-Ro 



m 



R 



10 3 pc 
n 



years 



-3 s 2 



(5.3) 



The additional factor (1-f^) - 2 comes from the fact that we are only considering 
interactions between single stars and the central number density of single stars 
is given by (l-f&)no, where f& is the binary fraction in the core (i.e. the fraction 
of objects that are binaries). For our chosen mass, we assu me a corresp onding 



average stellar radius using the relation M/M = R/R (jlben 



199 lh . The 



number of 1+1 collisions expected to have occurred in the last tbs years is 
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then approximated by: 



t~bs 

n+i 



(5.4) 



The rate of collisions between single stars and binaries, as well as 
between two binary pairs, can be roughly approximated in the same way 



as for single - single encounters ([Leonard 



Bacon et al. 



1996 



Fregeau et al. 



1989 



Sigurdsson fc Phinney 



1993 



20041 ) . We adopt the time-scales derived in 



Leigh fc Sillsl (120111 ) for the average times between 1+2 and 2+2 encounters. 



These are: 



^,4 Xl 0^W,)-^) 3 (^ 



3 x 2 



\5km/s 



O.5M \ flAU 



m 



(5.5) 



years 



and 



?2+2 = 1 

(- 



0.5M G) \ (\AU 



years, 



(5.6) 



\hkmjs) V m ) \ a 
where a is the average binary semi-major axis in the core in AU and we have 
assumed that the average binary mass is equal to twice the average single star 
mass. The numbers of 1+2 and 2+2 encounters expected to have occurred in 
the last r BS years are given by, respectively: 



N- 



1+2 



TBS 



(5.7) 



and 



N- 



2+2 



TBS 
T~2+2 



(5.1 



The outcomes of 1+2 and 2+2 encounters will ultimately contribute 



165 



Ph.D. Thesis — N. Leigh 



McMaster - Physics & Astronomy 



to the evolution of the binary fraction in the core. How and with what fre- 
quency binary hardening/softening as well as capture, exchange and ionization 



interactions affect the b inary fraction in the dense cores of GCs is current 



a subject of debate (e.g. 



Ivanova et al. 



2005 



Hurley. Aarseth fc Shara 



20071). 



y 



Observations are also lacking for binary fractions in the dense cores of GCs, 
however rough c onstraints suggest tha t they range from a few to a few tens 



of a percent (e.g. 



2008 



Davis et al 



Rubenstein fc Bailyn 



1997 



Cool et al. 



2002 



Sollima et al. 



20081 ). The situation is even worse for the distribution of 



binary orbital parameters observed in dense stellar environments. Our best 



constraints come from radial ve 



ters (ILatham 



2005 



Geller at al. 



ocity 



surveys of moderately dense open clus- 



20091 ) . however whether or not the properties 



of the binary populations in these clusters should differ significantly from those 
in the much denser cores of GCs is unclear. As an initial assumption, we as- 
sume a time- independent core binary fraction of 10% for all clusters, and an 
average semi-major axis of 2 AU. This semi-major axis corresponds roughly 
to the hard- soft boundary for most of the clusters in our sample, defined by 
setting the average binary orb ital energy equal to the kinetic energy of an 



average star in the cluster (e.g. 



Heggie fc Hut 



20031 ). We treat both the core 



binary fraction and the average binary semi-major axis as free parameters, 
and explore a range of possibilities using the available observations as a guide 
for realistic values. We will return to these assumptions in Section 15.51 

5.3.1.2 Binary Star Evolution 

Although we do not know the rate of BS formation from binary star evolution, 
we expect a general relation of the form = TBsl T mt for the number of 
BSs produced from binary mergers in the last r^s years, where r mf is the 
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average time between BS formation events due to binary star evolution. We 
can express the number of BSs formed from binary star evolution in the last 

Tbs years as: 



Nbin — fmtfbNc 



(5.9) 



where N core is the total number of objects (i.e. single and binary stars) in the 
core and i m t is the fraction of binary stars that evolved internally to form BSs 
within the last tbs years. We treat f m t as a free parameter since it is likely to 
depend on the mass-ratio, period and eccentricity distributions characteristic 
of the binary populations of evolved GC cores, for which data is scarce at best. 



5.3.1.3 Migration Into and Out of the Core 

Due to the relatively small sizes of the BS populations considered, the mi- 
gration of BSs into or out of the core is an important consideration when 
calculating the predicted numbers. In other words, we are dealing with rel- 
atively small number statistics and every blue straggler counts. In order to 
approximate the number of stars in the core as a function of time, two compet- 
ing effects need to be taken into account: (1) mass stratification/segregation 
(or, equivalently, dynamical friction) and (2) kicks experienced during dynam- 
ical interactions. These effects are accounted for with the variables N,-„ and 



N on4 in Equation 15. 1\ respectively. 

Blue stragglers are among the most massive s tars in clusters (e.g. 



1997 



van den Berg et al. 



2001 



Mathieu fc Geller 



Shara. Saffer fc Livio 



20091 ). so they should typ- 
ically be heayily mass segregated relative to other stellar populations (e.g. 



Spitzer 



1969 



Shara et al. 



1995 



King. Sosin &: Cool . 



19951 ) . The time-scale 



for this process to occur can be approximated using the dynamical friction 
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time-scale ( IBinney fc Tremaind . 119871 ) : 



0-r 



Td 



4/nAG 2 (2vr) 1 /2 m BS p{r) 



(5.10) 



where a(r) and p(r) are, respectively, the velocity dispersion and stellar mass 
density at the given distance from th e cluster centre r. Both a(r) a nd p(r) are 



found from single-mass King models (ISigurdsson fc Phinney 



fit to e ach cluster using the concentration parameters provided in 



19931). which are 



McLaughlin fc van den Marel 



( 120051 ) . The Coulomb logarithm is de noted by A, and we adopt a valu e of 



InA ~ 10 throughout this chapter (e.g. 



Spitzer 



1987 



Heggie fc Hut 



2003|). If 



T~dyn > t~bs at a given distance from the cluster centre, then any BSs formed 
at this radius in the last tbs years will not have had sufficient time to migrate 
into the core by the current cluster age. The maximum radius r max at which 
BSs can have formed in the last tbs years and still have time to migrate into 
the core via dynamical friction is given by setting Td yn = tbs- Therefore, N in 
depends only on the number of BSs formed in the last tbs years at a distance 
from the cluster centre smaller than r max . 

In order to estimate the contribution to Nbs m Equation 15. II from BSs 
formed outside the core, we calculate the number of BSs formed in radial shells 
between the cluster centre and v max . Each shell is taken to be one core radius 
thick, and we calculate the contribution from each formation mechanism in 
every shell. This is done by assuming a constant (average) density and veloc- 
ity dispersion in each shell. Specifically, we estimated the density and velocity 
dispersion at the half-way point in each shell using our single-mass King mod- 
els, and used these to set average values. The number of BSs expected to have 
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migrated into the core within the last tbs years can be written: 



N 

N in = ^ (fl+lN(l+l),i + /l+2-/V(l+2),i + /2+2A r (2+2),' 
i=2 



(5.1i; 



+f mt N { 



(bin),i 



X 1 



T (dyn),i 



where i — 1 refers to the core, i = 2 refers to the shell immediately outside the 
core, etc. and N is taken to be the integer nearest to r max /r c . We let the terms 
with N(i + i) )i , N(i +2 ),i, N( 2 +2),j and represent the number of BSs formed 

in shell % from single-single collisions, single-binary collisions, binary-binary 
collisions and binary star evolution, respectively. The time-scale for dynamical 
friction in shell i is denoted by T(d yn ),i, and the factor (1 - T(dyn),i/ T Bs) is 
included to account for the fact that we are assuming a constant formation 
rate for BSs, so that not every BS formed in shells outside the core will have 
sufficient time to fall in by the current cluster age. 

It is typically the least mass ive stars that are ejected fro m 1+2 and 



2+2 interactions as single stars (e.g. 



Sigurdsson fc Phinney 



19931). Combined 



with conservation of momentum, this suggests that BSs are the least likely to 
be ejected from dynamical encounters with velocities higher than the central 
velocity dispersion due to their large masses. This h as been confirmed by sev- 



eral s tudies of numerical scattering experiments (e.g. 



Hut 



1983 



Fregeau et al. 



20041 ). Based on this, we expect that N out should be very small and SO ; cLS cL 
first approximation, we take N out = 0. However, we also explore the effects of 
a non-zero N out by assigning a kick velocity to all BSs at birth. If dynamical 
interactions play a role in BS formation, we might naturally expect BSs to 
be imparted a recoil velocity at birth (or shortly before) due to momentum 
conservation. We will return to this assumption in Section 15.41 
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5.3.2 Statistical Comparison with Observations 

Our model contains 4 free parameters, which correspond to the fraction of 
outcomes that produce a blue straggler for each formation mechanism (1 + 1 
collisions, 1+2 collisions, 2+2 collisions, and binary star evolution). These are 
the / values described in the previous section: fi + i, /1+2, ^2+2, fmt- We assume 
that these values are constant throughout each cluster, and are also constant 
between clusters. By fitting the predictions of our model to the observational 
data, we can determine best-fit values for each of these / parameters, and 
therefore make predictions about which blue straggler formation processes are 
more important. 

In order to determine the best values for these / param eters, we need an 



appro priate statistical test. For this, we follow the approach of lVerbunt. Pooley fc Bassa 



(120081 ). The number of BSs seen in the core of a globular cluster can be de- 
scribed by Poisson statistics. In particular, the probability of observing iV 
sources when \i are expected is: 



u N 

P(N,») = ^e^ (5.12) 



We can calculate a probability for each cluster, and then calculate an overall 
probability P' for the model by multiplying the individual P values. We can 
then vary the / values to maximize this value. 

In practice, these P values are typically of order ten percent per clus- 
ter, and with 35 clusters, the value of P' quickly becomes extremely small. 
Therefore we chose to work with a modified version of this value: the deviance 
of our model to the saturated model. A saturated model is one in which the 
observed number of sources is exactly equal to the expected number in each 
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cluster. In other words, this is the best that we can possibly do. However, 
because of the nature of Poisson statistics, the probability P of such a model 
(calculated by setting N = \x in Equation 15 . 121) is not equal to 1, but in fact 
has some smaller value. For the numbers of blue stragglers in our clusters, the 
P values for the saturated model run from 0.044 to 0.149, and the value of P' 
is 2.08 x 10" 41 . 

The deviance of any model from the saturated model is given by 

D = 2.0(ln(P;„,) - ln(P; odd )) (5.13) 

The model which minimizes this quantity will be our best-fit model. Ideally, 
the scaled deviance (D/(Ndata — N parame t ers ) should be equal to 1 for a best fit. 
Given the simplicity of our model, we expect that our values will not provide 
this kind of agreement, and we simply look for the model which provides the 
minimum of the scaled deviance. 

5.4 Results 

In this section, we present the results of comparing our model predictions 
to the observations. After presenting the results for a constant core binary 
fraction for all clusters, we explore the implications of adopting a core binary 



fractio n th at depends on the c luster luminosity, as reported in 



Sollima et al. 



(120071 ) and 



Milone et al. 



fl2008h . 
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5.4.1 Initial Assumptions 

The predictions of our model for our initial choice of assumptions are shown 
in Figure 15.11 These numbers are plotted against the total stellar mass in the 
core along with the number of BSs observed in the core (blue triangles). We 
plot both the total number of BSs predicted to have formed within r max in the 
last tbs years (Nbs in Equation 15. lj green circles), as well as the total number 
formed only in the core (N coH + N fe j n ; red circles). Upon comparing ~N BS to the 
observed number of BSs in the core, the best-fitting model parameters predict 
that most BSs are formed from binary star evolution, with a small contribution 
from 2+2 collisions being needed in order to obtain the best possible match 
to the observations. The ideal contribution from 2+2 collisions constitutes at 
most a few percent of the predicted total for most of the clusters in our sample. 
Even for our best-fitting model parameters, our initial choice of assumptions 
predicts too few BSs in clusters with small core masses. 

5.4.2 Binary Fraction 

We tried changing our assumption of a constant f;, for all clusters to one for 
which the core binary fraction varies with the total cluster magnitude. First, 
we adopted a dependence of the form: 

/ 6 = 0.13M V + 1.07, (5.14) 

at ion comes from fitting 



where My is the total cluster V magni tude. This re 
a line of best-fit to the observations of Sollima et al 



(120071 ). who studied the 



binary fractions in a sample of 13 low-density GCs (we calculated an average 
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Figure 5.1 The predicted number of BSs plotted versus the total stellar mass 
in the core for the best-fitting model parameters found for our initial choice 
of assumptions. The filled blue triangles correspond to the observed numbers, 
the filled green circles to the number of BSs predicted to have formed within 
Tmax, and the filled red circles to the number predicted to have formed in only 
the core. The best-fitting model parameters used to calculate the predicted 
numbers are f 1+1 = 0, f 1+2 = 0, f 2+ 2 = 7.3 x 10~ 3 , and f mt = 1.7 x 10~ 3 . 
Estimates for the total stellar mass es in the core were obtained from single- 
mass King models, as described in iLeigh. Sills fc Kniggd (120 111 ). Error bars 
have been indicated for the observed numbers using Poisson statistics. 



of columns 3 and 4 in their Table 3 and used these binary fractions to obtain 
Equation I5.14p . In order to prevent negative binary fractions, we impose a 
minimum binary fraction of f^ m = 0.01. In other words, we set ft, = f^ m if 
Equation 15.141 gives a binary fraction less than f^ m . As before, we adopt an 
average semi- major axis of 2 AU. The results of this comparison are presented 
in Figure 15.21 As in Figure 15.11 both the numbers of observed (blue triangles) 
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and predicted (green circles) BSs in the core are plotted versus the total stellar 
mass in the core. Once again, the predicted numbers include all BSs formed 
within imax in the last tbs years. The best-fitting model parameters for this 
comparison suggest that both single-single collisions and binary star evolution 
are significant contributors to BS formation. Single-single collisions contribute 
up to several tens of a percent of the predicted total in several clusters. We 
obtain a deviance in this case that is significantly larger than that obtained for 
the best-fitting model parameters assuming a constant ff, of 10% for all clusters. 
Equation 15. 141 gives higher binary fractions in low- mass clusters relative to our 
initial assumption of a constant f&. This increases the number of BSs formed 
from binary star evolution and improves the agreement between our model 
predictions an d the observat i ons in low-mass clusters. This is consistent with 
the results of 



Sollima et al. 



(120081 ). However, adopting Equation 15.141 for f& 
also causes our model to under-predict the number of BSs in several high- 
mass clusters. 

The best fit to the observations is found by adopting the relation for 
fft provided in Equation 15.141 and setting f^ m = 0.1 (however we note that a 
comparably good agreement is found with a slightly lower f^ m = 0.05). This 
improves the agreement between our model predictions and the observations 
by increasing the number of BSs formed from binary star evolution in massive 
clusters. The result is a good agreement between our model predictions and 
the observations in both low- and high-mass cores, as shown in Figure 15.31 In 
this case, the best-fitting model parameters yield the lowest deviance of any 
of the assumptions so far considered. These best-fitting values suggest that 
most BSs are formed from binary star evolution, with a small contribution 
from 2+2 collisions being needed in order to obtain the best possible match 
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Figure 5.2 The number of BSs predicted in the cluster core using the binary 
fractions of lSollima et all (120071 ) with f^ m = 0.01 plotted against the total stel- 
lar mass contained within the core. The colours used to indicate the observed 
and predicted numbers are the same as in Figure 15.11 The predicted num- 
bers correspond to the best-fitting model parameters, which are fi + i = 0.41, 
fi+2 = 0, f 2+2 = 0, and f mi = 1.5 x 10~ 3 . 



to the observations. Similarly to what was found for our initial assumptions, 
the ideal contribution from 2+2 collisions constitutes at most a few percent 
of the predicted total for most of the clusters in our sample. On the other 
hand, if we change our imposed minimum binary fraction to f^ m = 0.05 we 
find that a non-negligible (i.e. up to a few tens of a percent) contribution 
from single-single collisions is needed in several clusters to obtain the best 
possible agreement with the observations (which is very nearly as good as 
was found using f^ m = 0.1). All of this shows that, although binary star 
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evolution consistently dominates BS formation in our best-fitting models, at 
least some contribution from collisions (whether it be 1+1 or 2+2 collisions, or 
some combination of 1+1, 1+2 and 2+2 collisions) also consistently improves 
the agreement with the observations. Moreover, it is interesting to note that 
an improved agreement with the observations could alternatively be obtained 
if we keep fj^ m = 0.01 but all or some of fi + i, fi + 2 and f2+2 increase with 
increasing cluster mass. This would also serve to improve the agreement at 
the high-mass end. We will return to this in Section 15.51 
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Figure 5.3 The number o f BSs predicted in the cluster core using the binary 
fractions of ISollima et al.l (120071 ) with f™ n = 0.1 plotted against the total stel- 
lar mass contained within the core. The colours used to indicate the observed 
and predicted numbers are the same as in Figure 15.11 The predicted numbers 
correspond to the best-fitting model parameters, which are fi+i = 0, fi + 2 = 0, 
f 2+2 = 3.6 x 10~ 3 , and f mi = 1.4 x 10~ 3 . 
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Fi nally, we also t ried a dopting; the observed dependence of f& on My re- 
ported in 



Miloneetal 



(120081 ) . who also found evidence for an anti-correlation 
between the core binary fraction and the total cluster mass. Despite this 
change, we consistently find that our results are the same as found when using 
the empirical binary fraction relation provided in Equation 15.141 

5.4.3 Average BS Lifetime 

We also tried changing our assumption for the average BS lifetime. We ex- 
plored a range of plausible lifetimes based on values found throughout the 
literature. Specifically, we explored the range 0.5-5 Gyrs. We find that at the 
low end of this range, our model fits become increasingly poor. This is be- 
cause lower values for tbs correspond to smaller values for r max and decrease 
the term (1 - tuyn),i/ T Bs) i n Equation 15.111 This reduces the contribution 
to the total predicted numbers from BSs formed outside the core that fall 
in via dynamical friction. Conversely, our model fits improve for tbs > 1-5 
Gyrs since this corresponds to a larger contribution to Ngs from N in . It is 
important to note, however, that this same effect can be had by increasing the 
number of BSs formed outside the core, since this would also serve to increase 
Ni n in Equation 15.11 This can be accomplished by, for instance, increasing 
the binary fraction outside the core (which would increase the number of BSs 
formed from binary star evolution outside the core that migrate in due to 
dynamical friction) relative to inside the core. This seems unlikely, however, 
given that observations of low-density globular clusters and open clusters sug- 
gest that their bina ry fractions tend to drop off rapidly outside the core (e.g. 



Sollima et al. 



20071 ) . We will return to this issue in Section 
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The best possible match to the observations is found by adopting an 
average BS lifetime of 5 Gyrs along with the relation for f?, provided in Equa- 
tion l5.14l with f^ m = 0.1. The predictions of our model are shown in Figure l5~4l 
for these best-fitting model parameters. As shown, the agreement between our 
model predictions and the observed numbers is excellent. The best agreement 
is found by adopting an average BS lifetime of 5 Gyrs, however the agreement 
is comparably excellent down to slightly less than tbs ~ 3 Gyrs. Although 
increasing tbs does contribute to improving the agreement between our model 
predictions and the observations, the effect is minor compared to the improve- 
ment that can be found by changing our assumption for the binary fraction. 
This is apparent upon comparing Figure 15.41 to Figure 15.31 



5.4.4 Migration 



In order to explore the sensitivity of our results to our assumption for r max , 
we also tried setting N in equal to the total number of BSs expected to form 
within 10 r c for all clusters. For comparison, for an average BS lifetime of 
tbs = 1.5 Gyrs, r max ranges from 2 - 15 r c for the clusters in our sample. 
Despite implementing this change, our results remained the same. This is 
because the largest contribution to the total number of BSs comes from those 
BSs formed in the first few shells immediately outside the core that migrate 
in due to dynamical friction. 

Several GCs have been reported to show evidence for a decrease in 



their binary 



Sollima et al. 



ractiqns with increasing distance from the cluster centre (e.g. 



2007 



Davis et al 



20081 ). This effect is often significant, with 



binary fractions decreasing by up to a factor of a few within only a few core 
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Figure 5.4 The predicted number of BSs plotted versus the total stellar mass 
in the core for the best-fitting model parameters found using an average BS 
lifetime of 5 Gyrs and the relation for the cluster binary fraction provided 
in Equation 15.141 with f™ n = 0.1. The colours used to indicate the observed 
and predicted numbers are the same as in Figure 15.11 The best-fitting model 
parameters used to calculate the predicted numbers are f i+ i = 0, f i+2 = 0, 
h+2 — 1-6 x 10~ 3 , and f mt = 9.9 x 10 -4 . The agreement with the observations 
is excellent for these best-fitting values. 



radii from the cluster centre. Based on this, our assumption that f& is indepen- 
dent of the distance from the cluster centre likely results in an over-estimate 
of the true binary fraction at large cluster radii. In order to quantify the pos- 
sible implications of this for our results, we tried setting 4 = for all shells 
outside the core. Although this assumption is certainly an under-estimate for 
the true binary fraction outside the core, our results remain the same (albeit 
the agreement with the observations is considerably worse than for most of 
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our previous model assumptions). Once again, the best-fitting model param- 
eters suggest that most BSs are formed from binary star evolution, with a 
non-negligible (i.e. up to a few tens of a percent in some clusters) contribution 
from binary-binary collisions. Our results indicate that, if the binary fraction 
is negligible outside the core, then the contribution from BSs that migrate into 
the core due to dynamical friction is also negligible. This is because the time 
between 1+1 collisions increases rapidly outside the core, and every other BS 
formation mechanism requires binary stars to operate. 

We also explored the effects of assuming a non-zero value for N OMf in 
Equation 15.11 by imparting a constant kick velocity to all BSs at birth. Using 
conservation of energy, we calculated the cluster radius to which BSs should 
be kicked upon formation, and used the time-scale for dynamical friction at 
the kick radius to calculate the fraction of BSs expected to migrate back into 
the core within a time tbs- Regardless of our assumption for the kick velocity, 
this did not improve the deviance for any of our best-fitting model parameters. 

5.4.5 Average Binary Semi-Major Axis 

We investigated the dependence of our results on our assumption for the av- 
erage binary semi-major axis. However, this had a very small effect on our 
results. This is because only N 1+2 and N 2 +2 depend on the average semi-major 
axis, and neither of these terms dominated BS production regardless of our 
model assumptions. Only f 2 +2 is non-zero for our best-fitting models however, 
as before, it consistently suggests that far fewer BSs should be formed from 
2+2 collisions than from binary star evolution. 
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5.5 Summary & Discussion 

In this chapter, we have presented an analytic model to investigate BS forma- 
tion in globular clusters. Our model predicts the number of BSs in the cluster 
core at the current cluster age by estimating the number that should have ei- 
ther formed there from stellar collisions and binary star evolution, or migrated 
in via dynamical friction after forming outside the core. We have compared the 
results of our model predictions for a variety of input p arameters to observed 



BS nu mbers in 35 GCs taken from the catalogue of 



Leigh. Sills &: Knigge 



20111). 

What has our model told us about BS formation in dense cluster en- 
vironments? The agreement between the predictions of our model and the 
observations is excellent if we assume that: 

• Binary star evolution dominates BS formation, however at least some 
contribution from 2+2 collisions (most of which occur in the core) must 
also be included in the total predicted numbers. Although it is clear 
that including a contribution from dynamical encounters gives the best 
possible match to the observations, it is not clear how exactly this is 
accomplished in real star clusters. Does the cluster dynamics increase 
BS numbers via direct collisions? Or do dynamical interactions somehow 
modify primordial binaries to initiate more mass-transfer events? We will 
return to this point below. 

• The binary fraction in the core is inversely correlated with the total cms 



ter lum ino sity similar to the e mpirical relations found by 



Sollima et al. 



(120071 ) and 



Milone et al. 



( 120081 ) . We also require a minimum core binary 
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fraction of 5 — 10%. The inverse dependence of f& on the total clus- 
ter mass contributes to a better agreement with the observations at the 
low-mass end of the distribution of cluster masses, whereas the imposed 
condition that f^ m = 0.05 — 0.1 contributes to improving the agreement 
at the high-mass end. 

• BSs formed outside the core that migrate in by the current cluster age 
contribute to the total predicted numbers. 

• The average BS lifetime is roughly a few (~ 3-5) Gyrs, since this increases 
the fraction of BSs formed outside the core that will have sufficient time 
to migrate in due to dynamical friction. 

Our model can only provide a reasonable fit to the observations for 
all cluster masses if we assume that the cluster binary fraction is inversely 
proportional to the total cluster mass. It is interesting to consider the possi- 
bility that such an inverse proportionality could arise as a result of the fact 
that the r ate of two-bod y relaxation is also inversely proportional to the clus- 
ter mass ( jSpitzerl . Il987l ). Consequently, since binaries tend to be the most 
massive objects in GCs, they should quickly migrate into the core in low- 
mass clusters, contributing to an increase in the core binary fraction over time 



( iFregeau. Ivanova fc Rasio 



20091 ). This process should operate on a consider- 
ably longer time-scale in very massive GCs since the time-scale for two-body 
relaxation is very long. Mass segregation could then contribute to the observed 
sub-linear dependence of BS numbers on the core masses by acting to prefer- 
entially migrate the binary star progenitors of BSs into the cores of low-mass 
clusters. This is one example of how a direct link could arise between the ob- 
served properties of BS populations and the dynamical histories of their host 
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clusters. Although this scenario is interesting to consider, we cannot rule out 
the possibility that an anti-correlation between the core binary fraction and 
the total cluster mass could be a primordial property characteristic of GCs at 
birth. 

When interpreting our results, it is important to bear in mind that 
binary star evolution and dynamical interactions involving binaries may not 
always contribute to BS formation independently. For example, dynamics 
could play an important role in changing the distribution of binary orbital 
parameters so that mass-transfer occurs more commonly in some clusters. 
One way to perhaps compensate for this effect would be to include a factor 
oil/ a (where a is the average binary semi-major axis) in Equation 15.91 This 
would serve to account for the fact that we might naively expect clusters 
populated by more close binaries to be more likely to have a larger fraction 
of their binary populations undergo mass-transfer. This does not, however, 
guarantee that more BSs will form since our poor understanding of binary 
star evolution prevents us from being able to predict the outcomes of these 
mass-transfer events, and whether or not they will form BSs. Moreover, little is 
known about the distribution of orbital parameters characteristic of the binary 
populations in globular clusters, and how they are typically modified by the 
cluster dynamics. For these reasons, the interpretation of our results must be 
done with care in order to ensure that reliable conclusions can be drawn. 

In general, our results suggest that binary stars play a crucial role 
in BS formation in dense GCs. In order to obtain the best possible agree- 
ment with the observations, an enhancement in BS formation from dynamical 
encounters is required in at least some clusters relative to what is expected 
by assuming a simple population of binaries evolving in isolation. It is not 
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clear from our results, however, how exactly this occurs in real star clus- 
ters. Dynamics could enhance BS formation directly by causing stellar col- 
lisions, or this could also occur indirectly if the cluster dynamics somehow 
induces episodes of mass-transfer by reducing the orbital separations of bina- 
ries. But in which clusters is BS formation the most strongly influenced by 
the cluster dynamics? Unfortunately, no clear trends have emerged from our 
analysis that provide a straight-forward answer to this question. However, 
our results are consistent with dynamical interactions playing a more signif- 
icant role in more massive clusters (although this does not imply that the 
cluster dynamics does not also contribute in low-mass clusters). This could 
be due to the fact that more massive c l usters also tend to have higher cen- 



tral densities (e.g. 



Djorgovski fc Meylan . 



19941 ). and therefore higher collision 



ra tes. This picture is, broadly 



of 



spea 



ring, roughly consistent with the results 



Davies. Piotto &: De Angelil ( 120041 ) . These authors considered the observed 



dependen ce (or lack thereof ) of BS numbers on the total cluster masses pre- 



sented in 



Piotto et al. 



(120041 ) . and suggested that primordial binary evolution 
and stellar collisions dominate BS production in low- and high-mass clusters, 
respectively. 

Our model neglects the dynamical evolution of GCs and the result- 
ing changes to their global properties, including the central density, velocity 
dispersion, core radius and binary fraction. As a young cluster evolves, dy- 
namical processes like mass segregation and stellar evaporation tend to result 
in a smaller, denser core. Within a matter of a few half-mass relaxation times, 
a gravothermal instability has set in and the collapse en sues on a time -scale 



determined by the rate of heat flow out of the core (e.g. 



Spitzer 



1983). We 



are focussing on the last tbs years of cluster evolution, a sufficiently late pe- 
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riod in the lives of most GCs that gravothermal collapse will have long since 
taken over as the primary driving force affecting the stellar concentration in 
the core. Most of the GCs in our sample should c urrently be in a phase 



of co re contraction (IFregeau. Ivanova &: Rasio 



2009 



Gieles. Heggie &: Zhao 



201ll ). and their central densities and core radii should have been steadily de- 
creasing over the last tbs years. Therefore, by using the currently observed 
central cluster parameters and assuming that they remained constant over the 
last tbs years, we have effectively calculated upper limits for the encounter 
rates. This could suggest that we have over-estimated the importance of dy- 
namical interactions for BS formation. On the other hand, some theoretical 
models of GC evolution suggest that the hard binary frac tion in the co r e of a 



dense stellar system will g e neral 



Fregeau. Ivanova fc Rasio 



Hurley et al. 



2005 



y increase with time (e.g 
20091 ). This can be understood as an imbalance 
between the migration of binaries into the core via mass segregation and the 
destruction of binaries in the core via both dynamical encounters and their 
internal evolution. This could suggest that our estimate for the number of 
BSs formed from binary star evolution should also be taken as an upper limit. 
The key point is that GC evolution can act to increase the number of BSs in 
the cluster core via several different channels. The effects we have discussed 



shou 



d typically be small, h owever, since tbs is much shorter than the cluster 



age (IDe Angeli et al. 



20051 ) and any changes to global cluster properties that 



occur during this time will often be small. 

Our model adopts the same values for all free parameters in all clusters. 
In particular, this is the case for several global cluster properties, including the 
average stellar mass, the average BS mass, the average BS lifetime, and the av- 
erage binary semi-major axis. With the exception of the average stellar mass, 
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there is no conclusive observational or theoretical evidence to indicate that 
these parameters should differ from cluster-to-cluster, although we cannot rule 
out this possibility. For instance, the distribution of binary orbital parameters 



could depen d on cluster properties 



persion (e.g. 



Sigurdsson fc Phinney 



te th e total mass, density or velocity dis- 



19931 ). In particul ar, the central velocity 



dispe rsion should be higher in more massive GCs (e.g. 



Djorgovski &: Meylan 



19941 ). which should correspond to a smaller binary orbital separation for the 
hard-soft boundary. This could contribute to massive GCs tending to have 
smaller average binary orbital separations sin ce soft binaries shou ld not sur- 



Heggie fc Hut 



2003|). In turn, 



vive for long in the dense cores of GCs (e.g. 
this could affect the occurrence of mass-transfer events, or of mergers during 
1+2 and 2+2 encounters. This last point follows from the fact that numerical 
scattering experiments have shown that the probability of mergers occurring 
during 1+2 a nd 2+2 interac t ions i ncreases with decreasing binary orbital sep- 



aration (e.g. 



Fregeau et al. 



20041 ). Both the average stellar mass and the 
average BS mass (an d hence lifetime) could also d epen d on the total cluster 



mass, as discussed in 



(1201 lh 



Leigh. Sills fc Kniggel ( I2009h and 



Leigh. Sills fc Knigge 



We have also neglected to cons ider the importance of t riples for BS 



formation throughout our analysis (e.g. 



Perets fc Fabrycky 



20091 ) since we are 



unaware of any observations of triples in GCs in the literature. Interestingly, 
however, our results for binary star evolution can be generalized to include 
the internal evolution of triples since they should have the same functional 
dependence on the core mass (i.e. N te oc ftM core , where N te is the number of 
BSs formed from triple star evolution and i t is the fraction of objects that are 
triples). 
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Finally, our model assumes that several parameters remain constant as 
a function of the distance from the cluster centre, including the binary fraction 
and the average semi-major axis. However, observations of G Cs suggest that 



their binary fracti o ns cou ld fall off rapidly outside the core (e.g. 



2007 



Sollima et al. 



Davis et al. 



20081 ). Our results suggest that, if the binary fraction is 
negligible outside the core, then the contribution from BSs that migrate into 
the core due to dynamical friction is also negligible. This is because the time 
between 1+1 collisions increases rapidly outside the core, and every other BS 
formation mechanism requires binary stars to operate. On the other hand, the 
presently observed binary fraction outside the core could be low as a result 
of bi naries having previously migrate d into the core due to dynamical friction 



(e.g. 



Fregeau. Ivanova fc Rasio 



20091 ). If these are the binary star progenitors 
of the BSs currently populating the core, then dynamical friction remains an 
important effect in determining the number of BSs currently populating the 
core. 

Despite all of these simplifying assumptions, we have shown that our 
model can reproduce the observations with remarkable accuracy. Notwith- 
standing, the effects we have discussed could be contributing to cluster-to- 
cluster differences in the observed BS numbers. Our model provides a well- 
suited resource for addressing the role played by these effects, however future 
observations will be needed in order to obtain the desired constraints (e.g. 
binary fractions, distributions of binary orbital parameters, etc.). 
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Chapter 



6 



A New Quantitative Method for 
Comparing the Stellar Mass 
Functions in a Large Sample of 
Star Clusters: Evidence for a 



Universal Initial Mass Function in 



Massive Star Clusters in the Early 
Universe 



6.1 Introduction 



It is now known t hat most, if not all, o 
star clusters (e.g. 



Lada fc Lada 



1995 



the stars in our Galaxy were 



2003 



McKee fc Ostriker 



20071 ) 



born m 
. And 
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yet, there remain several key details of the star formation process that are still 
not understood. Part of the problem lies in the fact that populations of young 
stars are typically hidden by a dense veil of optically-thick gas and dust. This 
prevents the escape of most of the light prod uced by infant stars, and often 



renders these regions difficult 



2005 



; o obs erve (e.g. 



Lada. Alves fc Lombardi 



Grenier. Casandijan fc Terrier 



20071 ) . Most of these cluste rs are 



ulated and are of relatively low mass (M < 10 4 M Q ) (e.g. 



Lada 



sparse ly pop- 



19851 ). They 



are also ver y young since these clusters are unlikely to survive for more than 



1 Gyr (e.g. 



Portegies Zwart. McMillan fc Gieles . 



20101 ). At the other end of 



the mass spectrum, most massive star clusters (M > 10 4 M ) in our Galaxy 
tend to be at least a few Gyrs old, and in many cases are nearly as old as the 



Universe itself (e.g. 



Harris 



1996. 2010 update 



De Angeli et al. 



20051 ). Un- 



fortunately, the co nditions present at t he time of the ir formation have been 



largely erased (e.g. 



Hurley et al. 



2005 



Murray 



20091 ). This presents a con- 



siderable challenge for studying star formation in the regime of cluster masses 
and metallicities that characterize Milky Way globular clusters. These old star 
clusters contain the fossil record of a very early episode of star formation in 
the Universe, and are the only means of studying it locally. 

One of the primary observational tests for star formation theories is 
the stellar initial mass function (IMF). Current observational evidence sug- 
gests that the IMF is very similar in several diff erent regions of ou r Galaxy, 



including the disk and young star clusters (e.g. 



Elmegreenl. 



this is still being debated throughout the literature (e.g. 



Scalo 



1999). h owever 



1998 



argua bly "standard" IMF to come from these observations is that of 



). The 



Kroupa 



(120011 ). who fit a three-part power-law with breaks at 0.08 M and 0.5 M Q . 
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This can be expressed as: 



dN 
dim 



/3m 



(6.1) 



where a is a constant given by a = 2.3 for 0.5 < m/M Q < 50, a = 1.3 for 0.08 
< m/M Q < 0.5, and a = 0.3 for 0.01 < m/M < 0.08, and (3 is a constant 
determined by t he total cluster mass. By considering the mass function up 



to only ~ 1 M G 



Miller fc Scald (119791 ) found a good fit to the observed mass 



distribution using a log-normal functional form: 



dN 
dlnm 



oc exp 



(Inm — lnm c 
2^2 



(6.2) 



where m c ~ 0.2 M Q and a ~ 0.55 ( iChabrier 



20051). 



Different star formation theories tend to predict different IMFs. These 
tend to va ry with the prop e rties of the gas clouds from which the stars are 



born (e.g. 



Elmegreen 



2001 



Bonnell. Larson fc Zinnecker 



20071 ). Given the 



sensitive nature of the observations, a large sample of IMFs spanning the en- 
tire range of cluster properties exhibited by star clusters in the Milky Way, 
including total mass and chemical composition, has yet to be compiled. This 
is a sorely needed step in order to advance our understanding of star for- 
mation by providing direct comparisons for theoretical predictions. This is 
especially true of massive (M > 10 4 M Q ), metal-poor star clusters since we 
are particularly lacking observations of IMFs in thi s regime of cluster masses 



and metallicites, especially in our own Ga laxy (e.g. 



McKee fc Ostriker 



2007 



Portegies Zwart. McMillan fc Gieles 



mm- 



For the very first time, the ACS Survey for Globular Clusters has pro- 
vided photometry that is nearly reliable all the way down to the hydrogen 
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burning limit in a large sample of Milky Way globular clusters. This offers a 
large sample of current stellar mass functions spanning the stellar mass range 
~ 0.2 — 0.8 M . All of the clusters are massive and of very old age, with total 



masses and ages ranging from ~ 10 4 - 10 6 A/ 



([Harris 



1996. 2010 update 



De Angeli et al. 



Sl _and ~ 10-12 Gyrs, respectively 



20051 ). This has allowed signif- 



icant time for their stellar mass functions to have been modified from their 
primordial forms due to both stellar evolution and stellar dynamics. However, 
most of the processes responsible for this evolution are now largely under- 
stood. Therefore, in principle, it is possible to use current observations of old 
star clusters together with theoretical models for their evolution to indirectly 
probe their IMFs. 

For most of the life of a star cluster, two- body relaxation is the domi- 
nant physical mechanism driving its evolution (e.g. 



Heggie fc Hut 



2003 



Gieles. Heggie fc Zhao 



201ll ). The term describes the cumulative effects of long-range gravitational 
interactions that occur between pairs of stars, which act to alter their orbits 
within the cluster. Two-body relaxation also acts to slowly modify the dis- 
tribution of stellar masses within clusters. Among other things, it results in 
a phenomenon known as mass segregation. This is the tendency for heavier 
stars to accumulate in the central cluster regions and low-mass stars to be 
dispersed to wider orbits. This same process also causes the continual escape 
of preferentially low-mass stars from the clust er, which has been confirmed to 



occur observationally in real star clusters (e.g . Ivon Hippel fc Sarajedini 



1998 



De Marchi. Paresce fc Portegies Zwart 



2010) 



The time-scale for two-body relaxation to operate can range anywhere 
from several million years to the age of the Universe or longer. The rate at 
which it occurs can be roughly approximated using the half-mass relaxation 
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time, w 



r ich pr ovides a rough average for the entire cluster. This is given by 



(jSpitzer 



19871): 



t rh = 1.7 x ltf>[r h (pc)f l2 N l/2 [m/M Q }- l/2 yeari 



(6.3) 



where is the half-mass radius (i.e. the radius enclosing half the mass of the 
cluster), N is the total number of stars within and m is the average stellar 
mass. The half-mass radii of MW GCs are remarkably similar independent of 



few over the course of a cluster's lifetime ( iHenon 



f Henon . 


1973; 


Murrav . 


2009) 



GCs that comprise the ACS sample show a range of masses spanning roughly 



3 orde rs of magnitude 



Gyrs) (IDe Angeli et al. 



■0 4 -1 6 M©), and have comparably old ages (~ 10-12 



20051 ). Therefore, Equation [63] suggests that the total 
cluster mass provides a rough proxy for the degree of dynamical evolution (due 
to two-body relaxation). In other words, the effects of two-body relaxation on 
the evolution of the stellar mass function should be the most pronounced in 
the least massive clusters in the ACS sample. Said another way, dynamical 
age increases with decreasing cluster mass. 

In this chapter, we present a new technique to quantify cluster-to- 
cluster variations in the observed stellar mass functions of a large sample 
of clusters spanning a diverse range of properties. Our sample consists of 33 
Milky Way globu l ar clu sters taken from the ACS Survey for Globular Clusters 



( iSarajedini et al 



20071 ) . With it, we constrain the universality of the IMF 



in a large sample of old, metal-poor star clusters spanning a wider range of 



masses than ever before considered. In Section 16.21 we present our sample of 
stellar mass functions and describe our technique. The results of our analysis 
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of the ACS observations are presented in Section 16. 31 Finally, we discuss the 
implications of our results for the universality of the stellar IMF in Section 
and describe how our results can be compared to theoretical models for GC 
evolution. 



6.2 Method 

In this section, we describe how we acquire our sample of mass functions from 
the ACS data. 



6.2.1 The Data 



The data used in this study consists of a sample of 33 MW GCs taken from the 

' ,0 



ACS Survey for Globular Clusters (jSarajedini et al. 



2007M 1 I The ACS Survey 



provides unprecedented deep photometry in the F606W (~ V) and F814W 
(~ I) filters that is nearly complete down to ~ 0.2 M . In other words, the 
colour-magnitude diagrams (CMDs) extend reliably from the HB all the way 
down to about 7 magnitudes below the main-sequence turn-off (MSTO). 

Each cluster was centred in the ACS field, which extends out to several 



core radii from the cluste r centre in most c 
centres were taken from iGoldsbury et al. 



usters. Coordinates for the cluster 



( 120101 ) . These authors found their 



centres by fitting a series of ellipses to the density distributions within the 
inner 2' of the cl uster c e ntre, and computin g an average value. The core radii 



were taken from 



Harris! (11996. 2010 updateh . 



1 The data can be found at http://www.astro.ufl.edu/ ata/publicJistgc/. 
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6.2.2 Mass Bin Selection Criteria 



In order to select the number o f stars belong i ng; to each mass bin, we fit 



theoretical isochrones taken from 



Potter et al. 



(120071 ) to the CMDs of every 



cluster in our sample. Each isochrone was generated using the metallicity 
and age of the cluster, and fit to i ts CMP usin g the corresponding distance 



modulus and extinction provided in 



Potter et al. 



( 120101 ). The MSTO was then 



defined using our isochrone fits by selecting the bluest point along the MS. 

We considered five mass bins along the main-sequence. These ranged 
from 0.25 - 0.75 M in increments of 0.1 M . This range was chosen to ensure 
complete sampling in all bins since the lowest MSTO mass in our sample 
corresponds to ~ 0.75 M , and the photometric errors remain small (< 0.05 
mag) within the corresponding magnitude range for each cluster. We have 
obtained number counts for all mass bins within the core, as well as for within 
two and three core radii. We do not consider circles outside this since the 
spatial coverage becomes incomplete for several clusters. This greatly reduces 
our sample size and causes the statistical significance of our analysis to suffer. 

We have obtained completeness fractions for each mass bin in all three 



annuli for every cluster in our sample. This was done using the results of 

,0 



artificial star tests taken from 



Anderson et al 



(120081 ) J I Number counts for each 



mass bin were then multiplied by their corresponding completeness corrections. 
The field of view of the ACS images is about 200" on a side, which gives 
physical scales ranging between 1.5 and 16 pc (for the closest and furthest 
clusters in our sample). Based on this, we expect foreground contamination 
by field stars to be negligible for most of the clusters in our sample given their 



2 Artificial star tests were obtained directly from Ata Sarajedini via private communica- 



tion. 
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Da Costa] (119821 ) considered star 



current locations in the Galaxy. For example, 
count data in a similar area and over a comparable range of stellar masses for 
three nearby globular clusters. The author found that the corrections resulting 
from field contamination were always less than 10% over nearly the entire range 
of stellar masses we are considering. 



6.2.3 Weighted Lines of Best-Fit 

In order to quantify cluster-to-cluster differences in the current stellar mass 
functions of the clusters in our sample, we have obtained lines of best-fit for 
(the logarithm of) the number of stars belonging to each mass bin versus (the 
logarithm of) the total number of stars spanning all five mass bins (which 
provides a rough proxy for the total cluster mass). These lines have been 
weighted by adopting uncertainties for the number of stars in each mass bin 
using Poisson statistics. We have performed this comparison for all three 
circles (i.e. within one, two and three core radii). Our motivation for adopting 
this technique is as follows. If the fraction of stars belonging to each mass 
bin (relative to the total number of stars in all five mass bins), which we 
denote by f mi _ m2 , is constant for all cluster masses, then we would expect 
the number of stars in each mass bin N mi _ m2 to scale linearly with the total 
number of stars spanning all five mass bins N tot (since f mi _ m2 = N mi _ m2 /N tot ). 
However, if there is any systematic dependence of f mi _ m2 on the total cluster 
mass, then we should find that N mi _ m2 does not scale linearly with N to4 . In 
particular, the power-law index for a given mass bin should be sub-linear 
if the fraction of stars belonging to that mass bin systematically decreases 
with increasing cluster mass. Conversely, the power-law index for a given 
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mass bin should be super-linear if the fraction of stars belonging to that mass 
bin systematically increases with increasing cluster mass. The slopes and y- 
intercepts corresponding to the weighted lines of best-fit for each mass bin 
provide a means of directly quantifying the number of stars belonging to each 
mass bin as a function of the total cluster mass. In other words, our method 
provides a means of quantifying cluster-to-cluster differences in the stellar mass 
function function of the total cluster mass. 

We have chosen to count the number of stars belonging to each mass 
bin directly from the observations in order to quantify the dependence of the 
fraction of stars belonging to each mass bin on the total cluster mass (or, more 
specifically, the total number of stars spanning all mass bins). We note that an 
alternative, albeit considerably less precise, way to go about this would be to 
characterize the stellar mass functions using a log-normal function of the form 
given by Equation 16.21 Let this function be represented by \l/(m) = dN/dlnm. 
We can normalize the distribution of stellar masses over the entire mass range 
of interest (0.25-0.75 M ) by setting: 



where, in our case, m min = 0.25 M and m max = 0.75 M . The fraction of 
stars belonging to a given mass bin would then be given by: 



where mi and m.2 are the lower and upper mass limits, respectively, of the 
mass bin under consideration. 




(6.4) 




(6.5) 
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The reason why our technique provides a means of quantifying the uni- 
versality of the IMF, in addition to the effects had by two-body relaxation 
in modifying it, can be understood as follows. If we assume that all things 
are equal, such as tidal effects from the Galaxy and the degree of primor- 
dial mass segregation, two-body relaxation operates the fastest in low-mass 
clusters. All of the clusters in our sample are of comparably old age, which 
implies that the dynamical age increases with decreasing cluster mass. With 
increasing dynamical age, the stellar mass function should be more severely 
depleted of preferentially low-mass stars due to stellar evaporation induced 
by two-body relaxation. Based on this, if all of the clusters in our sample 
were born with similar IMFs, we would expect the slopes of their lines of 
best-fit to systematically increase with decreasing mass bin. This is because 
low-mass clusters should be more depleted of their low-mass stars. Therefore, 
the weighted lines of best-fit and, in particular their corresponding uncertain- 
ties, provide a means of quantifying the universality of the stellar IMF in the 
range of cluster masses and chemical compositions characteristic of Milky Way 
globular clusters. Large uncertainties caused by a large degree of scatter could 
be interpreted as evidence against a universal IMF. This is because we would 
only expect a systematic dependence of the number of stars belonging to each 
mass bin on the total cluster mass, and therefore small uncertainties for the 
slopes and y-intercepts corresponding to their lines of best-fit, if all clusters 
began with very similar IMFs in the first place. We note that although our 
method constrains the universality of the stellar IMF, it does not constrain 
its precise functional form. The reason for this is that we do not know how 
much dynamical evolution has actually occurred in the clusters in our sample. 
Therefore, we do not know how much the stellar mass function has been mod- 
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ified by two-body relaxation. However, as we will describe in Section 16.41 our 
observational analysis is ideally suited for comparison to theoretical models for 
globular cluster evolution, and this will allow us to constrain the exact shape 
of the IMF, in addition to the degree of primordial mass segregation. 

Finally, we note that mass segregation should also contribute to the 
systematic dependence on cluster mass we have described in the previous para- 
graph for the slopes and y-intercepts corresponding to the lines of best-fit for 
each mass bin. This effect should be the most severe for small circles (cen- 
tred on the cluster centre), and should become less important as increasingly 
larger circles are considered. This is because as we consider progressively larger 
fractional areas of the cluster, the number counts for each mass bin become 
less sensitive to the stars' spatial distributions throughout the cluster, and 
therefore less sensitive to the effects of mass segregation. 

6.3 Results 

In Figure 16.11 Figure 16.21 and Figure 16.31 we plot the number of stars in each 
mass bin versus the total number of stars spanning all five mass bins. Slopes (a) 
and y-intercepts (b) for the weighted lines of best-fit performed for each of these 
relations are shown in Table 16.11 along with their corresponding uncertainties. 
These were found using a bootstrap methodology in which we generated 1,000 
fake data sets by randomly sampling (with replacement) number counts from 
the observations. We obtained lines of best fit for each fake data set, fit a 
Gaussian to the subsequent distribution and extracted its standard deviation. 
We have found that our results are insensitive to both changes in bin width 
and bin centering. Specifically, our slopes and y-intercepts remain consistent to 
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within one standard deviation upon adjusting either the bin width or centering 
by up to roughly half a bin width. 
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Figure 6.1 The logarithm of the number of stars belonging to each mass bin as 
a function of the logarithm of the total number of stars spanning all five mass 
bins in the core. MSI corresponds to the mass range 0.65 - 0.75 M , MS2 to 
0.55 - 0.65 M , MS3 to 0.45 - 0.55 M , MS4 to 0.35 - 0.45 M , and MS5 to 
0.25 - 0.35 M . Lines of best fit are shown for each mass bin. 

As shown in Table 16.11 the slopes tend to systematically increase with 
decreasing mass bin. For the comparison within two core radii (shown in 
Figure 16721) . this trend applies to all mass bins. For the comparisons within 
the core and within three core radii (shown in Figure 16.11 and Figure 16.31 
respectively), however, only the highest three mass bins (MSI, MS2, and MS3) 
follow this trend of increasing slope with decreasing mass bin. After the third 
mass bin, the slopes for the lowest two mass bins remain about the same. Note, 
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Table 6.1 Lines of Best Fit for log N MS = (a ± Aa)log (N tof /10 3 ) + (b ± Ab) 





MSI (0.65-0.75 M Q ) 


MS2 (0.55-0.65 Mq) 


MS3 (0.45-0.55 M ) 


MS 4 (0.35-0.45 Mq) 


MS5 (0.25-0.35 M ) 


Circle 


a Aa b Ab 


a Aa b Ab 


a Aa b Ab 


a Aa b Ab 


a Aa b Ab 


< r c 

< 2r c 

< 3r e 


0.70 0.07 2.85 0.09 
0.81 0.08 2.74 0.12 
0.91 0.12 2.50 0.19 


0.86 0.05 2.56 0.05 
0.90 0.06 2.55 0.09 
1.04 0.11 2.30 0.17 


1.08 0.05 2.22 0.07 
0.99 0.03 2.35 0.05 
1.08 0.08 2.23 0.12 


1.08 0.09 2.17 0.12 
1.07 0.06 2.17 0.10 
1.02 0.05 2.31 0.10 


1.00 0.20 2.24 0.25 
1.19 0.10 1.91 0.18 
1.15 0.10 2.04 0.16 



Ph.D. Thesis — N. Leigh 



McMaster - Physics & Astronomy 











- 



- 



- 



- 


, 




— A* 


— 




0.5 


1 


1.5 




2 






— • — • 


, 




— >•» 






0.5 


1 


1.5 




2 

■ ■ — 




1 1 i i 1 i 1 i 


0.5 


1 

log N < Zr 

1 


1.5 




2 






^ , r " »~ 

, 










0.5 


1 

log N< B,, 

1 


1.5 




2 






p , 










0.5 


1 


1.5 




2 





Figure 6.2 The logarithm of the number of stars belonging to each mass bin 
as a function of the logarithm of the total number of stars spanning all five 
mass bins within two core radii from the cluster centre. The mass bins are the 
same as in Figure 16.11 

however, that the uncertainties on the slopes and y-intercepts are also higher 
for the lowest mass bins (MS4 and MS5). This is likely to be due to the fact 
that the photometric errors are the highest at these dim magnitudes, however 
they remain at most ~ 10% of the width in magnitude of their corresponding 
mass bins. The completeness corrections are also the largest for MS4 and MS5, 
and this also introduces additional uncertainty. Although we have taken the 
appropriate measures, it is important that these effects be properly quantified 
in order to reliably use our technique to constrain the degree of universality 
of the IMF. We will return to this in Section 16.41 
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Figure 6.3 The logarithm of the number of stars belonging to each mass bin 
as a function of the logarithm of the total number of stars spanning all five 
mass bins within three core radii from the cluster centre. The mass bins are 
the same as in Figure I6.ll 

On the other hand, when we consider only the largest three mass bins 
(MSI, MS2, and MS3) for which the photometric errors and level of incom- 
pleteness are the lowest, the difference in slopes and y- intercepts between ad- 
jacent mass bins can differ at better than the 2 — a confidence level. Moreover, 
if we compare non- adjacent mass bins, then this trend is typically significant 
at the 3 — cr confidence level. In order to improve upon these statistics, we 
have also calculated reduced chi-squared values with added intrinsic disper- 
sion for the relations for each mass bin. That is, for each mass bin we added 
a constant term to the uncertainty for each data point, found the value that 
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yielded a reduced chi-squared of one, and looked at the subsequent effects on 
the uncertainties for the line of best-fit. Based on this, we appear to be slightly 
over-estimating the uncertainties for the MSI, MS2 and MS3 mass bins using 
our bootstrap approach, and slightly under-estimating them for the MS4 and 
MS5 bins. This suggests that the slopes and y-intercepts differ at nearly the 
2 — a confidence level for the MSI, MS2 and MS3 mass bins, but are consistent 
to within one standard deviation for the MS4 and MS5 bins. We conclude that 
our results are the most reliable (at nearly the 2 — a confidence level) in the 
mass range 0.45-0.75 M . 

The change in the distribution of stellar masses as a function of the total 
cluster mass can be illustrated using pie charts, as shown in Figure 16.41 Using 
the slopes and y-intercepts provided in Table 16.11 for the comparison within 
two core radii, we have generated pie charts for three total numbers of stars 
(spanning all five mass bins), namely N toi = 10 3 , 10 4 , 10 5 . As is clear, low-mass 
stars are preferentially depleted, and this effect becomes increasingly severe 
with decreasing total cluster mass (or, equivalently, increasing dynamical age). 
From right to left, what the pie charts are showing are the mass functions of 
progressively dynamically older clusters. If this is indeed the cause of the 
observed depletion of low-mass stars in low-mass clusters, we are effectively 
looking at the evolution of the stellar mass function in time. 

6.4 Summary & Discussion 

In this chapter, we have obtained completeness-corrected stellar mass func- 
tions in the range 0.25-0.75 M for a sample of 33 globular clusters using 
data taken from the ACS Survey for Globular Clusters. We have presented a 
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Figure 6.4 Stellar mass functions depicted in pie chart form. The total area 
of each circle corresponds to the total number of stars spanning all five mass 
bins, and each pie slice shows the fraction of this total corresponding to each 
mass bin. Each of these fractions was calculated using the weighted lines of 
best-fit provided in Table IBTTl From left to right, the total number of stars used 
to generate each pie chart was 10 3 , 10 4 , and 10 5 . Dark blue corresponds to 
MSI, red to MS2, green to MS3, purple to MS4, and light blue to MS5. From 
right to left, the pie charts effectively show the mass functions of progressively 
dynamically older clusters. 



new technique to quantify cluster-to-cluster variations in the observed stellar 
mass functions. We have shown how our method can be used to quantify the 
universality of the IMF for a large sample of clusters and, when used in con- 
junction with theoretical models for globular cluster evolution, can be used to 
constrain the degree to which two-body relaxation has modified the currently 
observed stellar mass functions from their primordial form. To this end, we 
have obtained weighted lines of best-fit by comparing the number of stars in 
five different mass bins to the total number of stars spanning the entire mass 
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range within several circles centred on the cluster centre. 

As shown in Table I6.1[ our results suggest that the slopes for the lines 
of best-fit tend to systematically increase with decreasing mass bin. Assuming 
all of the clusters in our sample were born with similar IMFs, this is precisely 
what we expect from two-body relaxation. That is, we expect low-mass stars 
to become preferentially depleted via stellar evaporation, and we expect this 
effect to be the most severe in low-mass clusters. This trend is clearly seen 
in the observations, as demonstrated in Figure 16.41 Interestingly, the fraction 
of stars belonging to each mass bin is remarkably constant for all mass bins 
at the high-mass end of our sample. The total masses of these clusters are 
sufficientl y high t hat their half-mass relaxation times are comparable to their 



ages (e.g. 



Harris 



1996. 2010 update! ) . Therefore, we expect that the highest 
mass clusters in our sample should have been relatively unaffected by stellar 
evaporation due to two-body relaxation. It follows that their current mass 
functions should be relatively unchanged from their primordial values (ignoring 
the effects of stellar evolution). If true, our results could suggest that the 
slope IMF in Equation 16.11 was consistent with zero in the mass range 0.25- 
0.75 M for all of the clusters in our sample. This would be inconsistent with 
a "standard" Kroupa IMF. This hypothesis can be tested by comparing the 
results of the observational analysis presented in this chapter to simulations of 
globular cluster evolution. These can be used to directly quantify the effects 
of two-body relaxation on the expected slopes and y-intercepts for each mass 
bin given any set of initial conditions and IMFs. By comparing the results of 
these simulations to the results of our observational analysis, the set of initial 
conditions that best reproduces the observations can be found. This will tell 
us about the precise functional form of the IMF for the clusters in our sample, 
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in addition to the degree of primordial mass segregation. 

The uncertainties for the slopes and y-intercepts are sufficiently large 
that they are often consistent with those of their adjacent mass bins to within 
one standard deviation. However, when we consider only the largest three 
mass bins for which the photometric errors and level of incompleteness are the 
lowest, the difference in slopes and y-intercepts between adjacent mass bins 
can differ at better than the 2 — a confidence level. Moreover, the trend of 
increasing slope with decreasing mass is typically significant at the 3 — a confi- 
dence level if we compare non- adjacent mass bins. Based on this, we conclude 
that our results are consistent with a universal IMF for the clusters in our 
sample, and that the currently observed mass functions have primarily been 
modified by internal two-body relaxation. Having said that, it is important 
to note that a higher degree of primordial mass segregation effectively acts 
to in crease the rate of dy namical evolution and therefore stellar evaporation 



(e.g. 



Heggie fc Hut 



20031 ). Consequently, our results could also be consistent 
with a non-universal IMF that depends on the total cluster mass, provided the 
degree of primordial mass segregation also depends systematically on the total 
cluster mass. We also wish to point out here that the large uncertainties found 
for the slopes and y-intercepts, particularly for the two lowest mass bins, are 
primarily due to the relatively large photometric errors at these dim magni- 
tudes and incompleteness resulting from crowding. Despite the high quality 
of the data used in this study, these issues are currently unavoidable given the 
nature of the observations. This will be a key challenge for future studies to 
resolve, however the method we have presented in this chapter offers a robust 
means of performing future analyses. 

In terms of addressing the degree of universality of the IMF, our tech- 
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nique offers a considerable advantage over the standard power-law and log- 
normal forms used to characterize the stellar mass function shown in Equa- 
tion |6J] and Equation 16.21 This is the case when comparing the mass functions 
of a large sample of clusters. The reason for this is that our method segments 
the mass function into mass bins, which characterizes cluster-to-cluster dif- 
ferences in the stellar mass function for each mass bin individually ( 1.6. cLS cL 
function of stellar mass). This is done using the slopes and y- intercepts of the 
relations found by correlating the number of stars in each mass bin with the 
total number of stars spanning all mass bins. Conversely, standard Kroupa or 
Chabrier mass functions quantify the mass function much more globally via 
functional fits to the data over considerably larger mass ranges. Of course, 
the standard forms remain as useful as ever for characterizing the mass func- 
tions of individual clusters, whereas our method is not appropriate for this 



purpose. That bein g said, it is worth stating here that ou r resu 



sistent with those of 



De Marchi. Paresce &: Portegies Zwart 



ts are con- 



(120101 ). who fit a 



tapered power-law distribution function with an exponential truncation to the 
stellar mass functions of a sample of 30 clusters containing both young and 
old members. 

There are several possible sources of additional uncertainty that could 
have affected our analysis. For instance, we have not yet considered the role 
played by binaries. These tend to be unresolved in GCs, usually appearing 
as single objects located above the MS in the cluster CMD. Therefore, some 
of the objects in each mass bin are in fact binaries masquerading as single 
stars. Typical binary fractions in Milky Way G Cs are thought to be on the 



order of a 



Cool et al. 



ew t o a few tens o f a percent (e.g. 



2002 



Sollima et al. 



2008 



Rubenstein &: Bailyn 



Davis et al. 



1997 



20081 ). This suggests that 
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a non-negligible number of binaries could often be included in our number 
counts for each mass bin. However, there is no reason not to expect that 
binaries will contribute to each mass bin in roughly equal proportions. If the 
fraction of binaries is the same in all mass bins, then the presence of binaries 
should not have significantly affected the slopes. Although they will have 
affected the y-intercepts, this effect should not cause them to deviate from 
within the uncertainties and should not have affected the interpretation of our 
results. On the other hand, some observational evidence suggests that the 



binary fraction cou l d be inversely 



Sollima et al. 



2008 



Milone et al. 



proportional to the total c 



2008 



Knigge. Leigh fc Sills 



uster mass (e.g. 
. It is not 



clear how this might have affected our results since we do not know how each 
mass bin should be affected. It is our intent to address this issue in a future 
study once binary fractions become available for the majority of the clusters 
in our sample (Ata Sarajedini; private communication). 

Throughout our analysis, we have consistently compared two projected 
quantities. Therefore, effects related to projection should not have significantly 
affected our results. Moreover, these effects should become less severe upon 
considering progressively larger circles, and we have performed our analysis 
for circles with three different sizes. We note that this could perhaps help to 
explain why the scatter in the lowest mass bins appears to become reduced for 
the comparisons within two and three core radii (compared to the comparison 
within one core radius), although it is not clear why the agreement appears 
slightly better for the former comparison. 

Tidal effects from the Galaxy effectively a ct to reduce the time-scale 



for two-body relaxation (e.g. 



Heggie k, Hut 



20031 ) . The same effect can also 



be had by increasing the degree of primordial mass segregation (e.g. 



Spitzer 
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19871 ; iPortegies Zwart et all . 1200 ll ). Therefore, on average, we would expect 
clusters with smaller Galactocentric radii and higher initial concentrations to 
have experienced a higher degree of stellar evaporation. In an effort to quantify 
tidal e ffects from the Galaxy, we perf o rmed several cuts in perigalacticon dis- 



tance (IDinescu. Girard fc van Altena 



1999 



Casetti-Dinescu et al. 



20071 ) and 



re-performed our weighted lines of best-fit. Despite removing clusters from 



our sample with small perigalacticon distances for w hich it is typic a 



Heggie fc Hut 



ly ar - 



20031 ) 



gued that tidal effects should be the most severe (e.g. 
our slopes and y-intercepts, in addition to their uncertainties, remain more 
or less unchanged. This can be interpreted as rough evidence that tidal ef- 
fects from the Galaxy have not significantly affected our results. It is much 
less clear which clusters in our sample, if any, were born with high cen- 
tral concentrations. Based on current observations, there is no known rea- 
son not to expect a universal degree of primordial mass segregation (e.g. 



Portegies Zwart. McMillan fc Gieles 



20101 ). If this was the case, then our 



results could be interpreted as indicative of a universal IMF for Milky Way 
globular clusters. Unfortunately, obs ervational constraints for primordial bi- 
nary populations are also lacking (e.g. 



McKee fc Ostriker 



20071 ) . Binaries play 



an important role in star cluster evolution by, for instance, providing an energy 
source t hat serves to resist a cluster's tendency toward increasing its central 



density ( iHut 



19831 ) . They are therefore an important consideration in decid- 
ing the dynamical evolution of the stellar mass function, and their role will 
need to be addressed in future studies. 

It is interesting and even surprising that, despite all of the aforemen- 
tioned factors expected to affect the dynamical evolution of globular clusters, 
we observe relatively little scatter in the relations within two and three core 
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radii (Figure l6\2l and Figure [673]) . Moreover, the slopes systematically increase 
with decreasing mass bin. Based on this, our results are consistent with a sin- 
gle universal IMF for all of the clusters in our sample that has been modified 
via internal two-body relaxation by an amount determined by the total cluster 
mass. This is a new result given the old ages and therefore low metallicities 
([Fe/H] ~ -2.28 - (-0.37)) of the clusters that comprise our sample. However, 
the exact form of the IMF required to reproduce the current observations is 
still not clear, nor is the role of primordial mass segregation. To what degree 
have low-mass clusters been depleted of their low-mass stars? What combina- 
tion of IMFs, Galactocentric radii, initial concentrations and primordial binary 
fractions should evolve dynamically to best reproduce the currently observed 
mass functions? These questions can only be answered using simulations of 
globular cluster evolution that explore a range of initial conditions and IMFs. 
Our observational analysis of the ACS data is ideally suited for comparison to 
these types of models. This would allow us to constrain the precise shape of 
the IMF in the mass and metallicity range of interest, as well as to learn about 
primordial mass segregation. This will be the focus of a forthcoming study. 
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Chapter 



7 



Summary & Future Work 



Within the last few decades, theoretical work has painted a comprehensive 
picture for the various forms of gravitational interactions operating within 
star clusters. And yet, direct observational confirmation that many of these 
processes are actually occurring is still lacking. The results presented in this 
thesis have connected several of these processes to real observations of star 
clusters, in many cases for the first time. This has allowed us to directly link 
the observed properties of several peculiar stellar populations to the physical 
processes responsible for their origins. Given the importance of star clusters 
for star formation, this also represents a key step toward re-constructing the 
history of our Galaxy. 

In Chapter [2j we have presented a new adaptation of the mean free 
path approximation. Our technique provides analytic time-scales for the rate 
of close gravitational encounters between single, binary and triple stars in 
dense star systems. With it, we showed that encounters involving triple stars 
occur more commonly than any other type of dynamical interaction in at 
least some star clusters, and that these could be the dominant dynamical pro- 
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duction mechanism for stellar mergers. This is a new result with important 
implications for both star cluster evolution and the formation of several types 
of stellar exotica. Our method can be generalized for application to systems 
composed of several different types of particles that evolve under the influence 
of any force-mediated interactions. For example, our technique is well suited 



for application to co 
medium (ISM) (e.g. 



lisions between a toms and dust grains in the interstellar 



Spitzer . 



1941a 



19421 ). In this case, the relevant forces 



governing the dynamics are electromagnetic instead of gravitational, however 
the general form adopted for the mean free path approximation remains the 
same. Specifically, we must only replace the gravitationally-focused collisional 
cross-sections with their charge-focused analogues. Recent observations have 
provided estimates for the concentrations of the different constituents of the 
ISM, in addition to their temperatures (e.g 



del Burgo et al. 



2003 



Kiss et al. 



20081 ). These data provide the required ingredients to obtain analytic esti- 
mates for the rates of collisions between atoms, molecules, as well as both 
small and large dust grains in the ISM, and would allow us to make predic- 
tions for the evolution of their relative concentrations. This is an important 
step in the development of our understanding of dust grain coagulation and, by 
extension, the initial pha ses of star formation in dense molecular clouds (e.g. 



McKee fc Ostriker 



surrounding IS M (e.g. 



tary disks (e.g. 



20071). the interactio n of winds from evolved stars with the 



Glassgo 



Absil fc Mawet 



2010 . 



19961 ). and planet formation in protoplane- 



In Chapter [31 we introduced a statistical technique to compare the rel- 
ative sizes of different populations of stars in a large sample of star clusters. 
We refined this technique in Chapter HI and applied it to a large sample of 
clusters collected using state-of-the-art observations. Our results suggest that 
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dynamical effects do not significantly affect the relative sizes of the different 
stellar populations. This is the case for all stellar populations above the main- 
sequence turn-off in the cluster CMD. Blue stragglers alone present a possible 
exception to this, since we have identified a general trend in which more mas- 
sive clusters are home to proportionately smaller blue straggler populations. 
This provides compelling evidence in favour of a binary origin for blue strag- 



glers in even the dense cores of massive star cl usters where direct co 



lisions 



2009|). 



between stars are expected to occur frequently ( iKnigge. Leigh fc Sills 
Although we have applied these techniques to observations of dense star clus- 
ters throughout this thesis, they can be generalized for application to a number 
of other studies related to population statistics. For example, massive elliptical 
galaxies have been identified in Galaxy Clusters and Groups, and these could 
have formed from the mergers of smaller galaxies. Observational studies have 
now been performed that are allowing for the compilation of catalogues provid- 
ing population statistics for sever al different galaxy types in a large number o f 



Galaxy Clusters and Groups (e.g. 



Abell 



1958 



Abell. Corwin fc Olowinl . 



Our technique is well suited for application to these data, and would allow for 
a systematic comparison between the relative population sizes of the different 
galaxy types. Among other things, this would help to constrain the origins of 
massive elliptical galaxies in Galaxy Clusters in much the same way we have 
constrained the origins of BSs in globular clusters. 

In Chapter [5l we present an analytic model for blue straggler forma- 
tion in globular clusters. Our model considers the production of blue stragglers 
throughout the entire cluster, and tracks their accumulation in the core due 
to dynamical friction (or, equivalently, two-body relaxation). Our results sup- 
port the conclusion that blue stragglers are descended from binary stars, as 
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first reported in 



Knigge. Leigh fc Sills! ( 2009 ) using the technique presented 



in Chapter [3] and Chapter HI Our model is applicable to studying the radial 
distributions of other stellar and binary populations in globular clusters. This 
is a sorely needed theoretical tool that can be used to address a number of 
recent observational studies that reported peculiarities among t he radial dis 



tributions of seyeral different stellar and binary populations (e.g . 



Rood 



1973 



Fusi Pecci et al 



1993 



Ferraro et al. 



200 1 



Lanzoni et al. 



20071). 



Finally, in Chapter El we extend the statistical technique introduced in 
Chapter [3] and Chapter [4] to include the entire main-sequence. By applying our 
method to the ACS data, we have obtained the first quantitative constraints for 
the degree of universality of the stellar initial mass function for a large sample 
of star clusters spanning a wider range of masses and chemical compositions 
than ever before considered. Given the very old ages of the clusters in our 
sample, our results have important implications for our understanding of star 
formation in the early Universe. Our results are consistent with a remarkably 
universal IMF in old massive star clusters, and are well suited for comparison 
to theoretical simulations of globular cluster evolution. This would provide 
a simple means of directly quantifying the extent to which the stellar mass 
function has been modified from its primordial form by two-body relaxation. 
In future studies, this will allow us to obtain the very first constraints for 
the precise functional form of the IMF and the degree of primordial mass 
segregation in a large sample of old star clusters. 

The results presented in this thesis have significantly advanced our 
understanding of stellar dynamics, stellar evolution and, in particular, the in- 
terplay that occurs between the two in dense star clusters. But we are not 
yet done. Our results have raised several important questions related to these 
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topics and, by extension, the history of our Galaxy. Once again, we are left 
asking: Where do we go from here? For example, how else might we use 
observations of blue straggler populations to learn about the dynamical histo- 
ries of their host clusters? How can we use this information to constrain the 
interactions that occur in clusters between stellar dynamics and binary star 
evolution? How do these interactions affect the observed properties of their 
binary populations and the various types of stellar exotica they are home to? 
There is also the issue of the universality of the stellar initial mass function. 
What combination of initial concentrations and mass functions should have 
evolved dynamically over the lifetimes of star clusters to reproduce the cur- 
rently observed mass functions? What can this tell us about primordial mass 
segregation in massive star clusters, and the exact form of the stellar IMF? 
All of these issues relate to the over-arching questions: How do stars form?; 
and: How did the history of our Galaxy unfold? The next few years promise 
to hold exciting advances in the search for answers to these questions. 
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Appendix 



A 



Appendix A 



In this appendix, we present the collisional cross sections and time-scales used 



in Chapter [2j The gravitationally-focused cross sections for 1+1, 1 
+3, 2+3 and 3+3 collisions can be found using Equation 6 from 



2- 2+2. 



Leonard 



(119891 ). Neglecting the first term and assuming that binary and triple stars 
are on average twice and three times as massive as single stars, respectively, 
this gives for the various collisional cross sections: 

8ivGmR 

o-i+i ~ 2 ' ( A - X ) 

V rel 

MGma b 

Oi +2 ~ 2 , ( A -2) 

rel 

MGma b 

<T 2 +2 ~ 2 ' 1 A "JJ 

V rel 

4iiGma t 

~ 2 ' ( A - 4 ) 

V rel 

5irGm(a b + a t ) , . 

02+3 2 ' ( A - 5 ) 

V rel 

12nGma t 

03+3 ~ 2 • ( A - 6 ) 

V rel 
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Table A.l Pericenters Assumed for Each Encounter Type 

Encounter Type Pericenter 



1+1 2R 

1+2 a b /2 
1+3 ot/2 
2+2 a b 

2+3 (a b + a t )/2 
3+3 a t 



Values for the pericenters assumed for the various types of encounters are 
shown in Table IA.1| where R is the average stellar radius, a b is the average 
binary semi-major axis and at is the average semi-major axis of the outer orbits 
of triples. 

In general, the time between each of the different encounter types can 
be found using Equation 12.101 and the gravitationall y- focused cross sections 



Leonard! ( 119891 ). we can 



for collision given above. Following the derivation of 
write the encounter rate in the general form: 



r ^ x+y -^x^y^x+y^x+yy v^"'^) 

where N x and Tl y £1X6 the number and number density, respectively, of single, 
binary or triple stars and v x+y is the relative velocity at infinity between objects 
x and y. For instance, the time between binary-binary encounters in the core 
of a cluster is given by: 

r 2+2 = L3xl07r 2 (^) 3 ( i0V " V 

/ Vrms \ f 0-5M Q \ ( lA U \ 

\5km/sJ\<m>J\ a b ) 



n ° ' (A.8) 
years. 



Similarly, the times between 1+1, 1+2, 1+3, 2+3 and 3+3 encounters are 
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riven by: 

,/lpC\3 /10 3 pc" 3 \ 2 



and 



r 1+1 = l.lxl0 10 (l-/ & -/ 4 )- 2 (^) ( 



/ v rms \ /O.5M \ /0.5R 



n 
years, 



\5km/sJ \< m > / \ < R > 

nH = !Ux io'(i- A -/,)- 1 A- , (^)'(^ 

V r c / V n 
/ Vs \ / O.5M \ ( }AU \ earg 
\5km/sJ \< m >J V a& / 



n„ = 2.6 xio'(i- A- /,)-/,-' (^)\^f K2 

f Vrms \/ 0.5M Q \/ lAU \ er 

\5km/sJ V< m >J V a t / 

r 2+3 = 2.0xl0W(^) 3 (^ 

V r r / V ran 



0.5M m \ / lAf/ 



/ Vrms \ / U.OMq \ / 

\5km/sJ \< m >J Ka^ + a t 



" 3 \ 2 

years, 



r 3+3 = 8^3xl0Vr 2 (^) 3 ( 1 ^) 2 

/ Vrms \ / 0-5M Q \ / 1A?7 \ 

\5km/sJ\<m>/\ a t J 



(A.9) 



(A.10) 



(A.ll) 



(A.12) 



(A.13) 
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Appendix 



B 



Appendix B 



In this appendix, we present our selection criteria for BS, RGB, HB and 
M STO stars used in Chapter HI Our method is similar to that described 



in 



Leigh. Sills fc Knigge 



( 120071 ). and we have used this as a basis for the se- 
lection criteria presented in this chapter. First, we define a location for the 
MSTO in the (F606W-F814W)-F814W plane using our isochrone fits. The 
MSTO is chosen to be the bluest point along the MS of each isochrone, which 
we denote by ((V-T)msto^msto)- 111 order to distinguish BSs from MSTO 
stars, we impose the conditions: 

F8UW < m 1 (F606W - F8UW) + b u , (B.l) 
where the slope of this line is mi = —9 and its y-intercept is given by: 

bn = {hisTO ~ 0.10) - m x {{V - I) MSTO - 0.10) (B.2) 

Similarly, we distinguish BSs from HB stars by defining the following 
additional boundaries: 

FSUW > m^FWGW - F8UW) + b 12 (B.3) 
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F8UW > m 2 {F606W - F8UW) + b 21 


(B.4) 




F8UW < m 2 {F606W - F8UW) + b 22 


(B.5) 




F814H/ > m#£(F606Vv — F814VK) + 


(B.6) 


(F606W - 


F814W) > (V-/W 


(B.7) 




F8UW < Imsto, 


(B.8) 


where m 2 = 6, iuhb 


= -1.5 and (V - /) flfl = (V - I) MS to - 0.4. 


We also 


define: 






b\2 = 


{Imsto - 0.55) - mi((V - I)msto - 0.55) 


(B.9) 


621 = 


(Wo - 0.80) - m 2 ({V - I)msto + 0.10) 


(B.10) 


b 22 = 


{Imsto + 0.30) - m 2 {{V - I)msto - 0.20), 


(B.H) 



and bus = Ihb + 1-2, where Ihb roughly corresponds to the mid-point of 
points that populate the HB and is chosen by eye for each cluster so that our 
selection criteria best fits the HB in all of the CMDs in our sample. 

We apply a similar set of conditions to the RGB in order to select stars 
belonging to this stellar population. These boundary conditions are: 

F8UW > m HB (F606W - F8UW) + b HB (B.12) 

F8UW > m RGB (F606W - F8UW) + b 31 (B.13) 

F8UW < m RGB {F606W - F8UW) + b 32 (B.14) 

F8UW < I RGB , (B.15) 

where m R GB = —23, Irgb is defined as the F814W magnitude corresponding 
to a core helium mass of 0.08 M Q and: 

&31 = {Imsto — 0.60) — m,RGB{{V — I)msto + 0.05) (B.16) 
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b 3 2 = (Wo - 0.60) - m RGB {(V -I)msto + 0.25) (B.17) 

Core helium-burning stars, which we refer to as HB stars, are selected 
if they satisfy one of the following sets of criteria: 

F8UW > 
F8UW < 
(FQ06W-F8UW) < 

F8UW 
F8UW 
(F606W - F8UW) 

or 

F8UW 
F8UW 
(F606W - F8UW) 

We define (V — I)hb on a cluster-by-cluster basis in order to ensure that we 
do not over- or under-count the number of HB stars. This is because the 
precise value of (F606W-F814W) at which the HB becomes the RGB varies 
from cluster-to-cluster. In addition, the precise location of the transition in 
the cluster CMD between HB and EHB stars remains poorly understood. To 
avoid this ambiguity, we consider HB and EHB stars together throughout our 
analysis, and collectively refer to all core helium-burning stars as HB stars 
throughout this chapter. 



m HB (F606W - F8UW) + {b HB - 1.0)(B.18) 
m HB (F606W - F8UW) + b HB (B.19) 
(V - I)msto + (V- I) HB , (B.20) 

> m HB (F606W - F8UW) + b HB (B.21) 

< Wo + 2.5 (B.22) 

< (V-I)msto-OA, (B.23) 



< mi (F606W - F8UW) + b 12 (B.24) 

> m HB (F606W - F8UW) + b HB (B.25) 

> (V — I)msto ~ 0.4 (B.26) 
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Finally, MSTO stars are selected according to the following criteria: 

F8UW > I RGB (B.27) 
FZIAW > mi(F606W-F814W) + &n (B.28) 
F8UW < (V-I)msto (B.29) 
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